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This  manual  is  written  to  give  the  user  of  the  DRA 
model  the  specific  information  he  will  need  to: 

a)  set  up  the  input  file 

b)  run  the  program 

c)  interpret  the  results 


This  model  describes  the  interaction  between  the  driver  of  a 
vehicle  and  an  air  cushion  restraint  system.  The  air  cushion 
is  mounted  to  a steeri ng  wheel/s teering  column  assembly  which 
is  modeled  to  stroke  and  rotate  during  the  crash.  A previous 
model  called  DRAC  could  not  model  this  rotation;  hence  the 
descriptor  "R"  for  "rotation"  is  added  to  this  model  to  dist- 
inguish it  from  the  functionally  similar  but  less  powerful  DRAC 
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1.0  INTRODUCTION 


The  purpose  of  this  manual  is  to  give  to  the  potential  DRACR  user 
specific  information  to: 

a.  Set  up  the  input  file 

b.  Run  the  program 

c.  interpret  the  results 

Prior  to  discussing  these  things,  however,  some  background  information 
on  the  DRACR  program  will  be  presented. 

DRACR  is  an  acronym  for  "DRiver  Air  Cushion  - Rotation".  As  the 
title  indicates,  the  program  was  written  to  describe  the  interaction, 
during  a crash,  between  the  driver  of  a vehicle  and  his  airbag.  The 
airbag  is  mounted  to  a steering  wheel  and  steering  column  which  is  mod- 
eled to  not  only  stroke  but  also  rotate  during  the  crash.  A previous 
computer  program  (DRAC),  also  written  by  Fitzpatrick  Engineering,  does 
not  model  the  rotation  of  the  steering  wheel  and  column;  Hence,  the 
term  "rotation"  is  added  to  the  descriptor  of  the  DRACR  program  so  it 
can  be  distinguished  from  the  functionally  similar  but  less  powerful 
DRAC  program. 

Other  programs  have  been  written  to  describe  the  interaction  be- 
tween a driver  and  an  airbag  but  none  were  specifically  suitable  to 
the  characteristics  exhibited  by  certain  types  of  restraint  systems. 
Some  of  the  unique  characteristics  of  the  DRACR  program  which  exist  in 
combination  to  provide  the  greatest  versatility  and  accuracy  possible 
for  a hardware  design  type  of  computer  program  are: 
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1.  It  has  the  capability  of  simulating  an  airbag  shape  typical 
of  the  ellipsoidal  shape  which  almost  all  driver  bags  have. 

2.  The  program  is  able  to  independently  describe  the  simultan- 
eous effects  of  steering  column  crush,  steering  wheel  crush, 
steering  column  rotation  with  respect  to  the  compartment,  and 
steering  wheel  rotation  with  respect  to  the  column  and  then 
assess  their  effect  on  the  degree  of  injury  imparted  to  the 
driver . 

3.  The  program  is  set  up  to  be  a design  tool  oriented  to  the 
user  requirements  of  a typical  restraint  system  engineer  with 
both  the  formulation  and  the  input-output  in  units  commonly  used 
and  measured  in  normal  hardware  design  and  test  environments. 

4.  As  a design  tool,  the  program  is  oriented  toward  the  test  hard- 
ware actually  encountered  in  most  test  situations.  For  example, 
past  computer  programs  might  model  the  driver  very  accurately 
but  neglect  the  bag  shape  actually  used  and/or  the  column  bind- 
ing and  frictional  forces  which  are  almost  always  present  and  in- 
fluence the  results  greatly. 

5.  DRACR  is  derived  to  have  the  same  level  of  detail,  complexity 
and  accuracy  for  all  of  the  restraint  system  components.  For 
example,  there  are  several  different  types  of  steering  columns, 
all  of  which  behave  somewhat  differently.  It  does  little  good 
to  have  an  elaborate  airbag  or  driver  algorithm  and  then  to  in- 
clude it  with  a steering  column  algorithm  that  is  so  simple  it 
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is  described  by  only  a force-stroke  characteristic . To  be 
accurate  and  capable  of  describing  the  important  differences 
between  column  types  in  the  areas  of  frictional  and  binding 
effects  and  the  wheel  and  column  resistance  to  rotation,  a 
program  such  a DRACR  is  required. 


-3- 


2.0  PROGRAM  DESCRIPTION 


DRACR  is  a two-dimensional,  lumped  mass  computer  model  of  the 
driver  interacting  with  an  ellipsoidal  airbag  mounted  to  a steering 
column  and  wheel  one  or  both  of  which  may  stroke  and/or  rotate  de- 
pending on  the  relative  stiffness  and  strength  of  the  individual 
components.  The  model  includes  five  masses.  Three  represent  the 
driver  (i.e.,  the  head,  torso  and  lower  body)  while  two  describe  the 
restraint  system  (i.e.,  the  steering  wheel  and  the  steering  column). 

The  airbag  is  simulated  by  an  oblate  spheroid  into  which  a program- 
med amount  of  gas  flows.  By  adjusting  the  airbag  vent  size,  a sel- 
ected amount  of  gas  can  be  exhausted  during  the  crash  in  order  to 
attenuate  peak  head  and  chest  g's  and  the  rebound  velocity. 

2 . 1 PROGRAM  FORMULATION 

The  mathematical  formulation  of  the  equations  of  motion  follows  the 
classical  Lagrangian  derivation  (see  Appendix  A of  this  manual)  with 
body  pivot  points  at  A and  B of  Figure  1.  The  lower  body  mass  (hips 
and  legs)  is  constrained  to  move  horizontal ly . 

DRACR  uses  a fixed  time  step  integration  routine  to  solve  the 
differential  equations  of  motion  numerically.  The  integration  routine 
chosen  was  the  Adams-Moul ton  predictor  corrector  method,  with  the  fourth 
order  Runge-Kutta  method  employed  to  determine  the  first  four  solution 
points.  DRACR  has  been  written  in  FORTRAN  IV  and  was  developed  to  oper- 
ate in  an  interactive  time  share  mode  for  greatest  versatility  for  the 
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Figure  1.  Schematic  of  "DRACR"  Computer  Model. 


design  engineer.  The  program  is  self  contained,  in  that  no  external 
routines  are  required  for  execution.  It  is  also  modular  in  construct- 
ion, so  as  to  facilitate  the  addition  of  other  subroutines  at  a future 
time,  if  desired.  The  data  input  is  from  a previously  created  disk 
file  and  the  input  parameters  appear  directly  on  the  terminal,  immed- 
iately preceding  the  complete  program  output. 

As  mentioned  above,  DRACR  has  been  programmed  in  modular  format 
with  several  subroutines.  The  reason  for  this  is  to  enable  other  fea- 
tures to  be  easily  added  at  a future  time.  For  example,  consider  the 
tabular  data  input.  For  simple  data  where  the  particular  value  of  the 
dependent  variable  is  a function  only  of  the  value  of  the  independent 
variable,  a simple  table  look-up  and  interpolation  subroutine,  "LOOK- 
UP" is  provided.  Gas  flow  versus  time,  vehicle  acceleration  versus 
time,  neck  torque  versus  angle,  column  force  versus  crush,  and  steer- 
ing wheel  force  versus  crush  are  examples  of  this  method  of  data  retrie- 
val . 

However,  in  those  cases  in  which  the  dependent  variable  is  a funct- 
ion not  only  of  the  value  of  the  independant  variable  but  also  de- 
pends on  whether  the  independent  variable  is  increasing  or  decreas- 
ing, a different  subroutine,  "SPRING",  that  allows  for  plastic  behavior, 
is  used.  This  subroutine  is  used  principally  for  those  cases  in  which 
hysterises  (or  plastic  action  of  a deformable  member)  is  modeled.  In 
these  cases,  one  must  not  only  specify  the  values  for  the  dependent 
variable  for  different  values  of  the  independent  variable,  but  must  al- 
so specify  the  "unload  slope(s)"  for  those  conditions  in  which  the  member 
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is  undergoing  unloading  during  a lessening  of  the  degree  of  deformation. 
Knee  restraint  force  versus  crush,  and  seat  friction  force  versus  stroke 
are  handled  by  this  subroutine. 

2.2  AIRBAG 

Most  of  the  dri ver/airbag  interaction  models  being  used  today  rely 
on  relatively  simple  spherical  or  cylindrical  airbag  shapes  (in  which  the 
bags  exhibit  a constant  radius  of  curvature  regardless  of  impact  angle). 
Unfortunately,  these  simple  bag  shapes  do  not  adequately  describe  the 
shape  of  most  of  the  airbags  presently  being  used.  Of  all  the  geometric 
shapes  that  could  be  postulated  as  candidates  for  the  driver  airbag,  the 
ellipsoid  (oblate  spheroid  to  be  specific)  is  most  nearly  the  shape  of 
the  driver  airbag.  For  this  reason  we  chose  the  oblate  spheroid  as  the 
shape  upon  which  to  base  the  bag  shape  algorithm  (see  Appendix  B of  this 
manual).  Unfortunately,  the  oblate  spheroid  is  not  as  mathematically 
easy  to  describe  as  the  sphere  or  cylinder,  since  the  bag  radius  of  cur- 
vature and  the  intercepted  volume  of  the  airbag  are  not  independent  of 
the  angle  of  torso  inclination. 

Figure  2 shows  a simple  schematic  of  the  airbag,  along  with  the 
variables  necessary  to  describe  the  airbag  shape.  The  airbag  is  assumed 
to  be  symmetric  about  a line  coincident  with  the  steering  column's  long- 
itudinal axis. 

2 . 3 STEERING  COLUMN 

As  previously  mentioned,  the  proper  modeling  of  the  steering  col- 
umn is  very  important  to  the  overall  accuracy  of  the  program  in  repro- 
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Figure  2.  DRACR  Airbag  Geometry. 
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ducing  "real  world"  behavior  for  the  driver  restraint  system.  Most 
driver  restraint  system  models  specify  only  a force  versus  stroke 
characteristic  for  the  column.  In  many  cases  these  models  neglect 
column  angle,  column  mass,  column  frictional  and  binding  effects  and 
the  location  and  strength  of  the  specific  bracketry  that  supports  the 
column.  In  setting  up  the  column  for  this  program,  we  have  chosen  a 
generic  type  of  column  that  is  widely  used.  This  column  is  the  General 
Motors  "bal 1 -and-groove"  type  column  shown  in  Figure  3.  Other  types 
could  have  been  specified,  but,  because  of  the  modular  program  construc- 
tion, other  types  may  be  easily  added  as  required.  As  it  turns  out,  the 
GM  type  column  is  fairly  typical  of  a wide  variety  of  columns,  so  the 
program  can  be  used  "as  is"  for  a wide  variety  of  situations  by  properly 
adjusting  the  input  parameters  that  describe  this  steering  column. 

The  program  is  set  up  to  calculate  the  airbag  forces  and  pressure 
before  moving  on  to  the  column  force  calculations.  Therefore,  the  air- 
bag forces  are  resolved  into  a column  axial  force,  a column  normal  force 
and  a column  bending  moment  as  applied  loads  to  the  column.  This  is  done 
before  entering  the  column  stroke  and  rotational  dynamics  portion  of  the 
routine.  All  the  user  need  specify  is  the  column  mass,  the  column  angle, 
the  basic  force-stroke  properties  of  the  column,  the  coefficient  of 
friction  at  the  column's  internal  and  external  "rub"  surfaces,  the  perti- 
nent column  dimensions,  and  the  location  and  strength  of  the  column  sup- 
port bracketry  for  the  computer  to  calculate  the  complete  column  dynamics. 

One  might  wonder  why  it  is  so  important  to  accurately  describe  the 
steering  column  rotation.  The  main  reason  is  that  it  commonly  occurs,  and 
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Figure  3.  Geometry  of  DRACR  Steering  Column. 


when  it  does  it  almostuniversal ly  causes  higher  injury  measures  than 
would  occur  otherwise.  The  reason  for  the  higher  injury  measures  is 
that  when  the  brackety  collapse  or  yielding  causes  column  rotation, 
the  rotation  is  almost  always  in  an  upward  direction  so  that  bag 
placement  moves  higher  and  higher  on  the  driver  resulting  in  less  and 
less  effective  torso  mass  being  brought  to  bear  on  the  column.  When 
this  happens,  the  column  ceases  to  stroke,  chest  g's  increase  and  the 
driver  begins  to  submarine.  Using  DRACR  we  may,  therefore,  determine 
whether  bracket  strengths  are  adequate  to  prevent,  or  at  least,  control 
such  an  occurance.  In  a slightly  different  vein,  we  could  determine 
just  what  loads  the  column  brackets  must  be  designed  to  withstand  with- 
out catastrophic  collapse. 

Appendix  C contains  additional  details  on  the  steering  column  algo 
rithm  derivation. 

2.4  STEERING  WHEEL 

The  steering  wheel  as  modeled  by  DRACR  has  two  degrees  of  freedom. 
It  may  crush  axially  and/or  rotate  up  or  down  relative  to  the  column. 
Like  the  steering  column,  the  accurate  simulation  of  the  steering  wheel 
is  important  to  achieving  overall  simulation  realism. 

In  terms  of  wheel  crush,  recent  testing  has  shown  that  typical  air 
bag  systems  can  be  made  safer  and  less  complex  by  using  advantages  in- 
herent in  a crushable  steering  wheel.  The  overall  system  can  be  made 
safer  since  wheel  crush  adds  to  the  total  controlled  stroke  in  the  pass 
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enger  compartment  plus  the  steering  wheel,  as  opposed  to  the  steer- 
ing column,  is  of  relatively  low  mass  with  no  friction  resulting  in 
generally  lower  injury  measures  and  greater  reliability.  Additionally, 
wheel  crush  can  add  redundancy  to  the  total  system  by  adding  a second 
mechanical  method  by  which  energy  may  be  absorbed  in  case  the  column 
fails  to  stroke.  Given  a choice  between  designing  a mechanical  energy 
absorber  into  a steering  column  or  into  a crushable  steering  wheel 
(assuming  the  wheel  may  crush  as  much  as  the  column  under  consideration), 
it  is  much  less  complex,  less  costly  and  more  reliable  to  implement 
this  crush  into  the  steering  wheel  as  opposed  to  the  steering  column. 

Now,  in  terms  of  steering  wheel  rotation,  recent  testing  has  also 
shown  that  a "self  aligning"  steering  wheel  has  several  desirable  attri- 
butes, some  of  which  are  listed  below.  Thus,  the  self  aligning  (rotat- 
ing) wheel  results  in: 

a)  Lower  contact  forces  if  the  bag  "bottoms"  on  the  wheel  rim. 

b)  Better  bag  placement  on  the  chest  since  the  wheel  has  this 
tendency  to  align  itself  with  the  chest. 

c)  More  total  energy  absorbing  stroke. 

Advantages  of  being  able  to  simulate  the  steering  wheel  rotation 
with  DRACR  are: 

a)  A more  realistic  simulation  of  what  actually  happens  to  the 
steering  wheel  in  an  actual  crash. 

b)  Makes  it  possible  to  investigate  the  relative  importance  of 
factors  that  cause  steering  wheel  rotation  and  which  factors 
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cause  rotation  in  either  the  upward  or  downward  directions 
so  as  to  obtain  the  desired  effect  in  the  final  design. 

c)  Makes  it  possible  to  tailor  the  steering  wheel  rotation- 
al stiffness  to  obtain  lowest  injury  levels. 

In  order  to  obtain  an  accurate  simulation  of  steering  wheel  be- 
havior in  the  crash,  the  user  of  DRACR  must  specify  the  steering  wheel 
mass,  the  mass  of  the  airbag  and  inflator  package,  the  wheel  angle  from 
a vertical  reference  (normally  this  angle  is  the  same  as  the  column 
angle  initially,  but  can  be  different  if  a vehicle  with  a tilt  wheel 
is  being  simulated),  the  basic  force-crush  properties  of  the  wheel,  the 
pertinent  wheel  dimensions,  and  the  torsional  resistance  at  the  wheel/ 
column  interface  in  both  the  upward  and  downward  directions. 

Figure  4 shows  a schematic  of  the  original  and  final  positions  of 
the  column  and  wheel  for  a typical  crash  situation.  In  the  figure  all 
four  mechanisms  (wheel  crush,  wheel  rotation,  column  crush  and  column 
rotation)  are  evident. 

One  might  wonder  how  the  DRACR  program  is  set  up  to  determine  which 
of  the  four  modes  of  energy  absorption  will  actually  occur  in  a given 
situation.  The  method  is  explained  in  the  following. 

All  four  modes  are  assumed  possible  each  time  the  program  computes 
new  values  each  time  through  the  computation  loop.  For  a one  millisecond 
integration  interval  this  would  be  every  one  millisecond.  The  com- 

puter is  instructed  to  compare  the  applied  torque  or  force  with  the  re- 
sisting torque  or  force  for  each  of  the  four  possible  yield  modes.  For 
each  case  in  which  the  applied  torque  or  force  is  less  than  the  resisting 
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Figure  4.  Schematic  of  "DRACR"  Steering 
Column  and  Wheel  Model. 


torque  or  force,  the  computer  allows  no  further  crush  or  rotation  due 
to  this  applied  torque  or  force  during  this  time  interval. 

In  those  cases  where  more  than  one  yield  mode  is  possible  during 
the  time  interval,  the  one  that  is  actually  allowed  to  occur  is  deter- 
mined by  the  "minimum  energy"  principle.  That  is,  the  energy  absorption 
mode  that  requires  the  1 east  expenditure  of  energy  per  hori zontal  distance 
traveled  by  the  deformi nq  i tern  during  the  time  i nterval  is  the  mode  that 
is  allowed  to  occur.  Other  modes  are  suppressed  until  the  next  time 
through  the  loop  at  which  time  the  whole  process  begins  again. 

2 . 5 KNEE  RESTRAINT  AND  SEAT  FRICTION 

The  DRACR  computer  program  is  set  up  to  accept  tabular  input  for 
the  force-crush  properties  of  the  knee  restraint  and  the  force-displace- 
ment properties  of  seat  friction.  It  is  the  sum  of  these  two  forces 
which  primarily  determine  what  the  lower  body  restraint  for  the  lower 
torso,  hips  and  legs  will  be. 

The  user  specifies,  in  tabular  format,  the  values  these  forces 
and  displacements  are  to  take.  In  addition,  the  user  specifies  the 
"unload  slope",  so  the  program  can  compute  the  unloading  force  path  to 
be  taken  during  rebound  away  from  the  knee  bolster  or  movement  rear- 
ward across  the  seat. 

The  user  also  specifies  the  initial  angles  of  the  tibia  and  femur 
and  the  femur  length.  As  the  driver  translates  forward  during  the  crash 
these  angles  are  recalculated  so  as  to  obtain  the  proper  angle  at  which 
the  knee  loads  are  applied  to  the  femurs  and  lower  body. 

The  specifics  of  how  all  the  input  is  handled  is  given  in  Section  3.1. 
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2.6  THE  DRIVER  ALGORITHM 


The  driver  is  modeled  by  three  lumped  masses  (i.e.,  the  head, 
the  torso  and  the  lower  body)  which  pivot  at  points  A and  B in 
Figure  5.  This  figure  also  describes  the  driver  geometry  and  the 
location  of  the  airbag  and  the  driver  with  respect  to  the  compart- 
ment. Specific  details  needed  to  provide  driver  related  computer 
input  are  described  in  Section  3.1. 

A comment  is  necessary  on  the  resisting  torque  generated  by  neck 
muscular  resistance  and  the  anatomical  interferences  caused  by  the 
relative  displacement  between  the  head  and  torso.  These  values  are 
only  applied  if  certain  conditions  are  met.  For  example,  in  cases 
where  the  head  is  returning  to  become  more  nearly  in  line  with  the 
torso  (OH-0T  becoming  smaller)  the  torque  is  not  applied.  Only  when 
the  head  is  becoming  more  out  of  line  with  the  chest  (0^-Oy  becoming 


r) 

is  the  i 

neck 

torque  appl  i 

ed . 

I . 

e. , 

For 

0H 

-°T 

> 

0 

and 

6h-6t> 

0, 

T< 

0 

For 

% 

-0T 

> 

0 

and 

W 

0, 

T= 

0 

For 

°H 

-°T 
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0 

and 
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1 — 
•© 

1 

m 

•© 

0, 

T= 

0 

for 

°H 

-°T 

< 

0 

and 

W 

0, 

T> 

0 

Where, 

T = Neck  Resisting  Torque 

Ou  is  defined  in  Figure  5 

0-j.  is  defined  in  Figure  5 

0H-0-|-  is  defined  in  Figure  1 

Oh-6t  is  the  time  rate  of  change  of  O^-0-j. 
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Figure  5.  Driver  Nomenclature 
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In  addition  to  the  neck  resisting  torque  based  upon  the  angular 
displacement  of  the  head  relative  to  the  torso,  we  also  have  included 
a damping  coefficient  for  the  head/neck  based  upon  the  angular  veloci- 
ty of  the  head  relative  to  the  torso.  The  value  of  this  coefficient 
is  known  as  "DCN"  in  the  input  data. 

Thus,  the  overall  resistance  of  the  head/neck  to  rotation  with 
respect  to  the  torso  is  composed  of  two  terms  - one  term  based  upon  rel- 
ative angular  displacement  and  a second  term  based  upon  relative  angular 
velocity. 

The  injury  measures  computed  by  DRACR  for  the  driver  are  the 
fol lowing : 

a.  Head  Anterior-Posterior  g's 

b.  Head  Superior-Inferior  g's 

c.  Chest  Anterior-Posterior  g's 

d.  Chest  Superior-Inferior  g's 

e.  Femur  Loads 

f.  Head  Injury  Criterion  (HIC) 

Appendix  D,  of  this  manual  contains  a listing  of  the  DRACR  program 
and  the  logic  flow  diagrams. 
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3.0  EXAMPLE  OF  DRACR  PROGRAM  USE 


This  section  describes  the  use  of  the  DRACR  program  by  citing  an 
actual  sled  impact  test  case.  This  test  case  was  conducted  as  part  of 
the  test  effort  in  Contract  No.  DTNH22-82-C-071 32,  "Systems  Analysis 
Approach  to  Restraint  Systems  Integration",  and  is  documented  in  a test 
report^.  The  test  was  conducted  in  East  Liberty  Ohio  using  the  HYGE 
sled  test  facility  of  NHTSA-SRL/TRC . 

In  addition  to  the  above  sled  test,  vehicle  component  and  dummy 
tests  were  conducted  to  determine  the  performance  characteristics  of 
the  various  key  components  comprising  the  vehicle  restraint  and  dummy 
systems.  These  tests  included  the  following: 
o Steering  Wheel  Static  Crush  Tests 
o Steering  Column  Static  Crush  Tests 
o Knee  Bolster  Static  Crush  Tests 
o Dummy  Measurements 
o Dummy  Neck  Static  Resistance  Tests 
o Dummy  Neck  Dynamic  Resistance  Tests 
These  tests  were  conducted  at  the  NHTSA-SRL/TRC  test  facility  in  East 
Liberty  Ohio,  under  Contract  No.  DTNH22-82-P-07224,  "Sled  Testing  and 
Data  Collection  at  NHTSA-SRL,  Ohio  - DRACR  Validation". 


Fitzpatrick  Engineering,  "Systems  Analysis  Approach  to  Restraint  Sys- 
tems Integration  - Sled  Test  Support",  Contract  No.  DTNH22-82-C-071 32 , 
Interim  Report,  Vols.  I and  III. 
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3.1  CREATING  THE  DATA  INPUT  FILE 


The  first  thing  to  be  done  in  preparation  for  making  a DRACR  com- 
puter run  is  to  set  up  a data  input  file  containing  information  which 
the  program  requires  to  run.  Table  1 is  a list  of  parameters  that 
make  up  this  file. 

In  Table  1,  the  first  column  shows  the  data  file  line  number.  The 
second  column  shows  the  location  of  a particular  piece  of  data  in  the 
line  ("2",  for  example,  would  indicate  the  2n<^  position  in  the  line). 

The  third  column  gives  the  alpha-numeric  name  of  the  variable  at  that 
particular  point  in  the  line,  followed  by  the  units  the  variable  must 
have  in  the  data  file.  The  fourth  column  contains  a short  description 
of  what  the  variable  is.  The  fifth  column  gives  the  actual  value  of 
the  variable  used  in  the  sample  case.  The  file  that  results  from  this 
input  is  shown  in  Figure  6.  We  have  chosen  to  call  the  file  "LSRL5005". 

A word  about  the  subroutine  "SPRING"  is  in  order  here.  As  prev- 
iously mentioned,  "SPRING"  is  used  in  those  cases  where  one  wishes  to 
include  the  effects  of  hysterises  in  the  deforming  element.  Figure  7 
illustrates  the  logic  involved  in  modeling  this  type  of  behavior.  In 
general,  the  force  will  increase  in  a non-linear  way  as  shown  by  Path  1 
(Figure  7).  If  the  spring  begins  to  unload  at  point  A,  we  expect  a 
linear  unloading  curve  as  shown  by  path  2.  If  the  spring  continues  to 
expand  past  the  zero  force  (point  B)  we  might  expect  some  kind  of  ten- 
sion force  (negative)  to  build  up  as  in  path  3.  Reloading  of  the  spring 
again  would  cause  a force  build-up  as  in  Path  4.  If  the  previous  maximum 
compression  point  A were  then  exceeded,  a new  unloading  curve  would  re- 
sult as  in  Path  5,  parallel  to  the  old  one. 
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TABLE  1.  Input  File 


Line  Location  Variable  Name  Description  Value 


No. 

In  Line 

and  Units 

1 

1 

Y(4),  mph 

Vehicle  Impact  Velocity 

28. 

1 

2 

Y(6),  deg 

Head  Angle  (0H  - Fig.  1 ) 

-6. 

1 

3 

Y(7),  deg 

Torso  Angle  (Oj  - Fig.  1) 

-24. 

2 

1 

ZL,  lb 

Weight  of  Lower  Body  (Fig.  8) 

71  .04 

2 

2 

ZT,  lb 

Weight  of  Torso  (Fig.  8) 

58.4 

2 

3 

V lb 

Weight  of  Head  (Fig.  8) 

11.44 

2 

4 

Rr  in 

Distance  From  H-Point  to 
Torso  c.g.  (Figure  8) 

14. 

2 

5 

Rn>  in 

Distance  From  H-Point  to 
Neck  Pivot  (Figure  8) 

20.5 

2 

6 

Rh,  in 

Distance  From  Neck  Pivot 
to  Head  c.g.  When  0„  = 0T 
(Figure  8) 

4.75 

2 

7 

RT0PH,ln 

Distance  from  H-Point  to 
Top  of  Head  When  0H  = 0^ 

27.75 

3 

1 

NPN 

No.  of  Points  in  Neck 
Torque  vs  Angle  Curve 

4 

3 

2 

NKR 

No.  of  Points  in  Knee  Force 
vs  Crush  Curve 

7 

3 

3 

NV 

No.  of  Points  in  Vehicle  g's 
vs  Time  Curve 

12 

3 

4 

NSF 

No.  of  Points  in  Seat  Frict- 
ion vs  Displacement  Curve 

6 

3 

5 

NPG 

No.  of  Points  in  Gas  Flow 
vs  Time  Curve 

12 

3 

6 

NPC 

No.  of  Points  in  Column 
Force  vs  Crush  Curve 

12 
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TABLE  1 . Input  File 


Cont 'd 


Line 

No. 

Location 
In  Line 

Variable  Name 
and  Units 

Description 

Value 

3 

7 

NPW 

No.  of  Points  in  Wheel 
Force  vs  Crush  Curve 

11 

3 

8 

NPWRP 

No.  of  Points  in  Wheel  Tor- 
que vs  angle  (positive)  Curve 

6 

3 

9 

NPWRN 

No.  of  Points  in  Wheel  Tor- 
que vs  Angle  (negative)  Curve 

6 

3 

10 

NPCRP 

No.  of  Points  in  Column  Tor- 
que (or  force,  depending  on 
column  type)  vs  Angle  (positive) 
Curve 

4 

3 

11 

NPCRN 

No.  of  Points  in  Column  Tor- 
que (or  force,  depending  on 
column  type)  vs  Angle  (neg- 
ative) Curve 

4 

4 

1 

SUN,  Ib/in 

Seat  Friction  Unload  Slope 

5000. 

4 

2 

SKR,  lb/in 

Knee  Restraint  Unload  Slope 

5000. 

5 

all 

GEN ( 1 , K) , msec 

Gas  Flow  Time  Points 

(Fig. 

9) 

6 

all 

GEN(2,K) , lb/s 

Gas  Flow  Rate  Points 

(Fig. 

9) 

7 

all 

COL ( 1 ,K) , in 

Column  Stroke  Points 

(Fig. 

10) 

8 

all 

C0L(2,K) , lb 

Column  Force  Points 

(Fig. 

10) 

9 

1 

ATMOP,  psia 

Local  Atmospheric  Pressure 

14.7 

9 

2 

PGZ,  psig 

Initial  Airbag  Pressure 

-14.7 

9 

3 

GTZ,  deg  R 

Gas  Temp.  Entering  Airbag 

1660 

9 

4 

u,  in  ,bf 

lb  deg  R 
m 3 

Universal  Gas  Constant 

662 
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TABLE  1 . Input  File 


Cont'd 


Line 

No. 

Location 
In  Line 

Variable  Name 
and  Units 

Descri ption 

Value 

9 

5 

PN1 

Polytropic  Gas  Exponent, 
Flow 

1 .4 

9 

6 

PN2 

Polytropic  Gas  Exponent, 
Compression 

1.4 

9 

7 

PN3 

Polytropic  Gas  Exponent, 
Expansion 

1.4 

10 

1 

THETAW 

, deg 

Wheel  Angle  Measured  From 
the  Vertical 

15.25 

10 

2 

RIP, 

i n 

Radius  of  assumed  circular 
Inflator  and  Airbag  Package 

2.5 

10 

3 

WWH , 

lb 

Weight  of  Steering  Wheel 

3.0 

10 

4 

WIP, 

lb 

Weight  of  Airbag  and  Inf- 
lator Package 

2.0 

10 

5 

RCOL , 

in 

Distance  from  Column  Pivot 
Point  to  Center  of  Mass 

15.0 

10 

6 

DCN, 

ft-1 b-sec 

Head/Neck  Damping  Coeff. 

0.69 

11 

1 

VC1 

rad 

Vent  Discharge  Coefficient, 
Subsonic  Flow 

0.7 

11 

2 

VC2 

Vent  Discharge  Coefficient, 
Sonic  Flow 

0.7 

11 

3 

AV, 

sq.  in. 

Vent  Area 

2.07 

11 

4 

SA, 

in 

Airbag  Major  Axis  Length 

12.5 

11 

5 

SC, 

in 

Airbag  Minor  Axis  Length 

8.0 

11 

6 

X1Z, 

in 

Horizontal  Airbag  Reference 
Distance  (Figure  5) 

30.6 

11 

7 

viz. 

i n 

Vertical  Airbag  Reference 
Distance  (Figure  5) 

18.9 

12 

1 

THETAC 

, deg 

Column  Angle  Measured  From 
Horizontal 

15.25 

12 

2 

MU 

Coefficient  of  Friction 
(Figure  3) 

0.0 

12 

3 

LSCZ , 

in 

Column  Ref.  Dimension  to 
Shear  Capsule  (Figure  3) 

11.0 

12 

4 

LFWZ , 

in 

Column  Ref.  Dimension  to 
Firewall  (Figure  3) 

16.5 
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TABLE  1.  Input  File 


Cont 'd 


Li  ne 
No. 

Location 
In  Line 

Variable  Name 
and  Units 

Description 

Value 

12 

5 

LBAZ,  in 

Column  Ref.  Dimension  to 
Aft  Column  Bearing  (Fig.  3) 

15.0 

12 

6 

LBFZ,  in 

Column  Ref.  Dimension  to 
Fwd  Column  Bearing  (Fig.  3) 

19.0 

12 

7 

WC,  lb 

Column  Weight 

7.0 

13 

1 

LF,  in 

Femur  Length  (Figure  1) 

18.0 

13 

2 

THFO , deg 

Initial  Femur  Angle 
(Figure  1 ) 

12.0 

13 

3 

THLO , deg 

Initial  Tibia  Angle 
(Figure  1 ) 

35.0 

14 

1 

XWH , in 

X-Coordinate  of  Wheel  Flub 
(Pt.  "p"  in  Figure  2) 

-6.0 

14 

2 

RIMRAD,in 

Radius  of  Wheel  Rim 

7.38 

14 

3 

X2Z , in 

Horizontal  H-Point  Ref. 
Dimension  (Figure  5) 

34.5 

14 

4 

Y2Z,  in 

Vertical  H-Point  Ref. 
Dimension  (Figure  5) 

7.4 

14 

5 

WB,  in 

Width  of  Torso 

11.25 

14 

6 

WH,  in 

Width  of  Head 

5.0 

15 

all 

SFN ( 1 ,k) , in 

Seat  Friction  Displacement 
Points 

(Fig.  6) 

16 

all 

SFN( 2 , k ) , lb 

Seat  Friction  Force  Points 

(Fig.  6) 

17 

all 

FNECK ( 1 ,k) , 
deg 

Head/Torso  Relative  Angle 
(0..  - ©T  in  Figure  1) 

(Fig.  11) 

18 

all 

FNECK(2 , k) , 
ft-lb 

Neck  Rotational  Resistance 
Torque 

(Fig.  11) 

19 

all 

VEHGS( 1 ,k) , ms 

Crash  Pulse  Time  Points 

(Fig.  12) 

20 

all 

VEHGS ( 2 , k ) , g’s 

Crash  Pulse  g Points 

(Fig.  12) 

21 

all 

KRN ( 1 , k ) , in 

Knee  Displacement  Points 

(Fig.  13) 

22 

all 

KRN(2,k) , lb 

Knee  (femur)  Force  Points 
(sum  of  two  knees) 

(Fig  . 13) 

23 

all 

WHE ( 1 ,k) , in 

Wheel  Stroke  (crush)  Points 

(Fig  . 14) 

24 

all 

WHE ( 2 , k ) , lb 

Wheel  Force  Points 

(Fig.  14) 
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TABLE  1.  Input  File 


Cont 1 d 


Line 

No. 

Location 
In  Line 

Variable  Name 
and  Units 

Descri ption 

Value 

25 

all 

WHRP(1 ,k). 

deg 

Wheel  Rotation  Angle 
Points,  Upward  With 
Respect  to  Column 

(Fig.  15) 

26 

all 

WHRP(2,k) , 

in-lb 

Wheel  Rotation  Resist- 
ance Points,  Upward 
Rotation 

(Fig.  15) 

27 

all 

WHRN ( 1 ,k) , 

deg 

Wheel  Rotation  Angle 
Points,  Downward  With 
Respect  to  Column 

(Fig.  15) 

28 

all 

WHRN(2,k), 

in-lb 

Wheel  Rotation  Resist- 
ance Points,  Downward 
Rotation 

(Fig  15) 

29 

all 

COLRP ( 1 ,k) , 

deg 

Column  Rotation  Angle 
Points,  Upward  With 
Respect  to  Vehicle 

(Fig.  6) 

30 

all 

C0LRP(2,k) , 
lb  if  Type  1 
Col . , in-lb 
Otherwi se 

Column  Rotation  Resist- 
ance Points,  Upward 
Rotation 

(Fig.  6) 

31 

all 

C0LRN( 1 ,k) , 

deg 

Column  Rotation  Angle 
Points,  Downward  With 
Respect  to  Vehicle 

(Fig.  6) 

32 

all 

COLRN ( 2 , k ) , 
lb  if  Type  1 
Col . , in-lb 
Otherwise 

Column  Rotation  Resist- 
ance Points,  Downward 
Rotation 

(Fig.  6) 
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LSRL5005 


IV: 13EDT 


07/29/82 


28. . -6. .-24. 

71. 04. 58. 4. 11. 44. 14..  20. 5. 4. 75. 27. 75 
4,7.12.6.12.12,11.6,6.4.4 

5000..  5000. 

0. .17. 5. 20. 5. 26.. 31.. 36.. 41.. 46.. 51. .56.. 71.. 100. 

0.. 0.. 3. 4, 4. 09. 4. 26. 3. 6, 2. 45. 1.63.. 98.. 49. 0..0. 

-10., 0., .55. .95.1.55,1.6.2.34.2.57.3.05.3.85.4.13.5.7 
0..0.. 710., 900.. 730., 500.. 500., 930.. 680.. 860.. 1100.. 1510. 
14.7,-14.7.1660. .662.. 1.4. 1.4, 1.4 

15. 25. 2. 5. 3. . 2. .15. ..69 

.7. .7. 2. 07. 12. 5. 8. .30. 6. 18. 9 

15. 25.0.  .11. .16. 5. 15. .19. .7. 

18. . 12. .35. 

-6. .7.38.34.5,7.4.11.25.5. 

-50. .0. . 1 . . 14. . 15. ,50. 

0. .0. .250. ,250. .0. .0. 

-80. .18. .20.5,90. 

115. . 0.. -6. 7. -86. 7 

0.. 4.. 14.. 25.. 33.. 45., 57.. 67.. 73.. 102.. 118.. 200. 

0..1..10. ,11. 2. 8. 9. 13.. 22.. 15.. 17. 2. 8. 6. -.5. -.5 
-10.. 3. 15. 3. 25. 3. 58. 3. 68. 4. 13. 10. 

0. .0. .400. .900. .2100. .2500. .2500. 

-10.. 0... 23.. 5,. 91. 1.72. 1.95. 2. 56. 3.. 3. 66. 6. 

0. .0. .1075. . 1305. . 1165. .945. .990. .990. , 1 120. .680. . 1355. 
0., 1.6. 4. 9. 8. 9, 16. ,23. 3 
0 . . 3990 . , 3800 . . 3400 . . 3030 . . 2850 . 

0. .1.6. 4. 9, 8. 9, 16. .23. 3 
0 . , 3990 . , 3800 . . 3400 . . 3030 . . 2850 . 

0..5..6..7. 

0.. 50.. 2000.. 2000. 

0..5..6..7. 

0.. 50.. 2000., 2000. 


Figure  6.  Listing,  Data  File  LSRL5005 
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Figure  7.  "SPRING"  Nomenclature 
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SEGMENT  WEIGHTS  (50th  % 

Male)1 

1 . 

Head  ( incl . Neck)  . . . . 

12.54 

lbs 

2. 

Foot  (ea . ) 

2.66 

3. 

Lower  Leg  (ea . ) 

8.07 

4. 

Upper  Leg  (ea . ) 

12.31 

5. 

Pelvis  

38.35 

6. 

Thor+Shoul +Abdom  

50.58 

7. 

Upper  Arm  (ea . ) 

3.91 

8. 

Forearm  (ea.)  

3.00 

9. 

Hand  (ea.)  .'. 

1 .45 

TOT. 

164.27 

1 bs 

DRACR/PAC  WEIGHTS 

1 . 

Head  (incl.  Neck  - meas 

.)  .. 

11.44 

lbs 

2. 

Torso  (6  + 2*7) 

58.40 

3. 

Hip  (5  + 2*4  + 3)  . . . . 

71 .04 

TOT 

140.88 

lbs 

Width  = L1  = 5.0" 


Recommendations  from  Crash  Test  Dummy 
Task  Force. 


GEOMETRY  DATA 


Width  = L 


1 . 

rt 

= L1 9 

VO 

C\J 

—i 

I 

2. 

rn 

= L9  ' 

L26  = 

3. 

rtoph 

= L10 

' L26 

4. 

lflesh 

= L30“ 

L29  = 

5. 

lf 

= L32 

- (L29 

= 14.0" 


27.75" 


Figure  8.  Dummy  Measurements 
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= 3S%, 
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Figure  10.  Steering  Column  Stroking  Charac- 
teristics, Static  - DeLorean  Design. 


31 


ACCELERATION  (CO 

10.0  15.0  20.0  25.0  30.0  35. 


TRC060IF  TRC060  55 1055003  82069  14-RPn -82  15:49; 36 

FTUXG  FILTER  = BLF'F  100/  317/  -40  30.0  MPH 

MIN.  MAX  VALUES  = -2.  '49  c 157.000,  22.33  c 5".  325 


E 


I I I I L U I JLU 

S 5 1 055003  - FLQORPflN  RCC 
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DRACR  Input  Data  Points 


ioaa 


Adult  Female  Dummy 


Figure  13.  Knee  Bolster  Force  vs 
Displacement  Characteristics. 


34 


1979  Volvo  GT  Steering  Wh 
« DRACR  Input  Data  Points 


inarm**., 


Figure  15.  Steering  Wheel  Rotat 
ional  Crush  Characteristics,  S 
atic. 
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In  order  to  handle  this  type  of  situation,  certain  data  is  req- 
u i red : 

a)  We  need  to  know  the  non-linear  spring  characteri Stic  for  Path  1. 

b)  We  need  to  know  what  the  slope  of  the  unloading  curve  is  in 
paths  2 and  5. 

c)  If  we  want  to  allow  the  deforming  member  to  develop  negative 
forces  when  it  unloads  past  the  zero  force  point,  we  need  to 
know  what  the  slope  of  that  portion  of  the  unloading  curve  is 
in  Path  3 (Note:  based  on  experience,  this  slope  is  set  equal 
to  zero  in  DRACR). 

Three  subroutines  have  been  included  in  DRACR  to  model  this  plastic 
behavior.  The  subroutine  "SPRING"  models  a piecewise-1 inear  force  vs 
deflection  characteristic  stored  in  a two-dimensional  array  called  "F". 
SPRING  also  requires  the  maximum  previous  displacement,  the  current  dis- 
placement, the  two  unload  slopes,  and  the  number  of  points  in  the  array  F. 
With  these  data,  SPRING  calls  a second  subroutine  "LOOKUP"  which  uses 
the  array  information  and  the  current  spring  displacement  to  calculate  a 
force  as  if  the  spring  were  loading  up  along  the  initial  force  vs  def- 
lection curve.  SPRING  then  uses  its  built  in  logic  to  determine  the  path 
upon  which  the  member  is  loading  by  comparing  the  current  and  the  pre- 
vious maximum  member  displacements  and  then  calculating  the  correct  force 
value.  The  previous  maximum  member  displacement  is  kept  current  in  a 
third  subroutine  called  "UPDATE". 

3 . 2 RUNNING  THE  PROGRAM 

Once  the  input  file  has  been  created  and  saved,  as  shown  in  Figure  6, 
we  are  ready  to  run  DRACR.  At  this  point  the  user  accesses  DRACR  and  tell 
the  computer  to  execute  the  program.  The  computer  will  respond  by  asking 
the  user  to  name  the  input  file;  in  the  subject  case  we  respond  with  "LSRL 
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5005",  as  shown  in  Appendix  E of  this  manual.  The  computer  will  then  ask 
the  user  to  type  "1"  if  the  steering  column  is  supported  at  the  firewall 
(as  are  most  conventional  columns),  or  "2"  if  the  column  is  not  (as  is  the 
case  with  some  of  the  newer  vehicles  such  as  the  GM  X-Body  cars  and  the 
Chrysler  front  wheel  drive  cars  such  as  the  OMNI  and  Horizon  where  the  col- 
umn is  supported  further  aft  by  the  A-post  and  plenum  respectively). 

Figure  16  shows  a Type  1 column  while  Figure  17  shows  a Type  2 column. 
If  the  user  is  simulating  a Type  1 column,  the  support  bracket  strength  at 
the  shear  capsule  location  is  specified  in  Line  Nos.  30  and  32  of  the  input 
file  (Table  1)  in  units  of  pounds . Thus,  for  this  case  the  applied  moments 
are  reacted  by  a couple  between  the  firewall  and  the  shear  capsule.  If, 
however,  the  user  is  simulating  a Type  2 column  where  no  firewall  reaction 
exists,  the  support  bracket  strength  is  specified  in  inch-pounds . This  is 
because  all  of  the  applied  moment  to  the  column  is  assumed  to  be  reacted  by 
this  single  bracket.  Again,  Line  Nos.  30  and  32  are  where  this  bracket 
strength  is  specified  in  the  input  file  for  the  angular  displacement  values 
specified  in  Line  Nos.  29  and  31. 

Once  these  answers  are  given  by  the  user,  the  computer  will  begin  to 
print  the  input  data  followed  immediately  by  the  desired  output.  Let  us 
now  discuss  this  output.  Altogether,  there  are  eleven  blocks  of  output 
data,  each  block  consisting  of  the  amount  of  output  that  can  be  conveniently 
grouped  together  (in  terms  of  subject  matter)  on  a single  line.  The  con- 
tents of  each  block  will  now  be  described.  Appendix  E contains  the  actual 
output  for  the  sample  run. 

Block  1 is  basically  a summary  of  some  of  the  main  parameters  of  interest 
and  consists  of  the  following  items  as  we  proceed  left-to-right  across  the 
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TYPE  " 1 " STEERI NG  COLUMN 


Column  Supported  at  the  Firewall  and  some  Intermediate  Point. 


Figure  16.  Tvne  "1"  Steerino  Column  Configuration. 


Column  Supported  Entirely  at  some  One  Point. 


Figure  17.  Type  "2"  Steering  Column  Configuration 


-39- 


page:  Elapsed  Time,  Vehicle  G's,  Vehicle  Velocity,  Vehicle  Displacement 

(crush).  Body  G's  (chest  A-P  g's).  Column  Displacement  (stroke),  and  Air- 
bag Pressure. 

Block  2 Contains:  Elapsed  Time,  Femur  Angle  (from  horizontal,  positive 

up),  H-Point  Velocity  (with  respect  to  ground),  H-Point  Acceleration, 

Femur  Force  (per  femur).  Seat  Friction  Force,  H-Point  Displacement  Rela- 
tive to  the  compartment. 

Block  3 contains:  Elapsed  Time,  Torso  Displacement  (horizontal  with  re- 

spect to  ground).  Torso  Angle  (Oy  in  Figure  1,  positive  when  toward  the 
dash).  Torso  Angular  Velocity  (positive  toward  the  dash),  Torso  Angular 
Acceleration  (positive  toward  the  dash),  Torso  Horizontal  Displacement 
Relative  to  Compartment,  Torso  Horizontal  Velocity  Relative  to  Compart- 
ment. 

Block  4 contains:  Elapsed  Time,  Head  Displacement  (horizontal  with  re- 

spect to  ground).  Head  Angle  (Qu  in  Figure  1,  positive  toward  dash). 

Head  Angular  Velocity  (positive  toward  dash).  Head  Angular  Acceleration 
(positive  toward  dash).  Head  Horizontal  Displacement  Relative  to  Compart- 
ment, Head  Angle  Relative  to  Torso. 

Block  5 contains:  Elapsed  Time,  Axial  Force  Applied  to  Steering  Wheel, 

Normal  Force  Applied  to  Steering  Wheel,  Bending  Moment  Applied  to  Steering 
Wheel  (positive  when  in  a direction  that  would  rotate  top  of  wheel  toward 
windshield).  Steering  Wheel  Resistance  to  Crush  (in  axial  direction),  Steer- 
ing Wheel  Crush,  Steering  Wheel  Crushing  Velocity  (in  axial  direction,  Rel- 
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ative  to  compartment). 


Blj)cJ<  6 contains:  Elapsed  Time,  Axial  Force  Applied  to  Steering  Column, 

Normal  Force  Applied  to  Steering  Column,  Bending  Moment  Applied  to  Steer- 
ing Column  (positive  when  in  a direction  that  would  rotate  column  upward 
toward  windshield).  Steering  Column  Resistance  to  Stroke  (in  axial  dir- 
ection), Steering  Column  Stoke,  Steering  Column  Stroking  Velocity  (in  axial 
direction,  relative  to  compartment). 

Block  7 contains:  Elapsed  Time,  Bending  Moment  Applied  to  Steering  Wheel 
(positive  when  in  a direction  that  would  rotate  top  of  wheel  toward  wind- 
shield), Torsional  Resistance  to  Steering  Wheel  Rotation  (same  sign  con- 
vention), Steering  Wheel  Angle  (measured  from  vertical,  positive  when  to- 
ward the  windshield).  Steering  Wheel  Angular  Velocity  (positive  when  top 
of  wheel  moving  toward  windshield). 

Block  8 contains:  Elapsed  Time,  Rotational  Torque  Applied  to  Steering 

Column  about  Column  Pivot  Point  (positive  when  it  would  cause  an  upward 
rotation  of  the  column).  Rotational  Resistance  to  Steering  Column  Rotation 
(same  sign  convention).  Steering  Column  Angle  (measured  from  horizontal, 
positive  when  column  tilted  upward).  Steering  Column  Angular  Velocity  (pos- 
itive when  column  rotating  upward). 

Block  9 contains:  Elapsed  Time,  Airbag  Penetration  (measured  normal  to 

the  torso,  halfway  between  the  two  bag  intercept  points),  Airbag  Volume, 
Airbag  Pressure,  Airbag  Wraparound  (fabric  tension)  Force  (applied  to  the 
chest),  and  Airbag  Pressure  Force  (applied  to  the  chest). 
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Block  10  contains:  Elapsed  Time,  Chest  A-P  Acceleration,  Chest  S-I 


Acceleration,  Head  A-P  Acceleration,  Head  S-I  Acceleration. 

Block  11  contains:  Energy  Absorbed  if  Steering  Wheel  Crush  Occurs, 

Energy  Absorbed  if  Steering  Column  Crush  (stroke)  Occurs,  Energy  Absorb- 
ed if  Steering  Wheel  Rotation  Occurs,  Energy  Absorbed  if  Steering  Column 
Rotation  Occurs,  Energy  Actually  Absorbed  for  the  mode  that  did  occur. 

After  all  of  the  above  listed  output  data  has  been  printed,  the 
computer  will  ask  the  user  to  type  "1"  if  he  wishes  the  value  for  HIC 
to  be  computed.  Typing  "1"  will  cause  the  computer  to  proceed  and  cal- 
culate the  value  for  HIC.  Typing  any  other  value  or  typing  a carriage 
return  will  cause  the  computation  to  cease  without  computation  of  the 
value  for  HIC. 

This  completes  the  presentation  of  the  information  needed  to  use 
the  DRACR  program. 

3.3  CORRELATION  OF  DRACR  RESULTS  WITH  TEST  DATA 

This  section  compares  the  DRACR  results,  for  input  File  LSRL5005 
(Figure  6)  with  actual  sled  impact  test  data.  Data  File  LSRL5005  rep- 
resents the  vehicle,  restraint  and  dummy  data  for  sled  test  No.  S51055003, 
which  was  conducted  as  part  of  Contract  No.  DTNH22-82-C-071 32 J This  test 
was  conducted  at  30  mph  with  a 50th  percentile  adult  male  ATD  driver. 


1 

Ibid. , page  19. 


-42- 


Some  of  the  data  in  File  LSRL5005  are  based  on  static  test  data 
obtained  under  Contract  No.  DTNH22-82-P-07224,  "Sled  Testing  and  Data 
Collection  at  NHTSA-SRL,  Ohio".  The  data  from  these  tests  include: 

o Steering  Wheel  Static  Crush  Characteristics  - Axial  (Figure  14). 

o Steering  Wheel  Static  Crush  Characteristics  - Rotational  (Figure  15). 

o Steering  Column  Crush  Characteristics  - Axial  (Figure  10). 

o Knee  Bolster  Crush  Characteristics  (Figure  13). 
o Dummy  Measurements  (Figure  8) 
o Dummy  Neck  Stiffness  - Static  (Figure  11). 

Other  data  in  File  LSRL5005  were  obtained  from  the  actual  sled  impact  test. 
These  data  include: 

o Vehicle  Acceleration  History  (Figure  12). 
o Dummy  Pre-test  Measurement  (Figure  18). 
o Impact  Speed. 

One  point  that  is  worth  noting  here  is  that  in  making  the  DRACR 
simulation  no  attempt  was  made  to  "fine-tune"  the  input  data  to  achieve 
a better  correlation  with  the  test  data.  We  simply  used  the  available 
data  without  any  modifications.  The  differences  that  resulted  between 
the  DRACR  predictions  and  test  results  are  discussed  below.  Figures  19 
thru  27  compare  the  DRACR  predictions  with  the  test  results.  Table  2 
summarizes  the  DRACR  vs  test  results. 

The  differences  between  the  DRACR  predictions  and  the  test  results, 
as  summarized  in  Figures  19  thru  27,  have  been  attributed  to  the  following 
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Figure  18.  Pre-test  Measurements,  Test  #5. 
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Figure  19.  Airbag  Pressure 
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CELERRT I ON 
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Figure  23.  Head  S-I  Acceleration 
History. 
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RCCELERRTION 
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TABLE  2.  Summary  of  DRACR  v.s.  Test  Results 


I tern 

DRACR 

Test 

Pk.  Head  Resultant  q's  (-3  msec)  

23.0 

27.5 

Pk.  Torso  Resultant  g's  (-3  msec)  .... 

22.0 

28.96 

Pk.  Femur  Loads  (lbs) 

. Left  

1120 

1080 

Right  

1120 

1200 

HIC  

68.3 

115 

T1  (msec) 

50.0  1 

100.0 

49.0 

T2  (msec ) 

128.5 

Column  Displacement  (inch)  

2.92 

2.75 

Wheel  Displacement  (inch)  

.26 

.22  (average) 

Pk.  Bag  Pressure  (psiq)  

6.45 

6.0 

1 

Run  terminated,  torso  contacted  rim. 
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The  differences  between  the  DRACR  predictions  and  the  test  results, 
as  summarized  in  Figures  19  thru  27,  have  been  attributed  to  the  following 
factors : ^ 

1.  The  actual  shape  of  the  ACRS  bag  was  a "tear-drop"  rather  than 
the  idealized  ellipsoidal  bag  shape  assumed  in  DRACR.  The  actual 

tear-drop  shape  gives  more  head  support  and  less  initial  bag 
penetration  (by  the  torso)  than  the  idealized  DRACR  shape. 

2.  The  torso  of  the  dummy  appears  to  have  contacted  the  steering 
wheel  rim  approximately  80  msec  into  the  impact  event.  DRACR 
does  not  account  for  this  interaction. 

3.  The  "femur  loads"  calculated  in  DRACR  are  actually  knee  con- 
tact loads.  These  loads  are  different  than  the  measured  fe- 
mur loads  in  the  sled  test  in  that  they  include  the  inertial 
loads  corresponding  to  the  total  lower  body  mass  of  the  occu- 
pant. The  measured  femur  loads,  on  the  other  hand,  do  not 
include  the  inertial  loads  corresponding  to  the  lower  body 
mass  of  the  occupant  located  forward  of  the  load  cell.  Thus, 
there  will  always  be  some  discrepancy  between  the  DRACR  "femur" 
loads  and  the  measured  femur  loads. 

4.  DRACR  restricts  the  movement  of  the  ‘ H ’ point  to  the  horizontal 
plane.  In  the  actual  sled  tests  the  dummy's  * H * moved  down- 
ward approximately  1-2  inches.  This  can  have  an  effect  on  the 
torso  S-I  accelerations. 


Based  on  the  conclusions  from  the  research  work  under  Contract  No.  DTNH- 
22-82-P-07224 , "Sled  Testing  and  Data  Collection  at  NHTSA-SRL , Ohio". 
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The  differences  between  the  DRACR  predictions  and  the  HYGE  sled 
impact  test  results,  as  summarized  above,  suggest  that  additional 
testing  need  to  be  conducted  to  validate  the  DRACR  model.  We  recomm- 
end a test  program  specifically  designed  for  this  validation  effort 
(rather  than  sharing  tests  with  an  ongoing  development  effort,  as  was 
the  case  for  the  tests  referenced  herein).  In  this  validation  test 
program  we  would  recommend  additional  instrumentation  to  monitor  the 
steering  wheel  and  steering  column  dynamic  response  as  well  as  the 
hip  response  of  the  dummy. 
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APPENDIX  A 


DERIVATION  OF  THE  EQUATIONS  OF  MOTION 

The  derivation  of  the  equations  of  motion  will  be  formulated  utilizing  Lagran- 
gian  techniques  based  upon  the  geometrical  representation  in  Figure  A-l. 

Writing  an  expression  for  the  total  kinetic  energy  of  the  occupant,  we  have: 

(1)  T = % (*£  ♦ Yj)  ♦ Mt  (Xj  ♦ Y*)  + Mjc*] 

Note  that  Y.  = 0,  as  no  movement  normal  to  the  X-direction  is  allowed  for  the 
hip-leg  mass. 

= Head  mass 

Mp  = Torso  mass 

Mp  = Hip -leg  mass 

Xp  = Horizontal  translation  of  the  hip-leg  mass  with  respect  to  inertial 
reference  point  -which  is  positive  when  it  is  in  direction  shown. 

Xp  and  X^  are  similarly  defined 

= Vertical  distance  from  H -point  to  the  center  of  gravity  of  the  head 
Y,p  = Vertical  distance  from  H-point  to  the  center  of  gravity  of  the  torso 
Successive  dots  indicate  velocity  and  acceleration,  respectively. 

Writing  the  transformation  equations,  we  have: 
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*L  * rT  SIN9T 

(3) 

-T' 

ii 

rT  COS0T 

(4) 
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*L  * rN  SINeT  + rH  SIN6H 

(5) 

it 

rN  C0SeT  * rH  COS0H 

(6) 

n 

\ * rT  C0S6I  9T 

(7) 
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-r_  SIN0_  L 
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MATHEMATICAL  MODEL 


Figure  A-l.  Geometrical  Representation  of  the 
Driver-Airbag  Interaction 
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(8) 


*H  * *L  rN  C0SeT  *T  + rH  C0S9H 

(9)  Yh  = -rN  SIN0t0t  - rH  SIN0H  0^ 

where : 

r^.  = Distance  from  hip  H-point  to  torso  center  of  gravity  (The  H-point 
is  assumed  to  be  coincident  with  the  hip-leg  center  of  gravity.) 

r^  = Distance  from  H-point  to  neck  pivot  point 

r„  ■ Distance  from  the  neck  pivot  point  to  the  center  of  gravity  of  the 
H head 

0^  and  0^  are  as  defined  in  Figure  A-l. 

Substituting  Equations  6 through  9 into  Equation  1,  we  have: 

(10)  T = 4 {m^  * Hj  [x*  - 2 COS6T  9t  *46*]  ♦ ^ ♦ 

2 XL  (rN  C0S6t  5t  * rH  COS0H  9„)  - 2 rN  rH  (C0SeT  C0S6H  ^ e„  * 

SIN9-T  SIN8H  ®T  V * rN  ®T  * rH  ®h]| 

Hie  potential  energy  portion  of  the  Lagrangian  is: 

(11)  VT  - Mj.  g rT  COS0T 

U2)  VH  - ^ g (rH  C0S9h  * rN  COS0T) 

Note:  The  applied  forces  and  moments  will  be  treated  separately  later  on. 

Writing  the  Lagrangian,  we  have: 

(13)  L = T-V  = T - (VT  ♦ VH) , 

where  the  values  to  be  substituted  into  this  equation  are  given  by  Equations  10 
11  and  12. 
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The  basic  equation  in  Lagrangian  mechanics  is: 


(14) 


— ( ^ \ _ _ p 

dt  \ 9q~  / qi 


where : 


q^  = generalized  displacement  of  the  i"^1  mass 
= generalized  velocity  of  the  i1^1  mass 
F . = generalized  force  acting  on  the  i^  mass 


Taking  the  required  derivatives  from  Equation  13  for  substitution  into  Equa- 
tion 14,  we  obtain: 


(15)  aL 


3%  ’ * *T  * V *L  * *r  rT 


cos6T  eT 


+ “h  (rN  C0SeT8T 


+ rH  C0SeH  V 


(Ml  ♦ Mj.  + \ - Mj.r?  SIN0t  0*  - (rN  SIN0T  0*  + 

rH  SIN9h  0>)  ♦ rT  COS0T  ^ ^ (rN  COS0T  ^ * 

rH  COS0H  V 

(17)  3L_  . 

3XL  ' 

(18)  3L  . ^ ^ r COS0  * r>  6 ) ♦ ^ rN  COS0T  * rN  rH  (COS0T  COS0H  0R 

3fl_  L 

T 

+ SIN0t  SIN0h  0h)  + r*  0T1 


(16)  d / aL  \ 
dt  \ a^  / 
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(19) 


5F  (if)  - "r  <*L  rT  C0SeT  - \ rT  S1N6T  ®T  * *T  6t) 


+ ”h  [ *L  rN  COS0T  ' *L  rN  SIN6I  ®T  ' rN  rH  (SIN9T  9T  C0S6H  ®H 

* coseT  siNe  0*  - coseT  coseH  eH  - coseT  eT  • siNeH  eH 
- SIN0t  coseH  9*  - SIN9t  SIN6h  eH)  ♦ r*  9t  1 


(20)  3L 


39  ' ^ *L  rT  SIN9T  6T  - % rN  fSIN0T  0T  \ + rH  SIN9T  C0S9H  * 

T 

0T  0H  - rH  COS0T  SIN0h  0t  0H)  ♦ Mj.  g rT  SIN0T  + ^ g rN  SIN0T 
f211  3L  f • 

W~  = ^ \ rH  COS0H  + rN  rH  (COS0T  mSQH  0T  + SIN0T  SIN0H  V 


H 


+ r2  0 
H H 


(22)  d / 8L 


d?  ( ir)  ■ ^ *1  rH  C0SeH  - *L  rH  SIN9H  ®H  • rN  rH  ^SIN9T  C0S9H  9 
H 

+ COS0T  SI)N0H  0t  0h  - COS0T  COS0H  0T  - COS0T  SIN0h  0^ 


- SIN0t  COS0h  0?0h  - SIN0t  SIN0h  0 ) + r. 


H 9h] 


(23)  8L 


80 


H 


”h  rH  [SIN9H  9H  *L  * rN  tC0S6T  SIN9H  9T  9H  - SIN9I  C0S9H 
9T  V]*  ^ « rH  SIN6H 


Substituting  Equations  16  and  17  into  Equation  14,  we  have: 

(24)  (Mj_  + Mj.  ♦ 5^  - Mj.  rT  SIN0t  9*  - (rN  SIN0.J.0*  - rH  SIN9H  9*) 

+ ^ rT  COS9T  0 * ^ (rN  COS0T  0 * rH  COS6H  0^  - 


which  is  the  equation  of  motion  for  mass  M^. 
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Substituting  Equations  19  and  20  into  Equation  14,  we  have: 

Mj.  (XL  rT  COS0T  - y rT  SIN8T  0T  + r?2  ©T)  + XL  rN  COS0T 

- ^ rN  SIN9t  6t  - rN  rH  (SIN9T  9?  COS9H  9H  * COS9T  SIN9H  e’ 

- coseT  C0S9H  eH  - cose  eT  sin6h  eR  - sme  coseH  e* 

- SIN0t  SIN0h  0h)  ♦ r2  0T]  + Mr  ^ rT  SIN0T  ^ ^ ( 

SIN0T  eT  \ + rH  SIN0T  COS0h  0t  eH  - rH  COS0T  SIN0h  0t  0h) 

- Mj.  g rT  SIN0t  -y  rN  SIN0T  = FQT  . 

Rewriting  the  above  yields: 

(25)  Mj.  (Xl  rT  COS0T  + r2  0T)  + 1^  Xj^  rN  COS0T  - rN  rH  (COS0T  SIN0H  02 

- COS0T  COS0H  0H  - SIN0T  COS0h  02  - SIN0t  SIN0h  0 ) 

+ rN  8t]  - Wr  g rT  SIN0T  ‘ “h  g rN  SINeT  = FeT  ' 

which  is  the  equation  of  motion  of  the  torso  mass. 

Substituting  Equations  22  and  23  into  Equation  14,  we  have: 

MH  [*L  rH  COS0H  - *L  rH  SIN9H  9H  - rN  rH  (3IN6T  C0SSH  ®T 

+ cose  siNe  e eH  - coseT  coseH  eT  - coseT  siNeH  e* 

- SIN0X  coseH  0T  bH  - Sim7  simH  eT)  ♦ ^ eH  j ♦ ^ rH  [ 

sin0h  ®H  \ rN  (coseT  siNeH  eT  eH  - siNeT  coseH  y] 


% 8 rH  SIN0H  F0H 
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Rewriting  the  preceding  yields: 

(26)  ^ [XL  rH  COS0H  " rN  rH  (SIN0T  COS0H  0T  ' COS0T  COS0H  0T 
- COS0T  SIN6h  02  - SIN0t  SIN0h  0?)  + Qh  j 

" ^ 8 rH  SIN6H  = F0H  ’ 

which  is  the  equation  of  motion  for  the  head  mass. 

Writing  Equation  24  in  terms  of  XL,  we  have: 

(27)  - . i f fxl  * <*r  ri  + ^ V smeT  ^Vh  ®h  SIN6h 

W eT  C0SeT  - rH  eH  } - 

Writing  Equation  25  in  terms  of  9^,  we  have: 

0T  " rT  + NIh  rN  (p0T  " ^ rT  + ^ V \ COS0T 
- rN  rH  £ 0H  (COS0H  COS0T  + SIN0H  SIN0T) 

♦ e*  C-sin0h  COS0T  + COS0H  sin0t)J 

♦ Hj.  g rT  SIN0t  «■  Mn  g rN  SIN0T  1 . 
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Writing  Equation  26  in  terms  of  8„,  we  have: 


(29)  9h  ■=  0H 


XLC0SeH  -rN 


Vh 


H 


‘H 


(COS0H  COS0T  + SIN9h  SIN0t)  0t 


+ (SIN0h  COS0t  - COS0H  SIN0t)  02 


+ JL_  SIN0h 
rH 


where: 


F0H  = FH  rH  * T NECK 

F9T  ■ FH  C0S(9H  - 6I)  rN  * FI  rF  ' T ““ 
FXL  - FH  C0SeH  + FT  C0SeT  * FL  • 
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APPENDIX  B 


DERIVATION  OF  THE  AIRBAG  ALGORITHM 


The  angles  and  distances  described  in  this  appendix  will  be  those  depicted  in 
Figure  B-l. 


COORDINATE  TRANSFORMATION  EQUATIONS 

(1)  y = y2  ♦ y" 

(2)  x * + x" 

(3)  y = y + y'  COS0w  + X'  SIN0^ 

(4)  x = x + x’  COS0  - y'  SIN0 

1 vJ  1 W 


f 


To  obtain  transformation  equations  for  x"  and  y”  into  the  x’,y'  system, 
substitute  Equations  1 and  2 into  Equations  3 and  4 to  get: 


(5) 


x'  = 


X2 " xi + x" * y* 


COS0 


vv 


(6) 


y'  = 


y2  ’ yi  + yM  * x*  SINe 


w 


COS0 


w 


Assume  that  the  torso  may  be  represented  by  a plane  that  intersects  the  airbag 
at  line  A-B  on  the  plane  of  symmetry  of  the  airbag  (as  shown  in  Figure  B-l) . 
Assume  further  that  the  airbag  is  an  ellipsoid  whose  plane  of  symmetry  in  the 
X-Y  plane  is  as  shown  in  the  figure.  Our  job  now  will  be  to  derive  an  equation 
for  the  bag  intercept  points  in  the  x'-yr  coordinate  system. 

In  the  x"-y"  system  the  equation  for  line  A-B  is: 

(7)  y"  = mx"  + b 


Substituting  Equation  7 into  Equation  6 yields: 


(8)  y* 


y2  - y^  + mx"  + b - x1  SIN©W 


COS0, 


w 


B 


1 


COtPARTKf  (fT , BA6.  DRIVER  COORDINATE  SYSTEM 


Figure  B-l.  Geometrical  Representation  of  the  Airbag  Parameters 
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Substituting  in  Equation  8 for  x"  from  Equations  2 and  4: 


(9)  y'  = 


y - y,  + m (x'  COS0  - x_  + x,  - y*  SIN0  ) - x*  SIN0  + b 

7 2 7 1 ~ 2 1 7 w ** 


COS0 


w 


which  is  the  desired  equation  in  the  x',y'  system. 

Let  ^ ‘ + b ~ m (x^  - x^)  = B (a  constant)  and  solve  Equation  9 for  y* 

y'  (COS0^ + m SIN0  ) = B + x'  (m  C0S6^  - SIN0J 


(10) 


y.  = 


B + x*  (m  COS0  - SIN0  ) 

w Vv 

COS0  + m SIN0 

w w 


(equation  for  A-B  in  x’ ,y*  system) 


The  equation  for  the  airbag  in  the  x'-y'  system  is: 


(ID 


v»2 

— + £—  = ! 


2 2 
c a 


Substituting  Equations  IQ  Into  11  and  collecting  terms: 

xr2  [a2  (C0S8  + m SIN0  )2  + c2f-SIN0  + m COS0  )2] 

+ 2 Bc2x'  (m  COS0w  - SIN0  ) + B2  c2  - a2  c2  (COS0^  + m SIN0J2  = 0 

which  is  a quadratic  equation  in  terms  of  x'. 

Let  A = a2  (COS0  + mSIN0  )2  + c2(-SIN0  + m COS0  )2 
W vs / W w 

D = 2 B c2  (m  COS0  - SIN0  ) 

w w 

E = B2  c2  - a2  c2  (COS0  + m SIN0  )2 

W \A< 

Ax'2  + Dx'  + E = 0 


(12)  x'  = 


■D  + Y D - 4AE 
2A 


Values  for  x'  obtained  with  (12)  when  substituted  into  (10)  will  give  the  corres- 
ponding values  for  y’.  We  now  have  defined  the  line  of  intercept  (A-B)  of  the 
occupant’s  body  with  the  mid-plane  of  the  airbag. 

With  this  line  now  established,  we  can  begin  to  calculate  the  restraint  forces 
that  will  be  applied  to  the  driver. 

Forces  will  now  be  calculated  due  to  pressure  effects  (Figure  B-2) . The  force 

on  the  head  and  chest  are  composed  of  two  ccssponents  - a pressure  component  and 
a "wrap-around"  component  due  to  fabric  tension;  i.e. , 

(13)  F„_„_  = F_  + F_ 

v J CHEST  P FT 

c c 


(14)  Fir_._  = F„  + F_ 

^ ' HEAD  P„  FT 

H H 


The  pressure  forces  act  normal  to  the  head  and  chest: 


(15)  fp  ■ p wb  (rn  - \ag> 


% - W < “ 


= P w,  AB 
b 


tRN  - W i “ 


(i6)  F PwH  (XT0?U  - 
H 


(RT0PH  " ^AG^  < ^ 


PWH  ^ ' (RN  ' W3 


(RT0FH  ' ^AG^  - ^ 


where  the  pressure  P must  be  calculated  due  to  bag  volume  and  thermodynamic 
effects. 


The  fabric  tension  component  will  be  calculated  later.  Let  us  now  calculate 
the  body  moments  caused  by  these  forces.  Using  the  H- point  and  neck  pivots  as 
our  reference  points: 

(17)  F0t  = FGHEST  * HjT 

(18)  F - • RHeaj)  » 

H 

where  F and  F are  given  by  Equations  13  and  14.  We  will  now  evaluate 
CHEST  HEAD 

Rpr  RHEAD” 
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Figure  B-2.  Head-Bag  Interaction  Forces  Diagram 
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In  order  to  solve  for  R^,  we  must  derive  equations  for  the  H-point  location  in 
terms  of  the  x'-y'  coordinate  system  (see  Figure  B-3) . 


From  the  geometry  of  the  mid-plane  of  bag  impact,  the  H-point  coordinates 
are: 


(19)  x'H  = (y2  - yp  SINO^-f  ^ C0S9w 

(20)  y'H  = cos©^  ty2  - y1  - (X2  ■ *x)  tan©j 

The  equation  for  is: 

(21)  Vr  = /^x'ft  " XV  + ^'ft  ’ y 'h^  * 

where : 


(22)  x’ 


(x’a  + x W 


FT 


(23)  y* 


y ' + y 1 

7 A 7 NECK 


FT 


The  equation  for  is: 


(24) 


(25) 


_ ^TQPH  ~ % 
HEAD  = 2 

* y g - «N 

HEAD  2 


**  + “BAG  > ^ 

* 


This  derivation  completes  the  solution  for  terms  needed  for  pressure  force  and 
body  moment  computation.  We  must  now  derive  equations  for  the  fabric  tension 
component  of  bag  force  due  to  bag  wraparound  in  the  lateral  plane.  (No  wrap- 
around in  the  vertical  plane  is  considered,  since  the  body  is  generally  as  long 
as  the  bag  is  high,  so  no  wrap-around  will  occur.) 


B - & 


I 

i 


1 


B 
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Figure  B-3.  H-Point  Location  in  x'-y'  Coordinate  System 


The  figure  on  the  right  is  a 
view  looking  down  at  the 
deformed  bag. 


Let  us  now  consider  the  body  wraparound  forces  caused  by  fabric  tension.  This 
component  of  force  is  influenced  most  by  bag  pressure,  bag  penetration  and 
body  to  bag  width. 

The  force  N is  the  tensile  force  in  the  bag. 

<P 

Wb 

At  2 = — , N is  obtained  by  a force  balance. 

L (p 


N x perimeter  of  Section  A-A  = P x Area  A-A , 

<f> 

or 


, P x Area  A-A  (force  per  unit 

<|)  Perimeter  A-A  length  of  AB) 


w. 


Solve  for  a'  and  c'  @ z = 


We  may  write  for  the  y',z?  plane, 

J2  . /2  2 
(27)  y + z = a 


F Wb 

For  z 3 — 


(28)  y 


' = a'  * yj 


w 2 
2 Wb 


b - e 


For  the  x,z  plane: 


Assume  the  area  of  Section  A-A  varies  linearly  with  bag  penetration  from 
its  initial  value,  ira’c',  to  zero  when  fully  compressed  at  a bag  penetration 
of  GT.  Then: 


(3°)  Aa.a  = *a'C  (1  - - ) . 


where  BP  = bag  penetration  perpendicular  to  and  at  mid  point  of  torso 
GI  * length  across  bag  in  BP  direction. 

The  perimeter  is  given  by: 


(31) 


PER.  . ^ 2tt 
A-A  — 


I a.'^  + c' 

V 2 


(an  approximate  formula  with  accuracy 
sufficiently  close  to  the  exact  formula 
which  involves  elliptic  integrals) 


Substituting  (30)  and  (31)  into  (26) : 


(32)  N 


a'c'P 


(1 


Substituting  (28)  and  (29)  into  (32) : 
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which  is  the  equation  for  the  tension  force  in  the  ellipsoidal  airbag  in  terms 
of  the  bag  pressure,  the  lengths  of  the  major  and  minor  axes,  the  bag  penetra- 
tion, and  the  driver  body  widths; in  units  of  force  per  unit  length  SB. 


The  bag  perimeter  in  the  plane  at  the  mid- 
point of  AB  at  an  angle  to  to  the  x-axis  can 
be  found;  see  figure  at  right.  This 
perimeter  will  remain  constant  and  the  wrap- 
around configuration  must  maintain  this 
perimeter.  To  find  the  perimeter  of 

ellipse  cut  by  B-B,  one  must  follow  the 
steps : 

1)  Find  the  distance  GT  (which  will 
be  the  length  of  one  axis) 

2)  Find  the  mid-point  of  CT 

3)  Find  the  other  axis  length. 


Calculate  BP  and  GI 

To  derive  the  equation  for  CT: 

(34)  Slope  of  GI  = m^^  • 


The  point  it  goes  through  is  Xpp,  y^. 
Then,  writing  the  equation  for  ST, 

(35)  Cy'  - /„)  = "W*’  - XFTI  ’ 
or,  rewriting  the  equation, 

(36)  y*  = * XFT  * 


The  equation  for  the  mid-plane  of  the 
ellipse  is 


(37) 


Substituting  y ' into  the  above  and  collecting  terms,  one 


(38) 


+ 


( ^ ^ ) 


obtains 
(7FT  ‘ ^AB 


-1  = 0 


B 
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Let 


(39) 

A1  = 

i 

PAB2 

c 

2 

a 

B1  - 

2 m'pAB 

“7“ 

Opr  m PAB  XFTJ 

Cl  = 

(yFT  ' 

, ^2 
ra  PAB  XFTJ 

2 1 
a 

(40) 

*G  = 

-B1  + V 

fBl2  - 4A1C1 

(41) 

V S3 

-B1  - V 

fBl2  - 4A1C1 

XI 

~JKl 

(42) 

GI  = 

V(yi 

- yG)2  * (xr  - xG)2 

(43) 

BP  = 

t/(yFT 

- yG)  + (xpp  - xG) 

At  the  midpoint  of  ST, 

* _ *G+XI 

T4GI  2 

y - 1*1 
yMGI  2 

In  the  plane  x'  y* , 

,2  ,2 

c a 


for  x'  = 


*»CI 


we  have: 


yG  "^PAB  ^G  " XFT-)  + 7 FT 

yI  "'PAB  ^XI  " XFT^  + 7 FT 


y 
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Now,  in  the  plane  y-z,  away  from  the  y',  zr  axis, 


/ 2 

2 2 2 / , *MGI 

z + y = a II-  ■ ^ 


or,  for  y = y^j. 


2 

XMGI 

2 
c 

which  is  the  length  of  the  other  axis. 
Finally, 


(45) 


PerB-B  = 2lT 


V 


CGI/ 2) 2 + z2 


where  GI  and  z are  given  by  Equations  42  and  44,  respectively. 
For  constant  perimeter  for  section  B-B, 


SBR  = 2 


V(¥T* 


DBR 


PerBB  ' Wb ' SBR 
C46)  ■ R/ 


(RIM  » RIM  chord  length  at  y^ ; i.e. , RIM 


, where  R 


RIM  DIA  . 
— 2 
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t IM 


(47)  2R  SIN(180  - <f>/2)  = 2R  SIN<J>/2  = 
<P  <P 


GI  - BP  - DBR 
COSa 


The  fabric  tension  force  can  now  be  calculated. 
(50)  ~ 


F_  = 2N  AB  COS(180  - (0/2  + a)) 
r L <p 


= -2N  AB  COS(<p/2  +a)  • 

9 

= -2N  AB  COSB , 

<P 


Let  ft  = <\)/2+*( 


where  N is  given  by  Equation  33.  AB  is  given  by  AB  = -%/(x. 
<P  1 A 


V2  + *a  ' vz 


<p  is  given  by  Equation  49,  and  a is  given  by  Equation  48. 


This  completes  the  restraint  force  computation.  It  remains  to  compute  the 
volume  change  as  a function  of  bag  penetration. 


Volume  Analysis 

(51)  Area  of  Intercept  = AOI  = 


where, 


X ^BAG  " ^INE 
AB 

for  a given  y ' between  yfi  and  y^. 


B 
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For  x 


BAG1 


*RAG 


= c * 1 - 


y' 


(Equation  for 
bag  original  shape) 


BAG  = 


(52)  BAG  = 


I ’l  ‘ *' 


2a 


-y* 

l7B 


•’  V a2  - y'2  + a^SIN_i  , , 

a/Jy; 


For  XLINE  , using  analytic  geometry, 
AB 


XLINEAB  = x 


'X\  ’ *!» 


A \yA  ‘ yB  A 


‘A  ~ *B 

(yA  ■ yB5 


y’ 


Therefore, 


yB 


(53) 


(XA  • V 


• * + 


(XA  ' V 


y.  - y.  yA/y  2(y,  - y„) 


.,2 


A B 


A B 


CHORD 


Under  certain  conditions  (x^  and/or  x^  < 0),  we  have  to  add  "ACHORD": 


ACHORD 


-2 J xfiAG  dy?  = c/a  y’  ^ a2  - y’2  + a2SIN  * (y'/a) 


y' 

JB 


AO I = BAG  - LINE  + ACHORD 


B 
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We  now  have  the  terms  necessary  for  the  volume  of  intercept  calculation: 

(54)  VOI  = (AOI)W  . 

AVG 


where 

_ ^ ' ^BAG^b  * RHEAD  WH 

AVG  = "ft  ' Vg  + “head 


SIDE  VIEW  TOP  VIEW 


The  error  introduced  by  adding  in  the  volume  marked  * above  only  occurs  for 
AB  > and  will,  for  this  condition,  be  a very  small  percentage  of  the  total 
volune  which  will  be  compensated  somewhat  by  bag  stretch;  i.e.,  we  may 
approximate  the  missing  volume  * by  (we  assume  LINE  AB  perpendicular  to 
x-axis) : 


VMISSING 


W 

2hAVG  * T X ” 
2 


f 


B ~ 13 


where 


h * 

AVG 


Vx  + ^VlAX 


Wh 

c - x'  @ y'  = -j- 


m 


x' 


Zl 


= 1 , 


Since  we  are  interested  in  this  volume  at  the  worst  condition,  high 
penetrations,  let: 


AB  - 20  inches 
a'  = 13  inches 


c’  - 6 inches  (a  bag  26  inches  high  by  12  inches  deep) 

» 15  inches 


^MISSING 


13.7  cubic  inches 


BAG,  ORIG. 


= 4/3tt  a c * 4/3tt  (13)  (6)  = 4244  cubic  inches 


% ERROR  = 


13.7 

4244 


0.31,  a negligible  amount. 


We  therefore  conclude  that  we  may  safely  sidestep  the  computation  of  the 
small  volumes  marked  by  the  *. 
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The  volume  of  the  airbag  for  a given  intercept  AB  is  extremely  difficult  to 
compute  exactly  due  to  the  assymmetry  of  the  volumes  marked  by  + and  - on 
the  previous  page.  In  order  to  facilitate  this  computation  we  make  the 
following  assumption. 

Assumption:  The  bag  volume  for  a given  volume  of  intercept,  VOI,  is  in- 

dependent of  the  torso  inclination  if  the  deformed  periphery  is  the  same; 
i.e.,  the  crosshatched  volumes  are  equal  for  constant  periphery  and  constant 
volume  of  intercept  CVOI) . 


The  validity  of  this  assumption  must  be  checked  by  computer  results  versus 
test  results. 

For  subsequent  volume  computations  we  use  the  actual  volume  of  intercept 
for  a link  in  order  to  change  the  assymmetrical  problem  into  a symmetrical 
one. 


For  the  symmetrical  case  we  assume  a disc  penetrates  the  bag. 
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(58)  VOISYM  = -rra 


(59)  VOISYM  = -rra 


xSYM  - 


xSYM' 


xSYM3 


xSYM 


+ xSYM 


3c 


xSYM2 


3c 


This  cubic  equation  will  be  solved  in  the  computer  program  by  using  the 
Newton - Raphson  method. 

For  constant  periphery; 

— 

(60)  Ellipse  Perimeter  - 2 (OC)  = 2 BGSYM  + ABSYM 


Rewriting  (60) ; 


2 + 2 — 

2 it  \J  a— - 2(OC)  - 2 BCSYM  + ABSYM  where, 


x' 


ABSYM  = 2y'  8 x'  = c - xSYM  in  — - 


y'  = a 


4 


^ (c  - xSYM) 


v*  2 

* \ -l 
a 


(61)  ABSYM  = 2a 


f ^ 


- xSYM) 


Solving  (60)  for  BCSYM 


(62) 


27rt/a-  - 2 OC  - 2a 

BCSYM  = ’ 2 


OC  - 2a  ^ 1 - ~ 


/■ 


With  BCSYM  known  we  can  solve  for  RSYM  and  8SYM;  i.e. , 
(63)  RSYM  (0SYM)  = BCSYM 
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(64)  RSYM  SIN  ! 

(65)  BCSYM  = y (x£  - 


BCSYM 


where , 


x')2  + (7g  - y^)2  (known) 


We  now  have  2 equations  in  2 unknowns  (63)  5 (64) 


C66)  w SIN  ,esm 


BCSYM 


This  transcendental  equation,  like  Equation  25  will  be  solved  numerically 
on  the  computer  for  RSYM  6 0SYM. 


Now  to  calculate  the  bag  volume. 

/ 

X, 


(67)  V 


ABCD 


- , [r 

J c-xS 


dx' 


where. 


c-xSYM  = 


*bs 


(68)  y"'iT^(x'  -xcj  +yc 


Substituting  (68)  into  (67) : 


(69)  V = rr 
^ ABCD 


- yc 


2 x. 


*BS 


- x: 


*bs 

X 


dx’ 


where  is  the  volume  enclosed  by  the  frustum  ABCD. 
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(70)  VBC  = 2ttRS  f y RSYM2  (eSYM  - SIN(0SYM) 

where  is  the  volume  of  the  ring  around  the  volume 

RS  = Distance  to  centroid  of  the  segment  BC  • EC 


(71) 


(72) 


(73) 


(74)  V = £ / ^^substituting  69,  70,  and  74  into  75  will 

fo  \ A V '''yield  the  bag  volume,  V^y-  . 

1 U I/AJL 

(75)  V = V + V + V 

J TOTAL  ABCD  BC  BRIM 

(For  any  AB  resulting  from  a calculated  volume  of  intercept,  VOX ._) 


This  completes  the  derivations  necessary  for  programming  the  bag  forces  and 
geometry. 
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APPENDIX  C 


DERIVATION  OF  THE  STEERING  COLUMN  ALGORITHM 
COLUMN  FORCE  CALCULATIONS 


The  total  force  acting  axially  along  the  column  is  given  by: 

I 

(76)  fac  ■ W 035  (9c  ■ V + fhead  005  (ec  ■ V 

The  total  force  acting  normal  to  the  column  is  given  by: 

(77)  fnc  ’ W SIN  (9c  - V + fhead  sin  (9c  - V 

The  total  moment  acting  at  point  0 is  given  by: 

(78)  M = F^^  • t (t  is  shown  on  page  35  and  calculated  on  page  53. 

U Lrlhol 


C 


1 


CM  Type  Column 


= reaction  at  shear  capsule 
R^  = reaction  at  aft  bushing 
Rgp  = reaction  at  forward  bushing 


Solve  for  external  reactions: 


'FW 


SC 


'SC 


H™t 


NC. 


FW 


The  problem  is  statically  indeterminate;  however,  it  can  be  reduced  to: 


(79)  R__  » FVIO  + 

J SC  NC 


- CM  + F £ ) 

2 U 0 NC  SCJ 


"FW 


C - 2 


(80) 


Av 


;M  + F Z > 
3/0  NC  SC 


2 V Z 


FW 


<81>  % 


M + F Z 
0 NC  SC 


We 


may  now  solve  for  R^p  and 


_ L 3 /Mp  * fnc  £sc\ 

^Af  = MBf‘  ^ UBF  " W " |_FNC  + 2\  Z^  ) 


(£bf  - V 


+ F Z ♦ M = M + F Z 
NC  BF  0 0 NC  SC 


3 / M0  * FNC  £SC 
RgF  NC  2 l 


"bF 


< F... 


FW 


0 NC  SC 


SL 


t«BA  H 


NC 

A 


'SC 


'BA 


'BF 


Solving  for  R^, 


(82^  V 


/ M + F Z 
„ 3/0  NC  SC 

■ NC  2V  lF„ 


(£  - £ ) + F.t_  (£„_  - Znr.) 

BF  SCJ  NC  v SC  BF 


l - J l 
BF  BA 


^Aa  ° = _rbf  Af  " W 


(M  + F Z 
0 NC  SC’ 


^BA  ' £SC-) 


Sic  £ba  VV  fnc  ^sc 


c - 3 


Solving  for  Rg  , 


(83)  v 


F ♦ 3 

f Mo  + fnc  £sc\ 

FNC  2 

V ^ FW  A 

(tBA  ' £SC5  + FNC  (iSC  ' £BA^ 


BF  BA 


Note:  For  a pinned  end  at  the  firewall  the  3/2  factor  in  Equations  79,  80, 

82  and  83  is  equal  to  1.0  and  s 0. 

Solve  for  "t"  the  moment  arm  for  Fp. 


The  methodology  for  this  calculation  is  as  follows: 

1.  Find  equation  for  line  from  x^,y^  perpendicular  to  AB. 

2.  Find  distance  from  this  line  to  x^,y^,.  This  distance  is  "t". 

The  equation  for  line  AB  is  given  by  Equation  10,  i.e., 

B + x'  (m  C0S6W  - SIN6J 
7 * COS0w  + m SING^ 


where  B s Y2  ' - x ). 

The  slope  of  AB  is: 

m CQS6^  - SIN6W 
m'  = COS0W  + m SIN6W 


For  a line  perpendicular  to  this  line, 

“PAB  = " in' 

The  line  from  x^,y^  perpendicular  to  AB  is  given  by 
(y’  ' y0}  * mPAB  (x'  ' V 
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Rewriting  the  equation. 


^AB  X'  * y'  - y0  * ”£AB  Xc 


The  distance  t is  given  by 


0 


t > 0 for  x^,y^,  above  a parallel  line  through  point  P 
t < 0 for  bel°w  a parallel  line  through  point  P 


(page  B-7 ) . 
(page  B-7)  . 
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APPENDIX  D-l : PROGRAM  LISTING 


DF'.ACR 


07/28/8? 


10CC 

120C 

HOC 

HOC 

180C 

200C 

220C 

230C 

231C 

232C 

233C 

234C 

235C 

236C 

E 

237C 
238C 
239C 
240C 
241C 
242C 
300C 
320C 
340C 
360C 
38  OC 
400C 
420C 
440C 
460C 
480C 
500C 
520C 
540C 
560C 
580C 
600C 
620C 
640C 
660C 
680C 
700C 
720C 
740C 
760 
780 
75>. 
800 
820 
840 
360 
870 
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#«««  DRACR  **** 


THIS  PROGRAM  PREDICTS  THE  DRIVER  KINEMATICS  IN  A CRASH  SITUATION 
IN  WHICH  THE  DRIVER  IS  RESTRAINED  BY  AN  AIRBAG  AND  KNEE  RESTRAINT . 
THE  AIRBAG  IS  ATTACHED  TO  A STEERING  COLUMN  AND  WHEEL  THAT  COLLAPSE 
ACCORDING  TO  A PREDETERMINED  FORCE-CRUSH  CHARACTERISTIC. 

IN  ADDITION  TO  THE  SIMPLE  AXIAL  CRUSH  OF  THE  STEERING  WHEEL  AND 
STEERING  COLUMN.  THE  PROGRAM  HAS  THE  CAPABILITY  OF  INDIVIDUALLY 
SIMULATING  BOTH  STEERING  WHEEL  AND  STEERING  COLUMN  ROTATION.  THUS. 
FOUR  SEPARATE  ENERGY  ABSORBING  MODES  ARE  POSSIBLE  FOR  THE  WHEEL 
AND  COLUMN.  WHICH  ONE  OF  THESE  MODES  (IE  ANY)  WILL  ACTUALLY  OCCUR 
IN  A GIVEN  CRASH  SITUATION  DEPENDS  UPON  TWO  CONDITIONS.  FIRST.  THE 
FORCE  OR  MOMENT  MUST  BE  LARGE  ENOUGH  - THAT  IS  GREATE  R THAN  I HE  FORC 

OR  MOMENT  CAPABLE  OF  BEING  RESISTED  BY  THE  WHEEL  OR  COLUMN.  SECOND. 
IF  THE  FORCE  OR  MOMEMT  IS  LARGE  ENOUGH  FOR  PERMANENT  DEFORM- 
ATION IN  MORE  THAN  ONE  OF  THESE  FOUR  MODES  TO  OCCUR.  THE  PROGRAM 
USES  A MINIMUM  ENERGY  PRINCIPLE  TO  SELECT  JUST  WHICH  DEFORMATION 
MECHANISM  ACTUALLY  WILL  OCCUR.  THE  ONE  OF  THE  FOUR  REQUIRING  THE 
LEAST  ENERGY  IS  THE  ONE  SELECTED. 

THE  DRIVER  IS  MODELED  BY  THREE  MASSES-A  HEAD  MASS.  A TORSO  MASS  AND 
A LOWER  BODY  MASS.  THE  DRIVER  IS  CONSTRAINED  TO  HAVE  PLANAR  MOTION 
SO  THAT  THE  PROGRAM  IS  STRICTLY  APPLICABLE  ONLY  TO  FRONTAL  CRASH 
SITUATIONS. 

EVALUATIONS  OF  AIRBAG.  STEERING  WHEEL.  STEERING  COLUMN.  KNEE 
RESTRAINT  AND  VEHICLE  PERFORMANCE  CAN  BE  MADE  BY  APPROPRIATE  CHANGES 
IN  THE  DESIGN  PARAMETERS. 

TYPICAL  DESIGN  PARAMETERS  THAT  CAN  BE  EVALUATED  ARE  BAG  SIZE.  BAG 
SHAPE.  INFLATION  CHARACTERISTICS.  VENT  AREA.  STEERING  COLUMN  AND/OR 
STEERING  WHEEL  CRUSH  CHARACTERISTICS.  KNEE  RESTRAINT  CRUSH  CHAR- 
ACTERISTICS. STEERING  COLUMN  AND  WHEEL  SUPPORT  STRUCTURE  STIFFNESS. 
AS  WELL  AS  OTHER  SYSTEM  PARAMETERS. 

THIS  PROGRAM  IS  SELF  CONTAINED  IN  THAT  NO  EXTERNAL  FUNCTIONS  OR 
SUBROUTINES  ARE  REQUIRED. 

AUTHOR:  MICHAEL  FITZPATRICK 

FITZPATRICK  ENGINEERING 
WARSAW.  INDIANA  46580 
TEL.  (?19>-267-4437 
DEC . 1 . 1981 


FILENAME  INFILE 

CQMMON/OUT/NPD.T ( 175) . X0(6. 175) .XI (6. 175) . X2( 6. 175 ) ,X3(6, 1 75) . X4<  6 . 1 

1X5(6.175) ,X6( 6. 175) ,X7( 6,175 ) ,X8' 6.175 >.X10( 6. 175) 

COMMON/OUT 1/N9 . X9 ( 4 . 175 ) . T9 ( 1 75 ) 

COMMON/NAME/INF ILE 

COMMON/FLAG/NCOLTYPE 

COMMON/TIME /STEP. XST0P.PINT1.IL00P 


DRACR 


07/28/82 


871 

C0MM0N/HIC/THIC(  175)  ,HRGS(  175) 

880 

1050  PRINT  .“INPUT  FILE  NAME" 

900 

INPUT  .INFILE 

920 

PRINT. “ENTER  1 IF  COLUMN  IS  SUPPORTED 

A!'  THE  FIREWALL:  OTHERWISE 

940 

% ENTER  2“ 

960 

INPUT. NCOLTYPE 

980 

NPD=0 

1000 

N9=0 

1020 

CALL  SQLUE(8) 

1040 

IF  < NPD.GT . 175)NPD=175 

1060 

IF ( N9.GT . 175 )N9=175 

1080 

1120  FORMAT ( 1H- ) 

1100 

1125  FORMAT (F7 .4 , 2F6. 1 .2F7. 2.2F8. 3.F8. 1 .2F7. 2) 

1120 

1130  FORMAT ( IX ,7F1 1 .2) 

1140 

1140  FORMAT (0) 

1160 

1150  FORMAT ( IX. 7(4X. "=======")//) 

1180 

PRINT  1120 

1200 

1170  FORMAT ( IX."  TIME  OEM  G’S 

MEH  MEL 

OEH  DISP 

BODY 

1220 

% G'S  COL  DISP  BAG  PRESS“/1X.“ 

(MS) 

(G’S) 

(MPH 

1240 

1 (INCHES)  (G’S)  (INCHES) 

(PS1G)“) 

1260 

PRINT  1170 

1280 

PRINT  1150 

1300 

DO  1221  K=1 .NPD 

1320 

1221  PRINT  1130. T(K) , (X0( J.K ) , J=1 . 6 ) 

1340 

PRINT  1120 

1360 

1223  FORMAT ( IX. “ TIME  FEM  ANGLE 

H-P  MEL 

H-P  ACC 

FEM  FO 

1380 

&RCE  SEAT  FR.  H-P  R.D.“/1X.“ 

(MS) 

(DEG) 

(MPH' 

1400 

4 (G’S)  (LBS)  (LBS) 

(INCHES) 

“) 

1420 

PRINT  1223 

1440 

PRINT  1150 

1460 

DO  1230  K=1.NPD 

1480 

1230  PRINT  1130.T(K),(X1(J.K).J=1.6) 

1500 

PRINT  1120 

1520 

1250  FORMAT ( IX,"  TIME  TORSO  DISP 

TORSO  ANG 

TORSO  MEL 

TORSO 

1540 

l ACC  TORSO  R.D.  TORSO  R.U.“/1X.“ 

(MS) 

(INCHES) 

1560 

% (DEG)  (D/SEC)  (D/SEC*#2) 

(INCHES) 

( MPH ) 1 ' ) 

1580 

PRINT  1250 

1600 

PRINT  1150 

1620 

DO  1310  K=1 .NPD 

1640 

1310  PRINT  1130.T(K).(X2(J.K).J=1.6> 

1660 

PRINT  1120 

1680 

1330  FORMAT (IX."  TIME  HEAD  DISP 

HEAD  ANG 

HEAD  MEL 

HEAD 

1700 

l ACC  HEAD  R.D.  HEAD  R.ANG'VIX." 

(MS) 

(INCHES) 

1720 

& (DEG)  (D/SEC)  (D/SEC**2) 

( INCHES) 

(DEG)") 

1740 

PRINT  1330 

1760 

PRINT  1150 

1780 

DO  1380  K=1 .NPD 

1800 

1380  PRINT  1 130.T(K) , (X3( J.K) , J=1 .6) 

1820 

PRINT  1120 

1840 

1390  FORMAT ( IX. “ TIME  WH  AX  FOR 

WH  N FOR 

WH  MOMENT 

WH 

DRACR 

1860 

188Q 

1900 

1920 

1940 

1960 

1980 

2000 

2020 

2040 

2060 

2080 

2100 

2120 

2140 

2160 

2180 

2200 

2220 

2240 

2260 

2280 

2300 

2320 

2340 

2360 

2380 

2400 

2420 

2440 

2460 

2480 

2500 

2520 

2540 

2560 

2580 

2600 

2620 

2640 

2660 

2680 

2700 

2720 

2740 

2760 

2780 

2800 

2820 

2840 
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5 RESIST  WH  STROKE  WH  ST  0EL"/1X."  (MS)  (LBS) 

6 (LBS)  (IN-LBS)  (LBS)  (INCHES)  (IN/SEC)") 

PRINT  1390 

PRINT  1150 
DO  1395  K=1.NPD 

1395  PRINT  1130.T(K).(X4(J,K).J=1.6) 

PRINT  1120 

1400  FORMAT (IX."  TIME  COL  AX  FOR  COL  N FOR  COL  MOMENT  COL 

l RESIST  COL  STROKE  COL  ST  VEL"/1X."  (MS)  (LBS) 

X (LBS)  (IN-LBS)  (LBS)  (INCHES)  (IN/SEC)") 

PRINT  1400 
PRINT  1150 
00  1460  K=1 „NPD 

1460  PRINT  1130.T (K) . (X5( J.K) .0=1 .6) 

PRINT  1120 

1465  FORMAT (IX."  TIME  WH  AP  MOM  WH  RES  MOM  WH  ANGLE  WH 

& ANG  OEL'VIX."  (MS)  (IN-LBS)  (IN-LBS)  (BEG) 

& (DEG/SEC)") 

PRINT  1465 
PRINT  1150 
DO  1470  K=1.NPD 

1470  PRINT  1130.T(K) . (X6( J.K) , J=1 .4) 

PRINT  1120 

1473  FORMAT (IX."  TIME  COL  AP  MOM  COL  RES  MOM  COL  ANGLE  COL 

% ANG  VEL'VIX."  (MS)  (IN-LBS)  (IN-LBS)  (BEG) 

5 ( DEG/SEC >") 

PRINT  1473 
PRINT  1150 

DO  1477  K=1 .NPD 

1477  PRINT  1130.T (K) , <X7( J.K) , J=1 .4) 

PRINT  1120 

1480  FORMAT ( IX."  TIME  BAG  PEN.  BAG  OQL.  BAG  PRESS.  W/A 

% FORCE  P.  F0RCE"/1X."  (MS)  (INCHES)  (CU.IN.) 

6 (PSI.G)  (LBS)  (LBS)") 

PRINT  1480 

PRINT  1150 
DO  1540  K=1  .NF'D 

1540  PRINT  1130,T(K) . (X8( J.K) , J-l ,5) 

PRINT  1120 

1560  FORMAT (IX."  TIME  CHEST  AP  CHEST  SI  HEAD  AP 

%.  HEAD  SI"/1X."  (MS)  (G'S)  (G'S)  (C’S) 

& (G'S)") 

PRINT  1560 
PRINT  1150 
DO  1620  K=1.N9 

1620  PRINT  1130. T9(K> . <X9(  J.K) . J=.t . 4)  .-PRINT  1120 

162?  FORMAT (IX."  TIME  EAWC  LACC  EAWR  EACR 

l EA'  '/ 1 y " (MS)  (LBS)  (LBS)  (LBS' 

% (LBS)  (LBS)  ") 

PRINT  1622 


DRACR 

2860 

2880 

2900 

2920 

2940 

2960 

2980 

2982 

2983 

2984 

2986 

2988 

2990 

2992 

2994 

2996 

2998 

3000 

3002 

3004 

3006 

3008 

3010 

3012 

3014 

3016 

3018 

3020 

3022 

3030 

3035 

3040C 

3060C 

3080C 

3100 

3120 

3140 

3160 

3180 

3200 

3220 

3240 

3260 

3280 

3300 

3320 

3340 

3360 

3380 

3400 
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PRINT  1150 
DO  1625  K=1.NPL» 

1625  PRINT  1 1 30 „ T ( K ) . ( X 1 0 ( J . K ) , J= 1 , 5 ) 

PRINT  1120 

1630  PRINT. "ENTER  1 TO  CALCULATE  HiC" 

INPUT  . NRES 

IF (NRES.NE.l ) GO  TO  2000 
1640  PEAK=0. 

NSTQP=N9 

DO  1715  1=1 .NSTQP 
DO  1716  J=1.I 
L=I+1 
SUM=0. 

DO  1717  K=i . J 
L=L-1 

SUM=SUM+HRGS(L)*PINT1 

1717  CONTINUE 
DELT=THIC(K> 

CHECK=SUM/DELT 

IF (PEAK-CHECK)  1718.1716.1716 

1718  PEAK=CHECK 
TL0W=(L-1)*PINT1 
THIGH=1«PINT1 

1716  CONTINUE 
1715  CONTINUE 

HIC=PEAK**2.5 
PRINT. "THE  HIC  IS". HIC 
PRINT. "T1=".TL0W 
PRINT. "T2=". THIGH 
2000  STOP 
END 

THIS  SUBROUTINE  SETS  UP  THE  DIFFERENTIAL  EQUATIONS  THAT  DESCRIBE 
THE  DRIVER  KINEMATICS. 

SUBROUTINE  DIFEQ(T, Y.DY) 

CGMMQN/MANDAT/ZL . ZT , ZH , RT . RN . RH . RTUPH . X2Z . Y2Z . WB . WH 
DOUBLE  PRECISION  Y(8) 

DIMENSION  DY(8) 

CALL  FORCETH( Y.TNECK) 

CALL  DECEL(T.GS) 

CALL  BAGSUB (T. Y. TNECK .FTH.FX.FTT ,GS) 

SH=SIN(Y(6) ) 

ST=SIN( Y ( 7 ) ) 

CH=COS( Y ( 6 ) ) 

CT=COS( Y ( 7) ) 

DY < 1 ) = (FX-(ZT#RT+ZH*RN )#(CT*DY (3)-ST»Y(3)*Y(3>) 

&-ZH«RH« (CH*DY<2 )~SH«Y ( 2 )«Y (2  > ) )/( ZL+ZT+ZH ) 

DY ( 2>=(FTH-ZH*RH#DY ( 1 )*CH-ZH*RN«RH« (CT*CH#DY ( 3)+CT*SH#Y(3) 

S*Y (3)-ST#CH*Y ( 3 )*Y(3)+ST#SH#DY<3) )+ZH*32. 17*SH*RH ) 

&/(ZH*RH»RH) 
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3420 

3440 

3460 

3480 

3500 

3520 

3540 

3560 

3580 

3600 

3620C 

3640C 

3660C 

3680 

3700 

F 

3720 

3730 

3740 

3750 

3760 

3780 

3800 

3820 

3840 

3850 

3860 

3880 

3900 

3920 

3940 

3960 

3980 

4000 

4020 

4040 

4060 

4080 

4100 

4120 

4140 

4160 

4180 

4200 

4220 

4240 

4260 

4280 

4300 

4320 

4340 
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DY(3)=(FTT-(ZT*RT+ZH*RN)#DY(1)*CT-ZH*RN*RH*(CT*CH*DY(2)+Y<2) 
fc*Y(2)»(ST«CH-CT#SH)+ST*SH*DY(2) )+ZT*32. 17*RT*ST+ZH#32. 17*RN#ST ) 
t/(ZT#RT*RT+ZH«RN«RN) 

DY(4)=-GS 

DY(5)=Y(1) 

DY(6)-Y(2) 

DY(7>=Y<3) 

DY(8)=Y<4) 

RETURN 

END 

THIS  SUBROUTINE  READS  IN  THE  INPUT  DATA.  SETS  UP  THE  INPUT  DATA 
FOR  DISPLAY  AND  INITIALIZES  KEY  VARIABLES. 

SUBROUTINE  SETUP(X.Y) 

REAL  MU.MC.LSC.LFW.LBA.LBF.LSCZ.LFWZ.LBAZ.LBFZ.KRN.MIP.MWH.MFW.LT.L 


CQMMON/WROTP/NPWRP , WHRP (2.12), OHWO , THETAWO 
COMMQN/WROTN/NPWRN . WHRN (2.12) 

COMMON/CROTP/NPCRP . COLRP (2.12). OMCO . THET ACO 
COMMON/CROTN/NPCRN , COLRN (2.12) 

COMMON/SEATFRIC/NSF , SUN .SFUT . RELSF . SFN ( 2 , 12 ) 

COMMON/KNEEREST/NKR , SKR , RUT , RELKR . KRN ( 2 . 1 2 ) 

CQMMQN/NAME/ INFILE.OUTFILE 
FILENAME  INFILE.OUTFILE 

CQHHON/MANDAT/ZL . ZT . ZH , RT . RN . RH . RTOPH ,X2Z. Y2Z.WB.WH 
COMMON/MMANDAT/LT , LF . THFO . THLO . DELTTB . THEF 
COMMQN/NECK/NPN , FNECK ( 2 . 1 2 ) . DCN 
COMMON/ VEH/NV . VEHGS (2,30) 

COMMON/G ASFLO/NPG . GEN ( 2 . 1 2 ) 

COHMQN/WHFOR/NPW. WHE (2.12). VWHCQ . SWHCO 
CQMMON/COLFOR/NPC . COL ( 2 , 1 2 ) 

COMMON/GASDAT/ ATMOP .PGZ.GTZ.U.PNl . PN2 . PN3 
CQMMON/BAGDAT/VC 1 . VC2 . AV , SA . SC . X 1 . Y 1 

COMMON/COLDAT/THETAC . THET AC Z . MU . LSCZ . LF WZ . LBAZ .LBFZ.WC.RCQL 
COMMON/WHEEL/XWH .YWH.YO .RIMRAD . RIP . WWH . WIP . M I P . MWH 
COMMON/PARAM/RFT . THETATZ . THETAHZ . F AW . FNW . FRW . SWHC . VWHC . T AW . TRW 
COMMON/MMPAR AM/BP . VOL . FFT . FP . THETAW . THETAWZ .OMW.TAC.TRC. RBRKT . OMC 
COMMON/MI SC/VCOLCO , SCOLCO .PR8.BIG 

COMMON/MM ISC/FM2 , PA5 , GT . FPN , MC . VOLO .GW.LSC.LFW, LBA . LBF 
DOUBLE  PRECISION  Y(8) 


x=o. 

2070  FORMAT (V) 

2080  FORMAT (IX, "INITIAL  VELOCITY:  ".G10.3/1X, 

^''INITIAL  HEAD  ANGLE:  ".G10.3/1X. "INITIAL  TORSO  ANGLE: 

1G10.3) 

2100  FORMAT (IX. 

&"  MLEG  MTORSO  MHEAD  RT  RN 

l RT0PH"/1X. 7G10. 3) 

2120  FORMAT (IX,"  ATMOP  PGZ  GTZ  U 

& PN2  PN3"/1X,7G10.3) 

2122  FORMATdX."  THETAW  RIP  WWH  WIP 


I I 


RH 


PN1 

RCOL 
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l DCN"/1X.6G10.3) 


2130  FORMAT (IX." 

VC1 

UC2 

AV 

SA 

SC 

rvj 

yH 

X 

YIZ'VIX,  7G10.3) 

2132  FORMAT (IX." 

THETAC 

MU 

LSCZ 

LFWZ 

& LBAZ 

LBFZ 

UC'VIX, 

8G10.3) 

2133  FORMAT (IX," 

LF 

THFO 

THLO" 

&/1X.5G10.3) 

2134  FORMAT (IX," 

XWH 

RIMRAD 

X2Z 

Y2Z 

t WB 

WH"/1X 

.7G10.3) 

2135  FORMAT* IX, "GAS 

FLOW  TIME 

- MSEC"/ IX 

, 10G10.3) 

2136  FORMAT (IX. "GAS  FLOW  - LB/SEC'VIX.  10G10.3) 

2137  FORMAT* IX. "COLUMN  STROKE  - INCHES"/1X. 10G10.3) 

2138  FORMAT (IX. "COLUMN  FORCE  - LBS'VIX. 10G10.3) 

2140  FORMAT (IX. "NECK  ANGLE  - DEG'VIX. 10G10.3) 

2150  FORMAT (IX. "NECK  TORQUE  - FT-LBS"/1X. 10G10.3) 

2160  FORMAT ( IX,"  NPTS  NECK  NPTS  KR  NPTS  UEH  NPTS  SEAT  NPTS  GAS 
% NPTS  COL  NPTS  WHC71X.  10G10.3) 

2170  FORMAT (IX."  NP  WRP  NP  WRN  NP  CRP  NP  CRN'VIX. 10G10.3) 
2180  FORMAT ( IX ."  SL  SEAT  SL  KR"/1X.10G10.3) 

2250  FORMAT ( 1X,"VEH.  PULSE  TIME  - MSEC71X.  10G10.3) 

2260  FORMAT* IX. "VEH.  PULSE  DECELERATION  - G'S'VIX.  10G10.3) 

2265  FORMAT ( IX. "SEAT  FRICTION  DISPLACEMENT  - INCHES'VIX. 10G10.3) 

2267  FORMAT (IX, "SEAT  FRICTION  FORCE  - LBS"/1X. 10G10.3) 

2268  FORMAT ( IX. "KNEE  DISPLACEMENT  - INCHES'VIX. 10G10. 3) 

2269  FORMAT* IX. "KNEE  FORCE  - LBS'VIX, 10G10.3) 

2280  FORMAT* IX, "WHEEL  STROKE  - INCHES’VIX. 10G10. 3 ) 

2290  FORMAT* IX. "WHEEL  FORCE  - LBS'VIX. 10G10. 3) 

2300  FORMAT (IX, "WHEEL  ANGLE.  POS.  - DEG  "/IX. 10G10.3) 

2310  FORMAT ( IX, "WHEEL  TORQUE.  POS.  - IN-LBS' V1X, 10G10.3) 

2314  FORMAT (IX. "WHEEL  ANGLE.  NEG.  - DEG'VIX. 10G10. 3) 

2318  FORMAT ( IX ."WHEEL  TORQUE.  NEG.  - IN-LBS"/1X. 10G10.3) 

2320  FORMAT (IX. "COLUMN  ANGLE.  POS.  - DEG'VIX. 10G 10. 3) 

2330  FORMAT ( IX. "COLUMN  REACTION.  POS.  - IN-LBS  IF  2:  LBS  IF  1"/1X.10G10. 

2334  FORMAT* IX. "COLUMN  ANGLE.  NEG.  - DEG'VIX. 10G10. 3) 

2338  FORMAT* IX. "COLUMN  REACTION.  NEG.  - IN-LBS  IF  2:  LBS  IF  1"/1X.10G10. 

READ* INFILE,2070)Y (4).Y(6),Y(7) 

READ (INFILE. 2070 ) ZL , ZT . ZH . RT . RN . RH . RTOPH 

READ (INFILE. 2070  > NPN .NKR.NV. NSF . NPG . NPC . NPW . NPWRP . NPWRN . NPCRP . NPCRN 

READ(INFILE»2070)SUN,SKR 

READ* INFILE. 2070) (GEN(1.K),K=1. NPG) 

READ (INFILE. 2070) (GEN* 2, K).K=1. NPG) 

READ ( INFILE , 2070 ) ( COL ( 1 . K ) , K=1 . NPC ) 

READ (INFILE. 2070 ) ( COL ( 2 „ K ) , K= 1 . NPC ) 

READ( INFILE, 2070)ATM0P.PGZ.GTZ.U,PN1 .PN2.PN3 

READ ( I NF I LE . 2070 ) THET AW , R IP , WWH . W IP . RCOL . DCN 

READ( INFILE. 2070)VC1.VC2. AO. SA. SC. X1Z.Y1Z 

READ (INFILE, 2070 ) THETAC , MU , L SCZ , LFWZ . LBAZ . LBFZ . WC 

READ ( INF I LE . 2070 ) LF , THFO . THLO 

READ ( INFILE . 2070 ) XWH . RIMRAD . X2Z . Y2Z . WB . WH 

READ( INFILE. 2070) (SFN(1.K).K=1. NSF) 
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READ ( INFILE . 2070 ) ( SFN ( 2 . K ) . K= 1 , NSF ) 

READ (INFILE, 2070 ) (FNECK (l.K) ,K=1,NPN) 

READ < I NFI LE , 2070 ) < FNECK  < 2 . K > , K= 1 . NPN  > 

READ<  XNFILE.2070) ( VEHGS( 1 ,K) ,K=1 .NO) 

READ( INFILE, 2070) (VEHGS(2.K> ,K=1 , NO) 

READ( INFILE, 2070) (KRN(l.K) ,K=1.NKR) 

READ! INFILE. 2070) <KRN(2.K).K=1, NKR) 

READ( INFILE, 2070) (WHEQ  ,K>  ,K=1  .NPW) 

READ< INFILE. 2070) <UHE< 2. K).K=1.NPW> 

READ ( INFILE . 2070 ) ( WHRP (l.K) ,K=1 . NPWRP ) 

READ (INFILE, 2070 ) ( WHRP ( 2 . K ) . K= 1 . NPWRP ) 

READ (I NF I LE , 2070 ) ( WHRN (l.K) ,K=1 , NPWRN ) 

READ ( INF I LE , 2070 ) ( WHRN ( 2 , K ) . K= 1 , NPWRN ) 

READ ( I NF I LE , 2070 ) ( COLRP (l.K) ,K=1 . NPCRP ) 

READ (INFILE , 2070 ) ( COLRP  ( 2 . K ) . K=1 . NPCRP  > 

READ (INFILE, 2070 ) ( COLRN ( 1 . K ) , K= 1 , NPCRN ) 

READ( INFILE. 2070) (C0LRN(2,K) .K=l .NPCRN) 

PRINT  2490 

PRINT  ."INPUT  VALUES  — INPUT  UNITS(  MSEC.  MPH.  DEGREES. 
& INCHES,  LBS,  FT-LBS.  G'S)" 

2480  FORMAT ( 1X.10G10.3) 

2490  FORMAT (1H-) 

2520  PRINT  2080, Y(4) ,Y(6) ,Y(7) 

PRINT  2100 . ZL . ZT , ZH , RT . RN , RH . RTOPH 
PRINT  2 1 60 , NPN . NKR , NM , NSF , WPG . NPC . NP W 
PRINT  2170. NPWRP . NPWRN . NPCRP . NPCRN 
PRINT  2180, SUN, SKR 
PRINT  2135, (GEN (l.K) ,K=1,NPG) 

PRINT  2136, (GEN(2,K) ,K=1,NPG) 

PRINT  2137, (C0L(1,K),K=1, NPC) 

PRINT  2138, <C0L(2,K),K=1, NPC) 

PRINT  2265, (SFN(1 ,K) ,K=1 r NSF) 

PRINT  2267, (SFN(2,K) ,K=1 , NSF) 

PRINT  2140, <FNECK(1,K).K=1. NPN) 

PRINT  2150, (FMECK(2,K) ,K=1 „NPN) 

PRINT  2250 , ( VEHGS ( 1 , K ) , K= 1 . W ) 

PRINT  2260 , ( VEHGS ( 2 , K ) , K= 1 . NV ) 

PRINT  2268 , (KRN(i.K) ,K=1. NKR) 

PRINT  2269, (KRN(2,K) ,K=1 . NKR) 

PRINT  2280 , ( WHE ( 1 . K ) , K= 1 , NPW ) 

PRINT  2290. <WHE(2,K),K=1, NPW) 

PRINT  2300 , ( WHRP (l.K) ,K=1 . NPWRP ) 

PRINT  2310. ( WHRP ( 2 . K ) . K- 1 , NPWRP ) 

PRINT  2314, (WHRN(1.K),K*1. NPWRN) 

PRINT  23 1 8 . ( WHRN ( 2 . K ) . K= 1 . NPWRN ) 

PRINT  2320 , (COLRP (l.K) ,K=1 , NPCRP ) 

PRINT  2330 , ( COLRP (2.K) „K=1 .NPCRP ) 

PRINT  2334, ( COLRN ( l.K). K=l. NPCRN) 

PR I NT  2338 . ( COLRN ( 2 , K ) . K= 1 . NPCRN ) 

PRINT  2120.ATM0P.PGZ .GTZ.U.PNl .PN2.PN3 
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PRINT  2130.VC1 ,0C2, AV.SA.SC ,X1Z. Y1Z 

PRINT  2132 , THETAC . MU . LSCZ.LFWZ . LBAZ . LUF Z . WC 

PRINT  2i33eLF,THF0,THL0 

PRINT  2134 . XWH , RIMRAD . X2Z , Y2Z . W8 . WH 

PRINT  2122, THETAW. RIP. UUH.WIP.RCOL.DCN 

Y(2)=0. 

Y(3)=0. 

Y( 4 )=Y<  4 ) #1 .4666667 
Y(5>=0. 

Y(6)=Y(6)« ,01745329 
THETAHZ=Y(6) 

Y(7)=Y<7>*. 01745329 
THETATZ=Y<7> 

Y(8)=0» 

Y( 1 )=Y(4) 

THFO=THFO«. 01745329 

THLO=THLO«. 01745329 

LT=<Y2Z+LF*SIN(THF0) )/SIN(THLO) 

DELTTB=(X2Z-LT#C0S<THL0)-LF«C0S<THF0>>»C0S<THL0> 

ZLsZL/32. 17 

ZT =ZT/32 , 17 

ZH=ZH/32. 17 

RT=RT/12. 

RN=RN/12. 

RH=RH/12. 

RT0PH=RT0PH/12. 

SUN=SUN#12. 

SKR=SKR*12. 

TH£TAW=THETAW# .01745329 
THETAC=THETAC* .01745329 
ADD=SC«SQRT ( 1 . -RIMRAD**2/SA**2 ) 

X1=X1Z+ADD#CQS( THETAW) 

Y1=Y1Z+ADD#SIN< THETAW) 

0WHC0=0. 

SWHC0=0. 

VCOLCOsO„ 

SCOLCOO. 

0MH0=0. 

THETAWZ=THETAW 
THETAMO=THET  AW 
0MCQ=0. 

THETACZ=THETAC 

THETACO=THETAC 

PR8=<2./(PN1+1.))#*(PN1/(PN1-1„)) 

FM2=0C2*SQRT ( PR8*« ( 2 . /PN1 ) -PR8#* ( (PN1+ 1 . ) /PN1 ) ) 

PA=PGZ+ATMOP 

PA5=PA 

GT=GTZ 

FPN=PN2 

HWH=WWH/386. 


DRACR 

6960 

6980 

7000 

7020 

7040 

7060 

7080 

7100 

7120 

7180 

7200 

7220 

7240 

7250 

7260 

7270 

7280 

7300 

7320 

7340 

7360 

7380 

7400 

7460 

7470 

7480 

7540 

7560 

7580 

7600C 

7620C 

7640C 

7660 

7680 

7700 

7720 

7740 

7760 

7780 

7800 

7820 

7840C 

7860C 

7880C 

7900 

7920 

7940 

7960 

7970 

7980 


07/28/82 


- 9 - 


MlP=WIP/386„ 

MC=HC/386. 

0QLZ=4 . /3 . *3 . 14159#SA**2*SC 
V0LQ=0QLZ 

GW=  ( PA*0(3LZ ) / ( U«GTZ ) 

LSOLSCZ 

LFW=LFWZ 

LBA=LBAZ 

LBF=LBFZ 

DO  3010  J=1 , 12 

GEN! 1 . J)~GEN( 1. JJ/IOOO. 

FNECK(1.J)=FNECK<1.J>*. 01745329 
WHRP (1 , J > =WHRP <1 . J ) « . 01745329 
WHRN ( 1 . J ) =WHRN < 1 , J > * .01745329 
COLRP ( 1 , J > =CQIRP ( 1 , J ) * .01745329 
COLRN (1, J) -COLRN < 1 . J > * .01745329 
SFN<1» J)-SFN(1, J)/12. 

3010  KRN ( 1 , J > =KRH ( 1 . J ) / 1 2 . 

DO  3050  J=1.30 

UEHGS ( 1 » J ) =0EHGS ( 1 , J ) / 1 000 . 

3050  UEHGS (2, J) =UEHGS< 2 . J ) *32 . 1 7 
SFUT-O. 

RUT=0. 

Y0=0. 

YWH=0. 

BIG=1000000. 

CONTINUE 

RETURN 

END 

THIS  SUBROUTINE  IS  A GENERALIZED  TABLE  LOOKUP  AND  INTERPOLATION 
ROUTINE  WHICH  IS  CALLED  BY  OTHER  ROUTINES. 

SUBROUTINE  LOOKUP! A. FUN . NPTS . B ) 

D I HENS ION  FUN (2. 30) 

DO  3190  J=1.NPTS 
3190  IF ( FUN ( 1 , J ) . GT . A ) GGT03200 
3200  IF(J.EQ.1)J=2 
K=J~1 

B=( A-FUN ( 1 ,K) >* (FUN( 2. J)-FUN(2, K) )/(FUN( 1 . J)-FUN( 1 ,K) )+FUN( 2,K ) 

RETURN 

END 

THIS  SUBROUTINE  CALCULATES  THE  NECK  TORQUE  AS  A FUNCTION  OF  THE 
NECK  ANGLE . 

SUBROUTINE  FORCETH! Y.TNECK) 

COliHON/NECK/NPN . FNECK  ( 2 . 12 ) . DCN 
DOUBLE  PRECISION  Y<8> 

TNECK-O. 

TDANP*-DCN*(Y(2)-Y(3) ) 

0REL=Y(2)-Y(3) 
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TREL=Y(6)-Y<7> 

IF(TREL.GT„0.0oAND. VREL.LT . 0. )GQ  TO  10 
IF(TREL.LT„0.0. AND.VREL.GT.O. )GQ  TO  10 
CALL  LOOKUP (TREL .FNECK , NPN . TNECK ) 

10  TNECK=TNECK+TDAHP 
RETURN 
END 

THIS  SUBROUTINE  OBTAINS  THE  CRASH  PULSE  G'S  AS  A FUNCTION  OF  TIME. 
SUBROUTINE  DECELCT.GS) 

CQHMGN/VEH/MV . VEHGS (2.30) 

CALL  LOOKUP (T. VEHGS. NO „GS) 

RETURN 

END 

THIS  SUBROUTINE  COMPRISES  THE  MAJOR  PART  OF  THE  DRACR  PROGRAM.  IT 
EVALUATES  THE  BAG  SHAPE  AS  A FUNCTION  OF  BAG  PENETRATION  AND  TORSO 
ANGLE,  CALCULATES  THE  FORCES  THE  BAG  APPLIES  TO  THE  DRIVER.  CALC- 
ULATES THE  BAG  VOLUME  AND  PRESSURE,  DETERMINES  THE  GAS  GENERATOR 
FLOW  CHARACTERISTICS  AND  CALCULATES  THE  STEERING  COLUMN  AND  WHEEL 
FORCES.  MOMENTS,  CRUSH  AND  ROTATION. 

SUBROUTINE  BAGSUB (X.Y, TNECK , FTH , FX ,FTT , GS ) 

REAL  HFW.MU.LSC  LFW.LBA.LBF.LFWZ.LSCZ.LBAZ.LBFZ.MO.MC.MIP.MWH 

REAL  KRN.LT.LF 

COMMON/FLAG/NCOLTYPE 

CQMMON/yRQTP/NPURP . WHRP (2,12), OMWO . THET AWO 
COHMON/WROTN/NPWRN . HHRN (2,12) 

COMHON/CROTP/NPCRP . COLRP (2,12). OMCO . THET ACO 
COMMON/CROTN/NPCRN f COLRN (2,12) 

COMMON/SEATFRIC/NSF , SUN . SFUT , RELSF , SFN ( 2 . 1 2 ) 

CQHHON/MANDAT/ZL , ZT , ZH , RT . RN . RH , RTOPH . X2Z . Y2Z . WB . WH 
COHMON/MMANDAT/LT , LF , THFO , THLO , DELTTB . THEF 
CONHON/KNEEREST/NKR . SKR . RUT . RELKR . KRN( 2 . 1 2 ) 

CGMMON/GASFLQ/NPG . GEN ( 2 . 1 2 ) 

COMMON/WHFOR/NPW . WHE ( 2 , 1 2 ) . VWHCO , S WHCO 
COHMON/COLFOR/NPC . COL ( 2 . 1 2 ) 

CQMMON/GASDAT/ATMOP . PGZ , GTZ . U . PN 1 . PN2 . PN3 
C0MM0N/BAGDAT/VC1 , VC 2 , AV . SA . SC . XI , Y1 

COMMON/COLDAT/THETAC , THET ACZ . MU , LSC  Z , LF WZ . LBAZ .LBFZ.WC. RCOL 

COMMON/ WHEEL/ XWH , YUH , YO , R I MRAD , R IP , WWH , W IP . M I P . MWH 
COMMON/M ISC/VCOLCO, SCOLCO.PR8, BIG 

C0MM0N/MMISC/FM2 , PA5 , GT , FPN , MC . VOLO , GW , LSC . LF W , L BA , LBF 
COMHON/PARAM/RFT , THET ATZ , THETAHZ , FAW .FNW.FRW, SWHC . VWHC . T AW . TRW 
COHMQN/TXME/STEP . XSTOP .PINT1 , ILOOP 

COHMON/MPARAM/FACQL , SCOLC , PG1 , FKNEE , SF , FNCOL . MO .FRCOL . VCOLC . REACT 
COMMON/MMPARAM/BP , VOL . FFT . FP . THET AW . THETAWZ . OMW . TAC . TRC . RRRKT . OMC 
COMMON/EABS/EAWC . EACC , EAWR , E ACR . EA 
DOUBLE  PRECISION  Y(8) .B,A.D.E,A1.B1 .C1.X2.Y2 
5 FORMAT ( 1H-) 

DIMENSION  DY(8) 
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WBACT=WB 

WHACT=WH 

THETAT=Y(7) 

THETAH=Y<6> 

CHECK  TO  SEE  IF  DRIVER  SUBMARINING. 

IF(ABS(THETAW-THETAT) .GT.l .4)  GO  TO  500 
CALCULATE  THE  SLOPE  OF  THE  DRIVER  TORSO. 

0SL0PE=TAN(3. 14159/2. +THETAT) 

CALCULATE  THE  NEW  H-POINT  COORDINATES  AND  THE  X-COORDINATE  OF  THE 
POINT  WHERE  THE  RIM  INTERSECTS  THE  BAG. 

X2=X2Z-<Y<5)-Y<8))*12. 

Y2=Y2Z 

XC=-SC#SQRT ( 1 . -R1MRAD**2/SA**2> 

B<=Y2-Y1-QSL0PE»(X2-X1) 

A=SA**2* ( COS ( THETAW ) +0SL0PE#S IN ( THETAW ) ) **2 
&+SC*#2« ( SIN< THETAW > -OSLOPE«COS ( THETAW ) ) *«2 
D=2 . *B*SC#»2# <QSLOPE»COS ( THETAW ) -SIN( THETAW ) ) 
EcB##2#(SC»«2)-SA**2#(SC*»2)*(C0S< THETAW )+QSLOPE*SIN<  THETAW) )**2 
TEST  FOR  SIGN  OF  DISCRIMINATE 

6 IF  (D**2-4.#A*E>  80.7,7 

REAL  DISTINCT  ROOTS  (DEFINITE  TORSO  AND  BAG  CONTACT) 

7 DISC=(D«»2-4.*A#E>»*.5 

BAG  INTERCEPT  POINTS.  XA.XB  AND  YA.YB 
XA= ( -D-DISC ) / ( 2 . #A ) 

IF(SC.LE. ABS(XA) )XA=ABS(XA)/XA#(SC-.001 ) 

XB* ( -D+DISC ) / ( 2 . »A ) 

IF (SC.LE.ABS(XB) )XB=ABS(XB)/XB*(SC~.001) 

YA= ( B+XA# ( OSLOPE*COS ( THETAW  > -SIN ( THETAW )))/( COS ( THETAW ) +OSLOPE* 
ISIN(THETAW) ) 

IF(SA.LE.ABS(YA))YA=ABS(YA)/YA«(SA-.001) 

YB= ( B+XB* ( OSLOPE#COS ( THETAW ) -SI N < THETAW  >>>/ < COS ( THETAW ) +OSLOPE# 
SSIN( THETAW)) 

IF (SA.LE.ABS( YB) )YB=ABS( YB)/YB#(SA-.001 ) 

ABST=DISTANCE  FROM  POINT  A TO  POINT  B. 

ABST=SQRT ( ( XA-XB  > «#2+ ( YA-YB )*#2 ) 

X AND  Y COORD.  OF  H-POINT  IN  XPRIME.  YPRIME  COORD.  SYSTEM 
XH=< ( Y2-Y1 )#SIN( THETAW ) + ( X2-X1 ) »COS ( THETAW ) ) 

YH=( (Y2-Y1 )-(X2-Xl )#TAN( THETAW) )*<+COS(THETAW> ) 
IF(THETAT-THETAW)8.8, 9 

8 YP=YB 
YN=YA 
XP=XB 
XN=XA 

GO  TO  10 

9 YP=YA 
YN=*YB 
XP“=XA 

xn*xb 

10  RBAG=SQRT<(XN-XH)<»*2+<YN-YH>«*2) 

IF(ABST+RBAG-12.*RN) 11 . 11,12 
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11  XFT=(XA+XB)/2. 

YFT=< YA+YB)/2. 

GO  TO  13 

12  XNECK=XN-(RN*12.-RBAG)#SIN(THETAT-THETAW> 

YNECK=YN+<RN#12.-RBAG ) »COS ( THETAT-THET AW ) 

XFT=  < XN+XNECK ) /2 . 

YFT= ( YN+YNECK ) /2  a 

RFT=DISTANCE  FROM  H-POINT  TO  POINT  OF  FORCE  APPLICATION  ON  TORSO 
13  RFT=SQRT(<XH-XFT>*#2+<YH-YFT>*«2> 

SLOPE  OF  LINE  PERPENDICULAR  TO  AB 

PSLQPE*- ( COS ( TBETAU ) +QSLQPE*SIN  < THET AW ) ) / ( OSLOPE*COS ( THET  AW ) 
l-SIN(THETAW)) 

T=MOM£NT  ARH  OF  TORSO  FORCES 
XQ=-SC 

T = ( -XFT #PSLOPE+YFT- YWH+PSLOPE*XWH ) /SORT ( PSLQPE*#2+i . ) 

Al=l  ./SC**2+PSL0PE#*2/SA**2 

Bl=2 . #PSLQPE/SA#«2# ( YFT-PSLOPE*XFT ) 

Cl=(YFT-PSL0P£#XFT>»«2/SA«*2-l. 

POINTS  G AND  I ARE  THE  POINTS  WHERE  THE  LINE  OF  ACTION  OF  THE  FORC 

WOULD  INTERSECT  THE  BAG. 

XG=(-B1+SQRT(B1##2-4.#A1*C1 ) )/<2.»Al > 
XI=(-B1-SQRT(B1**2-4.«A1«C1>>/(2.*A1) 

YG=PSLOPE« ( XG-XFT ) ♦ YFT 
YI=PSLOPE*(XI-XFT)+YFT 
YRIM^PSLQPE# ( XC-XFT )+YFT 
GI=SQRT(<XI-XG)««2+(YI-YG)»»2) 

CALCULATE  THE  BAG  PENETRATION. 

BP=SQRT ( ( XFT-XG ) #*2+ ( YFT-YG ) #*2 ) 

DETERMINE  THE  MIDPOINT  OF  LINE  GI. 

XHGI-(XG+XI )/2. 

YMGI=(YG+YI )/2. 

CALCULATE  MAJOR  AXIS  LENGTH  OF  ELLIPSE  PERPENDICULAR  TO  TORSO 
ZP=SQRT  < SA#*2» ( 1-XMGI*»2/SC«#2)-YMGI##2> 

CALCULATE  PERIMETER  OF  ELLIPSE  PERPENDICULAR  TO  TORSO 
PERBB=2.»3. I4159»SQRT( ( (GI/2. >#*2+ZP«#2>/2. ) 

SOLVE  FOR  THE  ANGLE (BET A)  THAT  THE  FABRIC  TENSION  FORCE  COMPONENT 
MAKES  WITH  RESPECT  TO  A LINE  NORMAL  TO  THE  TORSO.  FIRST  SOLVE  FOR 
PHI , USING  THE  NEWTON  RAPHSON  METHOD  OF  SOLVING  TRANCENDENTAL  EONS 

LET  PHIO-AM  ESTIMATE  OF  THE  ROOT  PHI  AND  EPSLON=THE  DESIRED  ACC- 
URACY OF  THE  ROOT. 

PHIQ=3 .14 
EPSLON=. 00001 
PHI=PHIO 

IF(ABS(YRIM)-RIMRAD) 17, 16.16 

16  RIH=0. 

G0T018 

"RIM"  IS  THE  CHORD  LENGTH  OF  THE  RIM  AT  POINT  I. 

17  RIM=S0RT(RIMRAD**2-YRIM«*2>*2. 

FOR  BAG  PENETRATIONS  LESS  THAN  ONE-HALF  THE  CHEST  THICKNESS.  THE 
BODY  WIDTH  IN  CONTACT  WITH  THE  BAG  WILL  NOT  EXCEED  THE  LENGTH  OF 
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THE  BODY  IN  CONTACT  WITH  THE  BAG. 

IB  CHESTT=WBACT/2 . 5 

I F ( BP . LT . CHESTT/2 . . AND . WBACT . GT . ABST ) WB= ABST 
DBR=SQRT( (XC-XI >**2+< YRIM-YI >**2> 

SBR*SORT( .25*RXM*«2+DBR#*2>*2. 

21  FPHI=<  PERBB-WB-SBR ) *SIN  ( PHI/2 . ) /PHI- ( GI-BP-DBR ) / ( COS ( ATAN( ( WB-RIM ) 
&/<2.»(GI-BP-DBR>>>)> 

DFPHI=(PERBB-WB-SBR)*CQS(PHI/2. )/(2.#PHI )-( (PERBB-WB-SBR)«SIN(PH1/ 

8>/PHI##2) 

DEL--FPHI/DFPHI 

PHI=PHI+D£L 

IF( PHI .GT.2.#3. 141593. OR.PHI .LT .0. )G0  TO  520 
IF  ( ABS ( DEL ) . LE . EPSLON ) GO  TO  22 
GO  TO  21 

22  ALPHA=ATAN( (WB-RIM)/(2.*(GI-BP-DBR> ) ) 

BETA=PHI/2.+ALPMA 

SOLVE  FOR  THE  ARC  SIN  OF  YP/SA  AND  YN/SA. 

25  ASYPSA=ATAN ( ( YP/SA ) / ( SORT ( 1 . - ( YP/SA ) **2 ) ) ) 

ASYNSA-ATAN < ( YN/SA ) / ( SQRT ( 1 . - ( YN/SA ) ##2 ) ) ) 

BAG,  TLINE  AND  ACHORD  ARE  INTERMEDIATE  VALUES  REQUIRED  FOR  THE 
AREA  OF  INTERCEPT  CALCULATION. 

BAG=SC/<2.#SA)#< (YP*SQRT(SA**2-YP*#2)+SA#*2*ASYPSA)- 
8 ( YN*SQRT ( SA##2-YN##2 > +SA*«2*ASYNSA ) ) 

TLXNE=(XA- < XA-XB )/(YA-YB)*YA)#YP+< XA-XB)/ (2.#( YA-YF' )*YP*«2- 
8 ( XA- ( XA-XB ) / ( Y A-YB ) *YA ) #YN- ( XA-XB  >/(2.*<YA-YB))*  YN*#2 
ACH0RD=0. 

IF(XP.LT.O. )ACH0RD3SC/SA«(3. 141593#SA»*2/2.-( YP#SQRT(SA##2-YP««2) 
&+SA**2#ASYPSA)> 

IF(XN.LT.O. )ACHQRD=ACHQRD+SC/SA*( YN#SQRT (SA#»2-YN**2)+SA»«2#ASYNSA 
S+3.14159*SA»*2/2.) 

SOLVE  FOR  THE  AREA  OF  INTERCEPT. 

AOI=BAG-TLINE+ACHORD 
"VOI"=VOLUME  OF  BAG  INTERCEPT. 

WAVG= ( ( RFT-RBAG ) *WB+RHEAD*WH ) / ( RFT-RBAG+RHEAD ) 

VQI=UAVG*AOI 

THE  FOLLOWING  ROUTINE  USES  THE  NEWTON-RAPHSON  METHOD  TO  SOLVE 
A CUBIC  EQUATION  FOR  THE  BAG  PENETRATION  THATWOULD  EXIST  FOR 
THE  SYMMETRICAL  CASE  WITH  A GIVEN  VOI. 

R0QTG=O.6*BP 

IF(BP.G£.SC+ABS(XC) )ROOTG=SC 
RQOT=ROOTG 

73  FR00T=3 . 14159*SA**2# ( -R00T**3/ ( 3 . #SC**2 ) +RQ0T**2/SC ) -VOI 
DFROOT=3.14159*SA»#2#(-ROQT**2/SC*#2+2.*ROOT/SC> 

DELRT =-FROOT /DFROOT 
ROOT  =ROOT +DELRT 

IF ( ROOT . GT . 2 . #SC . OR . ROOT . LT . 0 . > GO  TO  518 
IF  ( ABS < DELRT ).LE. EPSLON)  GO  TO  74 
GO  TO  73 

74  PER=2.#3.14159#SQRT< (SA**2+SC**2>/2. ) 

BPSYM=RGQT 
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ABSYM=2 .#SA*SQRT ( 1 . - ( SC-BPSYM ) ##2/SC**2 ) 

YC=RXMRAD 

QOSQRT  ( ( XC-XO ) ««2+  ( YC-  Y0 ) »*2 ) 

BCSYHCM  PER-2 . #OOABSYM  > /2 . 

XBS=SC-BPSYM 

YBS=ABSYM/2. 

BEFORE  CALCULATING  THE  LENGTH  OF  THE  LINE  BC  FOR  THE  SYMMETRICAL 
CASE  CHECK  TO  SEE  THAT  THE  HEAD  AND  CHEST  HAVE  NOT  BOTTOMED 
OUT  ON  THE  WHEEL  RIM.  IF  THEY  HAVE.  STOP  THE  RUN. 
IF(XA.LT.XC.AND.YA.LT.O.)  GO  TO  510 
IF ( XA . LT  u XC . AND . YA . GT . 0 . > GO  TO  512 
BCSYMS=SQRT  < ( X8S-XC  ) #»2+ ( YBS-YC ) **2 ) 

THE  FOLLOWING  ROUTINE  USES  THE  NEWTON-RAPHSON  METHOD  TO  SOLVE  A 
TRANCENDENTAL  EQUATION  SO  THAT  THE  RADIUS  AND  ANGLE  OF  THE  VERTICA 

BAG  ENDS  CAN  BE  CALCULATED. 

THE0=3. 14 
THE^THEO 

75  FTHE~BCSYMC*SIN( THE/2 . ) /THE-BCSYMS/2 . 

DFTHE=BCSYMC»COS ( THE/2 . >/<2.*THE)-BCSYMC«SIN<THE/2. >/THE*#2 
BELTH=-FTKE/DFTHE 

THE=TH£+DELTH 

IF(THE.GT.2.#3. 141593.0R.THE.LT.0. )G0  TO  516 
IF  ( ABS ( DELTH ) . LE . EPSLON ) GO  TO  76 
GO  TO  75 

76  RSYM=BCSYMC/THE 

THE  FOLLOWING  STATEMENTS  ARE  USED  TO  CALCULATE  THE  AIRBAG  VOLUME. 
SLOPE=(YBS-YC)/(XBS-XC ) 

Vl=SL0PE»*2»(XC*«3/3.-XBS**3/3.+XBS**2*XC-XC**2*XBS) 

V2=YC**2#(XC-XBS> 

V3=2 . »YC#SLOPE* ( XBS*XC-XBS«*2/2 . -XC**2/2 . ) 

VABCD=ABS(V1+V2+V3>#3. 14159 
ALPHAS=ATAN< SLOPE) 

RC= YC-RC'YM«COS  ( THE/2 . +ALPHAS ) 

RS=RC+4 . /3 .* ( RSYH*S I N ( THE/2 . ) **3 ) *COS ( ALPHAS ) / ( THE-S I N ( THE ) ) 

VBC=2 .*3.1 4 159»RS* ( RSYM**2/2 . * ( THE-SIN ( THE ) ) ) 

VBRIM=3 . 1 4159/6 .#  < SC+XC ) * ( 3 . *YC**2+ ( SC+XC  > *#2  > 

VGL=VABCD+VBC+VBRIM 

CONFINE  BAG  VOLUME  TO  ORIGINAL  VOLUME  IF  THE  BAG  PRESSURE  IS  LESS 

THAN  AMBIENT. 

I F ( VOL . GE . VOLO . AND . PG 1 . LE . 0 . ) VOL=VOLO 
GO  TO  101 

COMPLEX  ROOTS  (NO  TORSO  AND  BAG  CONTACT) 

80  VOL = VOLO 
FT 1=0. 

FTH=0. 

COMPUTE  GAS  FLOW  INTO  BAG 
101  CALL  GASIN(X.QIN) 

SINCE  SUBROUTINE  ■'SOLVE"  CALLS  "BAGSUB"  TWICE  PFR  SOLUTION 
POINT  WE  MUST  DIVIDE  THE  TIME  STEP  BY  2. 

DELTAT=STEP 
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IF ( IL00P.NE. 1 )DELTAT=0 
GU1=GW+QIN»DELTAT 

TEST  TO  SEE  IF  BAG  DEPLOYED  YET: IS  BAG  PRESS. >ATMOS.  PRESS . 7 

IF(X.EQ.0.)DTIME=.25 

IF ( PA . GE . ATHOP . OR . X . GT . DT I ME ) GO  TO  107 

GW=GW1 

PA=GW*GTZ*U/VQL 

PG1*PA-ATM0P 

GT=GTZ 

FFT=0. 

FP=0. 

FTT=0. 

FHEAD=0. 

FTH=0. 

DTIM£=X+STEP 
GO  TO  490 

COMPUTE  NEW  TEMP.  AND  PRESS.  DUE  TO  NET  GAS  GAIN  IN  BAG 

107  GT7=(GW*GT+QIN#GTZ*DELTAT>/GU1 
PNUH=U»GT7#GW1 
PA7=PNUM/VOLO 

COMPUTE  NEW  GAS  PRESS.  AND  TEMP.  DUE  TO  POLYTROPIC  COMP.  OR  EXPANS 

PA8= (PNUM/VOL ) ##FPN/PA7»# ( FPN-1 . ) 

GT8=GT7# ( PA8/PA7 ) »*  < ( FPN- 1 . ) /FPN ) 

BAG  VENTING  COMPUTATIONS;  FIRST  CALC.  PRESS.  RATIO  ACROSS  VENT 
PR7=ATMQP/PA8 

TEST  FOR  CHOKED  FLOW;  ALSO.  IF  PR7>1. BYPASS  QEXH.&  SET  GU=GW1 
IF  (PR7.LT. PR8)  GO  TO  108 
IF  (PR7.GE.1.)  GO  TO  110 

FH1=VC1*SQRT ( PR7»* ( 2 . /PN1 ) -PR7** ( (PN1+ 1 . ) /PN1 ) ) 

GO  TO  109 

108  FM1=FM2 

COMPUTE  EXHAUST  FLOW  AND  RESIDUAL  GAS  WEIGHT 

109  QEXH=SQRT ( ( 772 . *PN1 ) / ( PN1- 1 . ) ) * AV*PA8«FM 1/SQRT ( U*GT8 ) 
GW=GW1-QEXH*DELT  AT 

GO  TO  111 

110  GW=GW1 

COMPUTE  PRESS.  AND  TEMP.  OF  GAS  AFTER  VENTING 

111  RATIQ=GW/GW1 
PA=PA8*RATI0#»PN1 
GT=GT8»RATI0*«(PN1-1) 

COMPUTE  PRESS.  RATIO  TO  DETERMINE  WHETHER  GAS  COMP.  OR  EXPANDED  TH 

TIME  THRU  LOOP;  THEN  SET  PROPER  POLYTROPIC  EXPONENT. 

PR6=PA8/PA5 

IF  (PR6J.T.  1.0001)  GO  TO  112 

FPN-PN2 

GO  TO  113 

112  FPN=PN3 

COMPUTE  BAG  PRESSURE. 

113  PG1=PA-ATM0P 

IF  THE  BAG  PRESSURE  IS  NEGATIVE.  CALL  IT  0.  FOR  BAG  FORCE  CALCS. 


DRACR 

13980 

14000C 

14020C 

14040 

14060C 

14080 

14100 

14120 

14140 

14160 

14180 

14200 

14220 

14240 

14260 

14280 

14300 

14320 

14340 

14360 

14380 

14400 

14420 

14440 

14460 

14480 

14500 

14520 

14540 

14560 

14580 

14600 

14620 

14640C 

14660C 

14680C 

14700 

14720 

14740 

14760 

14780 

14800- 

14820 

14840 

14850 

14855 

14860 

14900 

15000 

15020 


16  - 


07/28/82 


IF(PG1 .LT.O. )PG1=0. 

IF  THE  TORSO  IS  NOT  IN  CONTACT  WITH  THE  BAG.  SKIP  THE  BAG  FORCE 
CALCULATION. 

IF(D##2-4.#A#E)490, 115, 115 

THE  BAG  FORCES  ARE  CALCULATED  IN  THE  NEXT  SEVERAL  STATEMENTS. 
115  ENPHI=PGl*SQRT(SA**2-WB**2/4.>»(SQRT<SC*«2*<l-WB**2/<4.# 

ISA#*2) ) ) )»( 1-BP/GI )/ (2.#SQRT (SA**2-WB*«2/4.+SC##2#( 1-WB**2 
&/(4.*SA»#2>>>> 

FFT=-2 . »ENPHI»ABST*COS ( BETA ) 

IF  (BETA. LE. 1.5708)  FFT=0. 

FP=PG1«WB*ABST 

I F ( RBAG*ABST-RN*12. ) 140 . 1 40 . 1 35 

135  1F(ABST+RBAG.GT.RT0PH«12. ) GO  TO  136 
RHEAD= ( ABST +RBAG-RN* 12 . )/2. 

GO  TO  137 

136  RHE AD= ( ( RT OPH-RN ) # 1 2 . ) /2 . 

137  HEADT=WHACT 

I F ( BP . LT . HEADT/2 . . AND . WHACT . GT . 2 . »RHEAD ) WH=2 . #RHEAD 
FHEAD=2 . *RHEAD«WH*PGl+2 . #RHEAD*FFT/ABST 
IF(Y(6)-Y(7) .LT.O. )FHEAD=FHEAD#C0S(Y(6)-Y(7) ) 
FTH=TNECK-FHEAD#RHEAD/12. 

FP=FP* (RN«12. -RBAG ) /ABST 
FFT=FF  T* ( RN* 1 2 . -RBAG ) /ABST 
GO  TO  141 

140  FHEAD=0 . 

FTH=TNECK 

141  FCHEST=FFT+FP 

I F ( FHEAD«RHEAD-ABS < TNECK ) * 12 . >142.143.143 

142  TRANST0R=-FHEAD*RHEAD/12. 

GO  TO  144 

143  IF(TNECK.LE.O. )TRANSTOR=TNECK 
IF ( TNECK . GT . 0 . ) TRANSTOR=-TNECK 

144  FTT =-FHEAD#RN*COS ( Y ( 6 ) - Y( 7 ) ) +TRANST0R-RF T*FCHEST / 12 . 

THIS  PART  OF  THE  PROGRAM  COMPUTES  THE  DYNAMICS  OF  STEERING 
WHEEL  CRUSH. 

FAW=FCHEST*COS(THETAW-THETAT)+FHEAD*COS(THETAW-THETAH)+GS#WWH/ 
$32. 17#C0S( THETAW) 

FNW=FCHEST*SIN(THETAW-THETAT)+FHEAD*SIN(THETAW-THETAH)+GS 
&#WWH/32. 17#SIN ( THETAW ) 

CALL  WHF(SWHC.FRU) 

IFTFRW.EQ-O. )FRW=1 . 

DELFW=FAW-FRW 
IF(DELFW) 190, 192, 192 
190  EAWC=BIG 
GO  TO  195 

192  VWHC=VWHCQ+OELFW«DELTAT/MWH 

SMHC=SWHC0+VWHCQ*DELTAT+DELFW#DELTAT««2/ ( 2 . *MWH ) 

EAWC=FRW/C0S( THETAW) 

X1WC=X1- (SWHC-SWHCO ) *COS( THETAW ) 
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Y1 WC=Y 1 - ( SWHC-SWHCO ) #SIN ( THET  AW ) 

THIS  PART  OF  THE  PROGRAM  COMPUTES  THE  DYNAMICS  OF  STEERING 
COLUMN  CRUSH. 

COMPUTE  THE  BENDING  MOMENT  APPLIED  TO  THE  COLUMN  AT  THE  WHEEL  HUB. 
195  HO--*FCH£ST»T 

COMPUTE  THE  NORMAL  FORCE  APPLIED  TO  THE  COLUMN. 

FNC=FCHEST*SIN<  THETAC-THETAT  >+FHEAD*SIN<  THETAC-THETAH ) 
FNC0L=FNC+GS/32 . 17# ( WC+WWH ) #SIN ( THET AC ) 

COMPUTE  THE  REACTIONS  AT  THE  COLUMN  SUPPORT  POINTS. 

RSC=FNCOL+ ( MQ+FNCOL#LSC ) /LFW 

RBA= ( RSC* (LBF-LSC ) +FNCOL" ( LSC-LBF ) ) / ( LBF-LBA ) 

RBF= { RSC» ( LBA-LSC ) +FNCOL* ( LSC-LBA ) ) / ( LBF-LBA ) 

FACOL=FCHEST#COS ( THETAC-THETAT ) +FHEAD#COS ( THETAC-THETAH )+GS«< UC+ 
&WWH)/32.17»C0S(THETAC) 

CALL  COLF ( SCOLC , FCOL ) 

IF(FCOL.EQ.O.)FCOL=1. 

FRCOL=FCOL+MU#ABS ( RBF+RBA+RSC ) 

DELFC=F  ACOL-FRCOL 
IF ( DELFC ) 205 . 207  f 207 
205  EACC=BIG 
GO  TO  208 

207  VCQLC=MCOLCO+DELFC*DELT AT /MC 
SCOLC-SCOLCQ+VCOLCO*DELTAT+  DELFC«DELTAT**2/ ( 2 . «HC ) 

EACC=FRCOL/COS ( THET AC ) 

X1CC=X1- ( SCOLC-SCOLCO ) «COS  < THET AC ) 
Y1CC=Y1-(SC0LC-SC0LC0)*SIN(THETAC) 

208  IF ( EAWC-EACC ) 210 . 230 . 220 
210  EAC=EAWC 

X1C=X1WC 
Y1C=Y1WC 
GO  TO  240 
220  EAC=EACC 
X1C=X1CC 
Y1C=Y1CC 
GO  TO  240 
230  EAC*BIG 

THIS  PART  OF  THE  PROGRAM  COMPUTES  THE  DYNAMICS  OF  STEERING 
WHEEL  ROTATION. 

240  WIN=MWH*RIMRAD#*2/2.+MIP#RIP**2/4. 

TAW=FCNEST*T 
RELWR=THET  AW-THET  AC 
IF(RELWR.GE.O. )G0  TO  244 
RELWR=-RELWR 
CALL  WHROTN ( RELWR , TRW ) 

GO  TO  248 

244  CALL  WHROTPI RELWR, TRW) 

248  IF(TRW.E0.0.)TRW*1. 

TRW=- ( TAW+1 . ) /ABS ( (TAW+1. ) ) #TRW 
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IF(TAWoE0»0. ) TRW=+ ( 0MW0+ 1 . ) / APS ( OMWO+ 1 . )*TRW 
DELTW=TAW+TRW 

IF(ABS(TAW)-ABS(TRW>  >250,250,260 
250  EAWR=BIG 
GO  TO  261 

260  ALW=D£LTW/WIN 
OMW=OHWQ+ALW#DELT AT 

THETAW=THETAW0+0MW0«DELTAT+ALW«DELTAT**2/2. 

DTHW=THETAW-THETAWO 

X1WR*X1-XWH«(COS(THETAW)-COS(THETAWO)> 

Y1WR=Y1-XWH#(SIN(THETAW)-SIN(THETAW0>) 

IF (XI .EQ.X1WR)X1WR=X1+.0001 
EAWR=ABS(TRW*DTHW/(X1-X1WR> ) 

THIS  PART  OF  THE  PROGRAM  COMPUTES  THE  DYNAMICS  OF  STEERING 

COLUMN  ROTATION,, 

261  C0I=MC#RC0L»#2+ ( MWH+M IP ) * ( LFW+LSC ) #*2 
IF(SC0LC.LT.2. )G0  TO  270 
RELCR=THETAC-THETACZ 
IF(R£LCR.GE.O. )G0  TO  263 
RELCR=-RELCR 

CALL  CQLRQTN(RELCR, REACT) 

GO  TO  265 

263  CALL  COLROTP(RELCR, REACT) 

265  IF ( REACT „LQo0„ )REACT=1 . 

IF(NC0LTYPE-1 >268,267,268 

267  RBRKT=REACT 
TAC=FNCOL* (LFW+LSC ) +M0 
TRC=-(TAC+1. )/ABS( (TAC+1. ) >«RBRKT*LFW 
IF ( TAC . EO .0 . ) TRC=- ( OMCO+ 1 . ) /APS ( OMCO+ 1 . >«RBRKT*LFW 
DELTC-TAC+TRC 

IF ( APS ( TAC ) - ABS ( TRC ) ) 270 . 270 . 280 

268  TAC=MQ+FNCGL*LSC 
TRC=- ( TAC+1. )/ABS< (TAC+1. )>*REACT 
IF(TAC.EQ.O. )TRC=-(QMC0+1 . )/ABS(OMCO+l . )+REACT 
DELTC=TAC+TRC 

IF ( ABS ( TAC ) - ABS ( TRC ) ) 270 . 270 . 280 
270  EACR=BXG 
GO  TO  300 
280  ALC=DELTC/COI 

OMC=ONCO+ALC«DELTAT 

TH£TAC=TH£T  ACO+OMCO*DELTAT +ALC+DELTAT  *.*2/2 . 
DTHC=THETAC-THETACO 
RPIUOT=LFW+LSC+XWH«CQS(THETAW-THETAC> 
X1CR=X1+RPIV0T#(C0S(THETAC)-C0S(THETAC0>) 
Y1CR=Y1+RPIV0T«(SIN(THETAC)-SIN(THETAC0) ) 

IF (XI .EQ.X1CR)X1CR=X1+.0001 
EACR=ABS(TRC*DTHC/(X1-X1CR) ) 

300  IF ( EAWR-EACR ) 310 , 330 , 320 
310  EAR=EAWR 


(Note:  Statement  16310  added  for  the 
specialized  test  case  - not  in  the 
original  DRACR  program) 
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X1R=X1WR 
Y1R-Y1WR 
GO  TO  400 
320  EAR=EACR 
X1R-X1CR 
Y1R=Y1CR 
GO  TO  400 
330  EARWIG 

400  IF(EAR-EAC>410,450.430 
410  EA=EAR 
X1=X1R 
Y1=Y1R 

OWHC=OWHCO-FRW«DELT  AT/MUH 
IF ( VWHC.LT.O. )VWHC=0. 

SmC==SWHCO+OWHCQ*DELTAT-FRW«DEL  TAT**2/  ( 2 . *MWH ) 

IF  < SWHC . LT . SWHCO ) SWHC=SWHCO 
X 1 =X1- ( SWHC-SWHCO ) »COS ( THETAW ) 
Yl=Yl-(5yHC-SWHC0)*SIN<THETAU) 

VWHCO=VWHC 

5W»:C0=SWHC 

OCOLC=MCOLCO-FRCOL*DELTAT/MC 
IF (VCOLC.LT.O. >VCQLC=0. 

SCGLC=SC0LC0+VC0LC0*DELTAT-FRC0L«I>ELTAT**2/ ( 2 . *MC ) 
IF ( SCOLC . LT . SCOLCO )SCQLC=SCQLCO 
X1=X1- < SCOLC-SCOLCO ) »COS ( THETAC > 

Yl=Yl-( SCOLC-SCOLCO ) «SIN ( THETAC ) 

VCQLCO=VCOLC 

SCOLCO=SCOLC 

IF(EAWR-EACR>417,421,421 

417  IF(ABS(OMCO) >418,418.419 

418  QMC=OHCO 
THETAC«THETACO 
GO  TO  420 

419  OMC*GMCG+TRC/COI«DELTAT 
IF(GMC#OMCO.LT.O.)0MC=O. 

THET  AC=THET  ACQ-MDMCO»D£LT  AT +TRC/  ( 2 . «C0 1 ) *DELTAT#«2 
THETAW=THETAW+THETAC-THETACO 
X1-X1+RPIV0T# ( COS( THETAC  > -COS ( THETACO > > 
Y1=Y1+RPIV0T#( SIN < THETAC )-S IN (THETACO) ) 

OHCO=OMC 

THETACO=THETAC 

420  0NW0=0MW 
THETAWO=THETAW 
GO  TO  490 

421  IF (ABS(OHWO) >422.422.423 

422  QHW=OHWO 

THET  AW=THET  AWO+DTHC 
THETAWO=THETAW 
GO  TO  424 

423  OMW=OHWO«-TRW/WIN#DELTAT 
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17332 

IF  ( 0MW»0HW0  .LT.O. ) 0MW«0 . 

17335 

THET  AW*THET  A W^MWQ#DEIT  AT+TRW/  <2.#UIN)»DELTAT**2+DTHC 

17345 

XI *Xi-XWH» ( COS ( THE! AW ) -COS ( THET AWO ) ) 

17350 

YlsYl-XWM«(SXN(THETAW>-SIN<THETAWQ>  > 

17355 

ohwo^omw 

17360 

THETAWO«THETAW 

17365 

424 

OMCO»OMC 

17370 

THETACOsTHETAC 

17375 

GO  TO  490 

17380 

430 

EA«EAC 

17385 

Xl^XiC 

17390 

Y1SY1C 

17395 

I F ( ABS  ( OMWO ) >43 1 , 43 1 , 432 

17400 

431 

om«mm 

1740S 

TH£TAy=THETAWO 

17410 

GO  TO  433 

17415 

432 

OHW=OMWO+TRW/y IN»DELT AT 

17417 

XF(QHW#OHWO.LT.O. >OMW=0. 

17420 

THET  AW^THCT  AWO+OHWO»DELT  AT +TRW/ ( 2 . «WIN ) «DELTAT**2 

17425 

Xl^X! -XWH* ( COS ( THET AW )-COS ( THET AWO ) ) 

17430 

YlffiYl-XyH*(SIH(THETAW)-SIN(THETAWO) > 

17435 

Of1tfO«ONW 

17440 

THET  AWO*THET  AW 

17445 

433 

IF(ABS(OHCO) >434,434.435 

17450 

434 

OfCsOHCO 

17455 

THETAC®THETACO 

17460 

GO  TO  436 

17465 

435 

OMCsOHCO^TRC/COI«OELTAT 

17467 

IF(OMC*OMCQ.LT.O. >0MC=0. 

17470 

THETAC=THETACO+OHCO»DELTAT+TRC/( 2.«C0I >#DELTAT**2 

17472 

THETAy^THETAW^THETAC-THETACO 

17475 

Xl®Xl+RPIUOT# ( COS ( THETAC  > -COS ( THETACO ) > 

17480 

Yt®Yl^PXVOT«(SIN(THETAC>-SIN<THETAC0>> 

17485 

OMCO«OMC 

17490 

THET ACO® THET AC 

17492 

THET AWO® THET AW 

17495 

436 

IF  C EAWC-EACC >439,443,443 

17520 

439 

OCOLC®OCQLCO“FRCOL»OELTAT/MC 

17522 

IF  (VCOLC.LT.  Co  >VC0LC*0. 

17525 

SCQLJ>SC0LC0-H/C01C0#I>ELT  AT-FRCOL#I>£LT  AT*  «2/(2.  *MC ) 

17527 

IF ( SCQLC . LT .SCOLCO  >SCOLC=SCQLCQ 

17530- 

„ . . _ 

SCOLC-SCOLCO  > #COS ( THETAC  > 

17535 

r 1«Y 1 ■ - ( SCOLC-SCOLCO  > « S I N ( THET  AC ) 

17540 

VCOLCQ*VCOLC 

17545 

SCOLCO*SCOLC 

17550 

VyHCO»WHC 

17555 

swco*syHc 

17560 

GO  TO  490 

17585 

443 

VWHC«WHCO“FRW»I>ELTAT/HWH 

17587 

IF<VWHC.LT.0.)WHC«O. 
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SWHC“SWHC0+VWHC0*DELT  AT-FRW»DELTAT#*2/ ( 2 . *MWH ) 

I F ( SWHC . LT . 8WHC0 ) SUHOSUHCQ 
X 1«X1- ( SWHC-SWHCO ) *C0S ( THETAW ) 
Y1«Y1-<SWHC-8WHC0)*SIN(THETAW> 

VWHC0*VWHC 
SWHCO-SWHC 
VC0LC0»VC0LC 
SCOLCO-SCOLC 
CO  TO  490 
450  EA*BIG 

VWHC*VWHCO-FRW#DELTAT/HWH 
IF  (VWHC.LT.O.  )VWHC=*0. 

SWHC*SWHCO^HCO*DELTAT-FRW*DELTAT*«2/  ( 2 . *MWH ) 

IF ( SWHC . LT . SUHCO ) SWHC-SWHCQ 
X1»X1- ( SWHC-SWHCQ ) #COS ( THETAW ) 

Y 1*Y l-( SWHC-SWHCO >»S IN (THETAW) 

VWHCO*VWHC 

SWHCO«SWHC 

VCOLCaVCOLCO-FRCOL»DELT  AT /MC 
IFCVCOLC.LT.O. >VC0LC*0. 

SCOLC*SCOLCO+VCQLCO*D£LT  AT-FRCOL*DELT  AT  *#2/ ( 2 . *MC ) 
IF < SCOLC . LT . SCOLCO ) SCOLC»SCOLCO 
XI *X1 - ( SCOLC-SCOLCO ) *COS ( THET AC ) 

Y 1 =Y1- < SCOLC-SCOLCO ) #S IN ( THETAC ) 

VCOLCO=OCOLC 

SCOLCO=SCOLC 

IF( ABS ( OMWO ) ) 457 , 457 , 458 

457  QHW*OHWO 
THETAW*THETAWO 
GO  TO  459 

458  OMW*GHWO+TRW/WIN«DELTAT 
IF(OHW*ONWO.LT.O. )0MW=0. 

THET  AW«THETAWO+OMWO#DELTAT+TRW/  ( 2 .*WIN ) *DELTAT#*2 
X 1*X1-XWH« ( COS ( THETAW  > -COS ( THET A WO ) ) 
Y1*Y1-XWH«(SIN(THETAW)-SIN(THETAW0) ) 

OMWO-OMW 

THETAWO-THETAW 

459  IF<ABS(OHCO) >460,460,461 

460  OMC*OHCO 
THETAC-THETACO 
GO  TO  490 

461  OMC=GMCO+TRC/CQI*DELTAT 
IF(OHC*OHCO.LT.O. )0HC»0. 

THETAC=THETACO+OMCO*DELTAT+TRC/ ( 2 . *C0 I ) #DELTAT**2 
THET  AW*THETAW+TH£T  AC-THET  ACO 
X 1 *X1 ♦RPI VOT*  < COS ( THETAC ) -COS ( THET ACO ) ) 
Y1=Y1+RPIVOT«(SIN(THETAC)-SIN(THETACO) ) 

OHCO*OHC 

THETACO-THETAC 

THETAWO=THETAW 
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UPDATE  OLD  VALUES. 

490  VOLO*»VOL 
PA5*PA 

LBAmLBAZ”SCOLC 

IFCLBA.LE.O.)  GO  TO  514 

LRWHZ-SCOLC 

WB^JBACT 

WH»HHACT 

R£L8F*Y(5)-Y(8) 

CALL  SPRING ( SFN . SFUT . RELSF . SUN . 0 . 0 , SF . NSF) 

HPRB*(Y<5>-Y<8>>»12. 

IF  C HPRD-DELTTB/COS ( THLO )) 492, 492, 493 

492  THEF^THFO 
TH£LSTKLQ 
RELKR®Y(5)-Y<8) 

GO  TO  499 

493  IF(D£LTTB.LT.O. >GQ  TO  522 
NL00P=1 

R0TG»(Y2+LT)/2. 

ROT«R0TG 

495  FK8T*R0T» Y2+X2«SQRT ( LT**2-R0T*«2 ) - ( LT»*2-LF««2+X2**2+Y2**2 )/2. 
DFR0T=Y2-X2«R0T/SQRT ( LT*«2-R0T*#2  > 

DELR8T  E~FROT /BFROT 
ROT~ROT+D£LROT 

ZFIROT.GT.LDROT'CT-.Ol 
I F < A®S  < PELROT ) » LE . EPSLON ) GO  TO  496 

NLOOP*NLOOP+1 

IF  (SCOOP . GT.  20.  >G0  TO  524 

GO  TO  495 

496  THEF*ATAN ( ( ( ROT- Y2 >/LF)/( SQRT (l.-( (RQT-Y2)/LF)**2) ) ) 

THELSATAN ( ( ROT/LT ) / ( SQRT  <1.-<R0T/LT)*#2) ) ) 

R£LKR*LT/i2 . « ( COS ( THLO ) -COS ( THEL ) ) +DELTT  B/< 12. *COS ( THLO ) ) 

499  CALL  SPRING (KRN, RUT, RELKR.SKR.O.O.FKNEE.NKR) 

FX®- ( 8F +FKNEE+FCHEST#COS  ( Y ( 7 ) >+FHEAD«COS< Y<6>  > ) 

N®1 

GO  TO  540 

500  XSTOP®X 

PRINT  5 

IF (N.EQ.2)PRINT  ."DRIVER  SUBMARINING.  NOT  RECOVERABLE,  RUN  STOPPED 

GO  TO  540 
510  XSTOPsX 
N«#H 
PRINT  5 

IF(N.EQ.2)PRINT  ."CHEST  IMPACT  WITH  LOWER  WHEEL  RIM.  RUN  STOPPED." 

GO  TO  540 
512  XSTOP^X 
N=N*1 
PRINT  5 
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IF(N.EQ.2)PRINT  ."BODY  IMPACT  WITH  UPPER  WHEEL  RIM.  RUN  STOPPED." 
CO  TO  540 
514  XSTOP=X 
N*N+1 
PRINT  5 

IF(N.EQ.2)PRINT. "COL. STROKE  > SHEAR  CAP.  DIM. (LBAZ) . RUN  STOPPED" 
GO  TO  540 
516  XSTOP«X 
M=N+1 
PRINT  5 

IF(N.EQ.2)PRINT, "THETA  NOT  CONVERGING,  RUN  STOPPED." 

GO  TO  540 
518  XSTOP=X 
N«=N+1 
PRINT  5 

IF(N.E0.2)PRINT,"BAG  PENETRATION  SOLUTION  NOT  CONVERGING,  RUN 
ISTOPPED." 

GO  TO  540 
520  XSTOP=X 
N-N+l 
PRINT  5 

IF(N.EQ.2)PRINT,"PHI  NOT  CONVERGING,  RUN  STOPPED." 

GO  TO  540 
522  XSTQP*X 
N=N+1 
PRINT  5 

IF(N.EQ.2)PRINT,"X2  INPUT  NOT  CONSISTENT  WITH  OTHER  LEG  INPUT." 

GO  TO  540 
524  XSTOP=X 
N®N+1 
PRINT  5 

IF(N.EQ.2)PRINT."KNEE  ANGLE  SOLUTION  NOT  CONVERGING." 

540  RETURN 
END 

THIS  SUBROUTINE  COMPUTES  THE  RATE  THAT  GAS  ENTERS  THE  BAG. 
SUBROUTINE  GASIN(X.QIN) 

CQMMON/GASFLO/NPG , GEN ( 2 . 1 2 ) 

CALL  LOOKUP (X, GEN, NPG.QIN) 

RETURN 

END 

THIS  SUBROUTINE  COMPUTES  THE  WHEEL  FORCE  AS  A FCN  OF  STROKE. 
SUBROUTINE  WHF ( SWHC , FRW ) 

COMMON/UHFOR/NPW , WHE ( 2 , 1 2 ) , VWHCO , SWHCO 
CALL  LOOKUP (SWHC, WHE, NPW, FRW) 

RETURN 

END 


THIS  SUBROUTINE  COMPUTES  THE  COLUMN  FORCE  AS  A FCN  OF  STROKE. 
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SUBROUTINE  COLRSCQLC,FCQL) 

CCfSHON/CGLFOR/NPC , COL  ( 2 , 1 2 ) 

CALL  LOOKUP ( SCOLC , COL , MFC , FCOL ) 

RETURN 

Bffi 

THIS  SUBROUT ISC  COMPUTES  THE  WHEEL  MOMENT  AS  A FCN  OF  POSITIVE 

(UPWARD)  ROTATION. 

SUBROUTINE  UHROTP<RELyR»TRW) 

COMMON/WROTP/NPWRP  , WH8P  (2,12).  OMWO , THETAWO 

CALL  LOOKUP <RELHR,WHRP,NPMRP.TRM) 

RETURN 

EMB 

TIf®WiiiM§«TIWPUT£S  THE  WHEEL  MOMENT  AS  A FCN  OF  NEGATIVE 

SUBROUTINE  WHROTN ( RELUR , TRW  > 

COW0M/WTOTN/NWJRN , WHRN  (2.12) 

CALL  LOOKUP ( RELUR , WHRN , NPWRN . TRW ) 

RETURN 

END 

THIS  SUBROUTINE  COMPUTES  THE  COLUMN  MOMENT  (OR  FORCE)  AS  A FCN 

OF  POSITIVE  (UPWARD)  ROTATION. 

SUBROUTINE  COLROTP(RELCR, REACT) 

C0MMOM/CR0TP/NPCRP , COLRP (2,12) ,OMCO, THETACQ 
CALL  LOOKUP (RELCR. COLRP, NPCRP, REACT) 

RETURN 

EMB 

THIS  SUBROUTINE  COMPUTES  THE  COLUMN  MOMENT  (OR  FORCE)  AS  A FCN 

OF  NEGATIVE  (DOWNWARD)  ROTATION. 

SUBROUTINE  COLROTN( RELCR, REACT) 

COMMON/CROTN/NPCRN , COLRN (2.12) 

CALL  LOOKUP (RELCR , COLRN , NPCRN , REACT ) 

RETURN 

END 

THIS  SUBROUTINE  PLACES  CERTAIN  VALUES  IN  MATRIX  FORMAT  FOR  PRINTIN 

SUBROUTINE  PRINTICX.Y.BY) 

COMMON  /0UT1/N9.X9C4. 175) .19(175) 

DOUBLE  PRECISION  Y(8) 

DIMENSION  BY  (8) 

CQHHON/MANBAT/ZL , ZT , ZH , RT , RN , RH . RT  OPH , X2Z , Y2Z , WB . WH 

C0WOM/HIC/THIC  (175),  HRGS  ( 175 ) 

COPWON/TIME/STEP „ XST OP , P INT 1 , ILOOP 

n?*n?+i 

IF CN9.GT. 175) RETURN 
CH«C0S(Y(6)> 

CT»C0S(Y<7> ) 
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SH«8IN(Y<6)> 

ST=SIN(Y(7) 1 

X9(ltN9>«<DY<l>*CT+RT*DY(3))/32.17 
X9(2,N9)“-<DY(1)*ST-RT*Y(3)*Y(3) 1/32. 17 
X9 ( 3 f N9 1 B ( RH#DY ( 2 1 +DY ( 1 ) #CH+RN#  < DY ( 3 1 » < CH«CT +SH#ST 1 + Y < 3 1 
&« Y ( 3 ) « ( SH*CT-CH»ST 1)1/32.17 

X9<4,N?)»<RH*Y<2)«Y(2)+RN#<Y(3)»Y<3)«<SH«ST+CH*CT)-DY<3> 

S* ( SH#CT-CH#ST) 1 -DY  < 1 1 «SH 1/32.17 
T9(H9)*X«1000. 

THIC<N9)=<N9»PIMT1 )**0.6 
HRGS^=SQRT(X9(3,N9)**2+X9<4,N9>«#2> 

RETURN 

END 

THIS  SUBROUTINE  SOLVES  THE  DIFFERENTIAL  EQUATIONS  THAT  DETERMINE 
THE  DRIVER  KINEMATICS  AND  VEHICLE  MOTION.  THE  FOURTH  ORDER  RUNGE- 
KUTTA  METHOD  IS  USED  TO  START  THE  INTEGRATION.  BUT  ONCE  THE  FIRST 
FOUR  POINTS  ARE  OBTAINED  WE  SWITCH  TO  THE  MORE  ECONOMICAL  FOURTH 
ORDER  ADAMS-KOULTON  PREDICTOR-CORRECTOR  METHOD. 

SUBROUTINE  SOLVE(N) 

COMMON/TIME/STEP, XST0P.PINT1 . ILOOP 
DOUBLE  PRECISION  Y<8),YT<8) 

REAL  DY 

DOUBLE  PRECISION  B270.B19.B251 
DIMENSION  DY(8)fF<3,8)fJ<3),8<8) 

B251=251 . 

B270=270. 

B19=19. 

PRI1=0. 

PR I 2=0. 

XSTGP=. 120 
STEP*. 001 
CALL  SETUP <X,Y) 

CALL  DIFEQ<X,Y.DY) 

PINT 1». 005 
PINT2®.005 
CALL  PRINT1(X,Y,DY) 

CALL  PRINT2(X,YfDY) 

I F C XSTOP . CT . . 25 1 XSTOP= . 25 

START  OF  INTEGRATION  ROUTINE 

RUMGE  KUTTA  START  UP 

J(l)=l 

J<2)=2 

J(3)a3 

DO  7270  K=1.N 
7270  F(3.K)*DY(K) 

DO  7550  JK=1 ,3 
DO  7300  K=1,N 
7300  S(K1=DY(K)«STEP 
XN«X+STEP/2. 
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DO  7330  K=1,N 

2oee© 

7330 

YT(IO»Y<KHS(K)/2. 

20900 

CALL  MFEQ(XN„YT»DY> 

20920 

M 7360  K®1,N 

20940 

736© 

S(K>=S<KH2a«DY<K)*STEP 

20960 

DO  7380  KsleM 

20980 

7380 

YT(K)®Y(K)+STEP#DY(K)/2. 

21W® 

CALL  MFEQCXN.YT.DY) 

21020 

DO  741©  K«1  ,N 

21040 

7410 

SCK)=S(K>-«'2ctt0Y(K)»STEP 

21060 

DO  7430  K~ioH 

2ioeo 

7430 

YT(K)”YCK)+DY(K)»STEP 

21100 

X^X-^STEP 

2112© 

CALL  DIFEO<X,YT,DY) 

21140 

DO  7470  K*1,N 

21160 

7470 

Y {K ) =Y ( K > +( S ( K ) +DY ( K ) «STEP ) /6 . 

21180 

CALL  MFEQ<X,Y„DY) 

21200 

GOTO  < 7500 , 7530 , 7550  > , JK 

21220 

7500 

DO  7510  K»1,N 

21240 

7510 

F < 2 , 10 ®DY ( K ) 

21260 

GO  TO  7550 

21280 

7530 

00  7540  K=1»N 

21300 

7540 

F(lffK)KDY(K) 

21320 

7550 

FRli=*X 

21340 

FRI2-X 

21360C 

PREDICTOR-CORRECTOR  SECTION 

21360C 

PREDICTOR 

21400 

7590 

&9  760©  K=1,N 

21420 

7600 

YT<K)sY(K)^STEP*(55.*DYCK)-59.»F( J< 1 ) ,K)+37.#F( J(2) ,K) 

2144© 

»-9.*F(J(3).K))/24. 

21460C 

SAVE  DY?S 

21480 

DO  7640  K*1,N 

215#© 

7640 

FCJ(3> „K)SDY(K> 

21520C 

EVALUATE  STEP 

21540 

X*X*STEP 

21S50 

I LOOP*! 

21560 

CALL  DIFEp(X,YT,DY> 

215SOC 

ROTATE  VECTOR  POINTER 

21600 

UTSJC3) 

21620 

J(3>»J<2> 

21640-  - 

■ - . 

21640 

J(1)«UT 

21480C 

CORRECTOR 

21700 

DO  7750  K»lfN 

21720 

7750 

Y(K)»Y(K)+STEP»(9.*DY(K)+19.«F(J(1),K)-5.«F(J(2).K)+F(J(3) 

21740 

UC)>/24. 

21780C 

ADDITION  OF  ERROR  TERM 

21800 

10  780©  K-l.N 

21820 

780© 

Y(K>=<B251*Y<K>+B19#YT(K>)/B270 

21840C 

SEC QM  EVALUATION  STEP 

o 

-7 
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21050 

21060 

21870 

21880C 

21900 

21920 

21940 

21960 

21900 

22000 

22020 

22040 

22060 

22120 

22140 

221 60C 

2210OC 

22200C 

22220 

22240 

22250 

22260 

22200 

22300 

22320 

22340 

22360 

22300 

22400 

22420 

22440 

22460 

22400 

22490 

22500 

22520 

22540 

22560 

22580 

22600 

22610 

22620 

22640 

22660 

22680 

22700 

22720 

22740 

22760 

22780 
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IL00P»2 

CALL  DIFEQ(X,Y.DY> 

CALL  UPDATE <X,Y,DY) 

PRINTING  SECTION 

PRI1=PRI1+STEP 

PRI2»PRI2+STEP 

IF(PRI1 .LT.PINT1 )G0T07890 

PRI1=PRI1-PINT1 

CALL  PRINTHX.Y.DY) 

7890  IF(PRI2.LT.PINT2)G0T0792O 
PRI2*PRI2-PINT2 
CALL  PRINT2(X,Y,DY) 

7920  IF(X.LT.XST0P)G0T07590 
RETURN 
END 

THIS  SUBROUTINE  COMPUTES  THE  KNEE  RESTRAINT  CRUSH  FORCE  AND  THE 
SEAT  FRICTION  FORCE.  HYSTERESIS  EFFECTS  CAN  BE  INCLUDED. 
SUBROUTINE  SPRI NG ( F , DELTA , DI ST , SLOPE 1 . SL0PE2 , FORCE , NPTS  > 
DIMENSION  F(2,12) 

R°0. 

IF(DIST.GE.DELTA)GO  TO  8340 
M=2 

UT=D1ST 
GO  TO  8360 
8110  M=3 

UT=DELTA 
GO  TO  8360 
8140  F(1,K)*DELTA 
F(2,K)=F0RCE 
IF(K.GT.3)G0T08180 
GO  TO  ( 8240 ,8240, 8300 ).K 
8100  KK»K-3 
Ral. 

DO  8210  L=2,NPTS 
F(t,L)»F(l,L+KK> 

8210  F(2,L)=F(2,L+KK> 

NPTS*NPTS-KK 
GO  TO  8300 
8240  KK*3-K 
R=2. 

DO  8280  LL=1 ,NPTS 

L=NPTS+1-LL 

F(1,L+KK)*F<1,L) 

8280  F(2,L+KK)*F<2,L) 

NPTS*NPTS+KK 

8300  F(1,2)BF(1,3)-F(2,3) /SL0PE1 
F(2,2)=0. 

F ( 2 , 1 ) =-SL0PE2«F (1,2) 

F<1,1)«0. 
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22800 

22820 

22840 

22860 

22880 

22900 

22910 

22911 

22920 

22940 

22960 

22980 

23000 

23020 

23040 

23060C 

23080C 

G. 

23100 

23120 

23140 

75) 

23160 

23180 

23200 

23220 

23230 

23240 

23260 

23280 

23300 

23320 

23340 

23360 

23380 

23400 

23420 

23440 

23460 

23480 

23500 

23520 

23540 

23560 

23580 

23600 

23620 

23640 

23660 

23680 

23700 

23720 
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8340  M=1 

UT=D1ST 

8360  00  8370  J=2.NPTS 

8370  IF(F(1, J).GT.UT)G0T08380 

8380  K=J-1 

8390  FORCE® (UT-F (1,K))«(F(2,J)-F(2,K))/(F(1.J)-F(1,K))+F(2,K) 

IF(R.EQ.l. )NPTS=NPTS+KK 
IF(R.EQ.2. )NPTS=NPTS-KK 
GO  TO  <8450, 8410, 8140), M 
8410  IF < FORCE. LE.O. )G0T08450 

IF(ABS<  <F<2. J)-F(2,K) >/(F(l , J)-F ( 1 ,K) )-SLOPEl ) .LT. .01 ) 

&G0T08450 
GO  TO  8110 
8450  RETURN 
END 

THIS  SUBROUTINE  PLACES  CERTAIN  VALUES  IN  MATRIX  FORMAT  FOR  PRINTIN 

SUBROUTINE  PRINT2(X, Y.DY) 

REAL  MO 

COMHON/OUT/N,T( 175) ,X0(6, 175) ,X1 (6, 175) ,X2(6,175) ,X3<6, 175) ,X4<6,1 

&,X5(6, 175) ,X6(6, 175) ,X7(6, 175) ,X8(6,175> ,X10(6, 175) 

DOUBLE  PRECISION  Y<8) 

DIMENSION  DY<8) 

COMMON/MANDAT/ZL . ZT , ZH , RT , RN , RH . RTOPH , X2Z , Y2Z , WB , WH 
COMMON/MMANDAT/LT , LF . THFO , THLO , DELTTB, THEF 

COMMON/PARAM/RFT . THETATZ , THETAHZ , FAW , FNW , FRM , SWHC , VWHC , TAW . TRW 
COMMON/HPARAM/F ACOL , SCOLC , PG 1 , FKNEE , SF , FNCOL , MO , FRCOL , VCOLC , REACT 
COMMON/MMPARAM/BP , VOL . FFT , FP . THETAW , THETAWZ , OMW . TAC , TRC , RBRKT , OMC 
CQMMON/EABS/EAWC , EACC , E AWR , EACR , EA 

COMMQN/COLDAT/THETAC , THETACZ .MU.LSCZ, LF WZ , LBAZ , LBFZ, WC . RCOL 

COMMON/M ISC/VCOLCO , SCQLCO . PR8 , BI G 

DATA  F/ 12 . / , V/ . 68181818/ ,G/32 . 17/ , D/57 . 295780/ 

N=N+1 

IF ( N . GT . 175 ) RETURN 
CT=CQS<Y<7) ) 

T(M)=X*100O. 

X0< 1 ,N)®-DY<4)/G 
X0<2,N)»Y<4)*V 
XO  < 3 > »Y  ( 8 > »F 

X0<4,N)*<DY<1)*CT+RT#DY(3))/G 

X0(5,N)=SC0LC 

X0(6,N)=PG1 

X1Q,N)»TH£F#D 

X1(2,N)=Y ( 1 )»V 

X1(3,N)®“DY( 1)/G 

X 1 < 4 , N ) ®FKNEE*COS  < THEF ) /2 . 

XK5,N)=SF 

XH6,N)»<Y(5)-Y<8>>*F 

X2<1,N)=(Y(5)+RT#<5IN<Y<7) )-SIN( THETATZ) ) )«F 
X2(2,N)=Y(7)«D 
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23740 

23760 

23780 

23800 

23820 

23840 

23860 

23880 

23900 

23920 

23940 

23960 

23980 

24000 

24020 

24040 

24060 

24080 

24100 

24120 

24140 

24160 

24180 

24200 

24220 

24240 

24260 

24280 

24300 

24320 

24340 

24360 

24380 

24400 

24420 

24440 

24460 

24480 

24500 

24520 

24540 

24560 

24580 

24600C 

24620 

24640 

24660 

24680 

24740 

24760 
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X2(3,N)=Y(3)»D 
X2(4,N)*DY(3)»D 
X2(5,N)*X2< 1 ,N)-Y(8)»F 
X2<6.N)=><RT*Y<3>*CT+Y<1)-Y<4>)*V 

X3(1,N)*F*(Y(5)+RN*(SIN(Y(7) )-SIN(THETATZ) )+RH#<SIN<Y(6) )- 
&5IN(THETAHZ) ) ) 

X3(2.N)*Y(6)*D 

X3(3.N)=Y(2)*D 

X3(4,N)®DY(2)»D 

X3(5,N>aX3( 1 ,N)-Y(8)»F 

X3(6,N)«(Y(6)-Y(7> )#D 

X4(1,N)=FAW 

X4(2,N)»FNW 

X4<3,N)=MG 

X4(4,N)=FRW 

X4(5,N)=SWHC 

X4<6,N)aVWHC 

X5U.N)»FAC0L 

X5(2,N)*FNC0L 

X5(3.N)“H0 

X5(4,N)»FRC0L 

X5(5,N)*SC0LC 

X5(6,N)»MC0LC 

X6(1,N)«=TAW 

X6<2,N)*TRW 

X6(3,N)=THETAW*D 

X6(4,N)*0MW#D 

X7(1.N)*TAC 

X7(2,N)=»TRC 

X7(3,N)-TH£TAC*D 

X7(4,N)*0MC#D 

X8(1,N)»BP 

X8(2,N)=V0L 

X8(3,N)=PG1 

XB(4,N)=FFT 

X8(5,N)*FP 

X10( 1 ,N)*EAWC 

X10(2,N)aEACC 

X10<3,N)*£AWR 

X10(4,N)*EACR 

X10(5,N)3EA 

RETURN 

END 

SUBROUTINE  UPDATE(X, Y, DY> 

DOUBLE  PRECISION  Y(8) 

DIMENSION  DY<8> 

COMMON/KNEEREST/NKR , SKR . RUT . RELKR 
COMMON/SEATFR IC/NSF . SUN . SFUT . RELSF 
SFUT=RELSF 
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24800 

24840 

24860 


RUT*RELKR 

RETURN 

END 
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19676 
CO LUO  I P , N 


FLOW  DIAGRAM  - BAGSUB 


5 


DEFINE 
VARIABLES 
8400  - 89SjO 


DRIVER 
SUBMARINES  ' 


CALC.  SLOPE 
OF  DRIVERS' 
TORSO 

- 


CALCULATE 
NEW  "H"  PT. 
COORD.  S 
X COORD.  @ 
RIM/ BAG 
INTERSECTION 

kP80-9260 


TEST  SIGN 
OF  DISCRIM- 
INATE 
9300 

IT 


20 


CALC.  BAG 
INTERCEPT 
POINTS 


jMM 


To  19 


10 


CALC . BAG 
PENETRATION 
10460 

ZZTZZr 

CALC.  MID- 
POINT OF 
GI 


CALC.  MAJOR 
AXIS  LENGTH 
OF  ELLIPSE 
PERPENDICULAR 
TO  TORSO 
10560 


To  11 


3 


I 


1 


CALC.  PERIMETER 
OF  ELLIPSE 
PERPENDICULAR 
TO  TORSO 
10600 


f 


18 


12 


CONFINE  BAG 
VOLUME  TO 
ORIGINAL  VOL. 
IF  BAG  PRESS. 
BELOW  AMBIENT 
12660-12700 


SOLVE  FOR 
ANGLE  (BETA) 

OF  FABRIC 
TENSION  FORCE 
(NEWTON  RAPHSON 
METHOD) 
10720-11280 
1= 


19 


0-8* 

From  6 


FOR  REAL  RE- 
PEATED ROOTS 
DEFINE  XA,  XB 
12740-12780 


f 


CALC.  AREA  OF 
INTERCEPT 
11340-11540 


From  6 


CALC.  VOLUME  OF| 
BAG  INTERCEPT 


1560-11600 

T 


20_ 

FOR  COMPLEX 
ROOTS  (NO 
TORSO/ BAG 
CONTACT)  SET 
VOL=VOLO 
FTT=0 
FTH=0 

12820-12880 
21-  "~^T 


CALC.  BAG  PENET. 
FOR  SYMMETRICAL 
CASE  WITH  GIVEN 
VO  I - USE  NEWTON 
RAPHSON  METHOD 
TO  SOLVE  CUBIC 
EQUATION 

11620-12140 


From  19 


CALC.  GAS 
FLOW  INTO 
BAG  (CALL 
GAS IN) 
12900-1328D 


HK.  FOR  BOTTOMING 
OUT  OF  HEAD  6 CHEST 


HD. 


JL L 


or  so 


CALC.  RAD  & ANGLE 
OF  VERT.  BAG  ENDS 
USING  NEWTON  RAPHSONl 
METHOD  TO  SOLVE  A 
TRANSCENDENTAL  EQ. 
12160-12400 


Strike 
R-o  29 

Strike 

29 


:hk.  BAG 
DEPLOYMENT 
P G.T.  ATM 
OR 

vX  G.T.  DTI  ME 


To  21 


22 


N 

FFT=0 

FHEAD=0 

FP=0 

FTH=0 

FTT=0 

To  22 


-HD 

To  29 


CALC.  AIRBAG  VOL. 
12420-12640 

To 

? 

CALC.  NEW 
TEMP.  5 PRESS. 
DUE  TO  NET  GAS 
GAIN  IN  BAG 

-mam  - 


TORSO  IN  CONTACT 
WITH  BAG  ? 

IT 


CALC.  BAG  FORCES 
ON  BODY 
14080-14620 


I 


To  25 
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From 


From 

From 

From 

From 


24 


CALC.  DYNAMICS 
OF  STEERING 
WHEEL  CRUSH 
(CALL  WHF) 
14660-15060 

,26  I 

CALC.  DYNAMICS 
OF  STEERING 
COLUMN  CRUSH 
(CALL  COLF) 
15080-15840 
27  I 

CALC.  DYNAMICS 
OF  STEERING 
WHEEL  ROTATION 
(CALL  SPRING) 
15860-16240 


28 


T 


CALC.  DYNAMICS 
OF  STEERING 
COLUMN  ROTATION 
(CALL  SPRING) 
16260-18200 
"29  [ 


3 O— *•  UPDATE  OLD 

iso— Values 

(CALL  SPRING) 
18220-18755 


RETURN/ END 
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INPUT  VUE  NAME  7L  SRI. !>00*J> 

ENTER  1 IE  COLUMN  IS  SUPPORTED  AT  THE  FIREWALL:  OTHERWISE  ENTER  271 


INPUT  VALUES  — INPUT  UNITS<  MSEC.  MPH.  DEGREES.  INCHES.  LBS.  FT-LBS.  G'S) 


INITIAL  VELOCITY: 
INITIAL  HEAD  ANGLE: 
INITIAL  TORSO  ANGLE: 

28.0 

-6.00 

-24.0 

MLEG  MTORSO 

MHEAD 

RT 

RN 

RH 

RTOPH 

71.0  58.4 

11.4 

14.0 

20.5 

4.75 

27.7 

NPTS  NECK  NPTS  KR 

NPTS  VEH 

NPTS  SEAT 

NPTS  GAS 

NPTS  COL 

NPTS  WHC 

4 7 

NP  WRP  NP  WRN 

6 6 

SL  SEAT  SL  KR 
0.500E+04  0.500E+04 
GAS  FLOW  TIME  - MSEC 

12 

NP  CRP 
4 

6 

NP  CRN 
4 

12 

12 

11 

0.  17.5 

0 51.0  56.0 

GAS  FLOW  TIME  - MSEC 
71.0  100. 

GAS  FLOW  - LB/SEC 

20.5 

26.0 

31.0 

36.0 

41.0 

46. 

0.  0. 

3 0.980  0.490 

GAS  FLOW  - LB/SEC 

0.  0. 

COLUMN  STROKE  - INCHES 

3.40 

4.09 

4.26 

3.60 

2.45 

1.6 

-10.0  0. 

7 3.05  3.85 

COLUMN  STROKE  - INCHES 
4.13  5.70 

COLUMN  FORCE  - LBS 

0.550 

0.950 

1.55 

1.60 

2.34 

2.5 

0.  0. 

680.  860. 

COLUMN  FORCE  - LBS 
0. 110E+04  0. 151E+04 

710. 

900. 

730. 

500. 

500. 

930 

SEAT  FRICTION  DISPLACEMENT  - INCHES 


-50.0 

0. 

1.00 

14.0 

15.0 

50.0 

SEAT 

FRICTION  FORCE  - LBS 

0. 

0. 

250. 

250. 

0. 

0. 

NECK 

ANGLE 

- DEG 

-80.0 

18.0 

20.5 

90.0 

NECK 

TORQUE 

- FT-LBS 

115. 

0. 

■6.70 

-86.7 

VEH. 

PULSE 

TIME  - MSEC 

0. 

4.00 

14.0 

25.0 

33.0 

45.0 

0 

73.0 

102. 

VEH. 

PULSE 

TIME  - MSEC 

118. 

200. 

VEH. 

PULSE 

DECELERATION  - 

G'S 

0. 

1.00 

10.0 

11.2 

8.90 

13.0 

0 

17.2 

8.60 

VEH. 

PULSE 

DECELERATION  - 

G'S 

-0.500 

-0.500 

KNEE 

DISPLACEMENT  - INCHES 

-10.0 

3.15 

3.25 

3.58 

3.68 

4.13 

KNEE  FORCE  - LBS 


0. 

0.  400. 

900. 

0.210E+04 

0.250E+04 

0.250E+04 

WHEEL  STROKE 

- INCHES 

-10.0 

0.  0.230 

0.500 

0.910 

1.72 

1.95 

6 3.00 

3.66 

WHEEL  STROKE 

- INCHES 

6.00 

WHEEL  FORCE  - 

LBS 

0. 

0.  0. 108E+04 

0. 131E+04 

0. 117E+04 

945. 

990. 

. 0. 112E+04 

680. 

WHEEL  FORCE  - 

LBS 

0. 136E+04 

WHEEL  ANGLE. 

POS.  - DEG 

0. 

1.60  4.90 

8.90 

16.0 

23.3 

WHEEL  TORQUE. 

POS.  - IN-LBS 

0.  0. 

399E+04  0.380E+04 

0.340E+04 

0.303E+04 

0.285E+04 

WHEEL  ANGLE. 

NEG.  - DEG 

0. 

1.60  4.90 

8.90 

16.0 

23.3 

WHEEL  TORQUE. 

NEG.  - IN-LBS 

0.  0. 

399E+04  0.380E+04 

0.340E+04 

0.303E+04 

0.285E+04 

COLUMN  ANGLE. 

POS.  - DEG 

0. 

5.00  6.00 

7.00 

COLUMN  REACTION.  POS.  - IN-LBS 

IF  2:  LBS 

IF  1 

0. 

50.0  0.200E+04 

0.200E+04 

COLUMN  ANGLE. 

NEG.  - DEG 

0. 

5.00  6.00 

7.00 

COLUMN  REACTION.  NEG.  - IN-LBS 

IF  2:  LBS 

IF  1 

0. 

5?.0  0.200E+04 

0.200E+04 

ATMOP 

PGZ  GTZ 

U 

PN1 

PN2 

PN3 

14.7 

-14.7  0. 166E+04 

662. 

1.40 

1.40 

1.40 

VC1 

VC2  AO 

SA 

SC 

X1Z 

Y1Z 

0.700 

0.700  2.07 

12.5 

8.00 

30.6 

18.9 

THETAC 

MU  LSCZ 

LFWZ 

LBAZ 

LBFZ 

WC 

15.2 

0.  11.0 

16.5 

15.0 

19.0 

7.00 

LF 

THFO  THLO 

18.0 

12.0  35.0 

XWH 

RIMRAD  X2Z 

Y2Z 

UB 

WH 

-6.00 

7.38  34.5 

7.40 

11.2 

5.00 

THETAW 

RIP  WWH 

WIP 

RCOL 

DCN 

15.2 

2.50  3.00 

2.00 

15.0 

0.690 

CHEST  IMPACT  WITH  LOWER  WHEEL  RIM.  RUN  STOPPED. 


E-2 


TIME. 

(MS) 


0. 

5.00 

10.00 

15.00 

20.00 

25.00 

30.00 

35.00 

40.00 

45.00 

50.00 

55.00 

60.00 

65.00 

70.00 

75.00 

80.00 

85.00 

90.00 

95.00 

100.00 


TIME 

(MS) 


0. 

5.00 

10.00 

15.00 

20.00 

25.00 

30.00 

35.00 

40.00 

45.00 

50.00 

55.00 

60.00 

65.00 

70.00 

75.00 

80.00 

85.00 

90.00 

95.00 

100.00 


OFH  G'S 

OEH  0E.L 

VEH  DISP 

(G'S) 

(MPH) 

(INCHES) 

zzmotss 

S3SC.X93 

sea  at  a a % 

0. 

28.00 

0. 

1.90 

27.93 

2.46 

6.40 

27.47 

4.90 

10.11 

26.53 

7.28 

10.65 

25.39 

9.57 

11.20 

24.19 

11.75 

9.76 

23.04 

13.83 

9.58 

22.02 

15.81 

11.29 

20.88 

17.70 

13.00 

19.55 

19.48 

16.75 

17.92 

21.13 

20.50 

15.87 

22.62 

19.90 

13.56 

23.91 

16.40 

11.57 

25.02 

16.10 

9.86 

25.96 

16.61 

8.02 

26.75 

15.12 

6.28 

27.37 

13.64 

4.70 

27.86 

12.16 

3.29 

28.21 

10.68 

2.04 

28.44 

9.19 

0.95 

28.57 

FEM  ANGLE 

H-P  OEL 

H-P  ACC 

(DEG) 

(MPH) 

(G'S) 

12.00 

28.00 

0. 

12.00 

28.03 

-0.31 

12.00 

28.06 

-0.23 

12.00 

28.08 

0.03 

12.00 

28.05 

0.57 

12.00 

27.94 

1.40 

12.00 

27.73 

2.50 

13.47 

27.44 

2.66 

15.10 

27.01 

4.40 

16.81 

26.46 

5.26 

18.62 

25.83 

6.11 

20.55 

25.12 

6.61 

22.59 

23.39 

30.37 

24.36 

19.75 

33.51 

25.68 

16.15 

32.17 

26.64 

12.65 

31.74 

27.28 

9.17 

31.80 

27.61 

5.66 

32.29 

27.59 

2.12 

30.67 

27.24 

-0.82 

22.52 

26.67 

-2.73 

11.11 

E-3 


BODY  G'S 

COL  DISP 

BAG  PRESS 

(G'S) 

(INCHES) 

(PSIG) 

2s.sic3ac 

cesxtarzsz 

ssczsas 

-0.35 

0. 

-14.70 

-0.46 

0. 

-14.70 

-0.47 

0. 

-14.70 

-0.49 

0. 

-14.70 

-0.52 

0. 

-13.56 

-0.57 

0. 

-9.40 

-0.63 

0. 

-4.66 

-0.62 

0. 

-0.01 

-10.55 

0.08 

3.22 

-16.30 

0.70 

5.12 

-18.69 

1.05 

6.18 

-17.95 

1.72 

6.45 

-17.59 

2.75 

6.07 

-14.58 

2.90 

5.39 

-11.53 

2.92 

4.57 

-10.04 

2.92 

4.09 

-9.43 

2.92 

3.82 

-9.15 

2.92 

3.58 

-8.95 

2.92 

3.38 

-8.40 

2.92 

3.22 

-7.39 

2.92 

3.08 

FEM  FORCE 

SEAT  FR. 

H-P  R.D. 

(LBS) 

(LBS) 

<I~6-0 

0. 

0. 

0. 

0. 

0.83 

0.00 

0. 

7.62 

0.03 

0. 

30.37 

0.12 

0. 

76.59 

0.31 

0. 

147.11 

0.59 

0. 

240.50 

0.96 

0. 

250.00 

1.41 

0. 

250.00 

1.92 

0. 

250.00 

2.49 

0. 

250.00 

3.14 

0. 

250.00 

3.89 

1002.77 

250.00 

4.75 

1138.74 

250.00 

5.55 

1126.51 

250.00 

6.18 

1117.30 

250.00 

6.66 

1110.95 

250.00 

6.99 

1107.69 

250.00 

7.17 

1063.03 

110.58 

7.16 

773.13 

0. 

6.97 

299.76 

0. 

6.68 

G 


nm 

TORSO  DISP 

TORSO  ANG 

TORSO  OEL 

TORSO  ACC 

TORSO  R.D. 

TORSO  R.V 

(MS) 

(INCHES) 

(DEG) 

(D/SEC) 

<D/SEC««2> 

(INCHES) 

(MPH) 

0. 

0.00 

-24.00 

0. 

-550.61 

0.00 

0. 

5.00 

2.46 

-24.01 

-5.64 

-1172.29 

0.00 

0.03 

10.00 

4.92 

-24.06 

-11.33 

-1071.77 

0.02 

0.45 

15.00 

7.38 

-24.13 

-16.00 

-723.00 

0.09 

1.35 

20.00 

9.83 

-24.21 

-18.07 

0.52 

0.26 

2.43 

25.00 

12.27 

-24.30 

-15.59 

1113.88 

0.52 

3.55 

30.00 

14.71 

-24.36 

-6.65 

2595.12 

0.88 

4.61 

35.00 

17.14 

-24.36 

7.42 

2855.39 

1.33 

5.52 

40.00 

19.54 

-24.33 

-7.84 

-10339.37 

1.84 

6.03 

45.00 

21.85 

-24.53 

-81.44 

-18186.32 

2.37 

5.89 

50.00 

24.01 

-25.18 

-180.12 

-20804.07 

2.88 

5.65 

55.00 

25.99 

-26.34 

-282.25 

-19005.27 

3.37 

5.74 

60.00 

27.81 

-27.89 

-311.86 

14628.15 

3.89 

6.00 

65.00 

29.42 

-29.21 

-210.99 

23172.70 

4.41 

5.62 

70.00 

30.84 

-29.97 

-88.82 

25819.48 

4.88 

5.22 

75.00 

32.09 

-30.08 

44.77 

27518.07 

5.34 

5.17 

80.00 

33.17 

-29.51 

185.39 

28831.09 

5.79 

5.13 

85.00 

34.10 

-28.22 

333.44 

30507.13 

6.24 

5.04 

90.00 

34.89 

-26.16 

488.40 

29359.21 

6.68 

4.92 

95.00 

35.55 

-23.39 

611.67 

19383.94 

7.11 

4.94 

100.00 

36.12 

-20.14 

675.87 

4808.77 

7.55 

5.14 

TIME 

HEAD  DISP 

HEAD  ANG 

HEAD  VEL 

HEAD  ACC 

HEAD  R.D. 

HEAD  R.AN 

(MS) 

(INCHES) 

(DEG) 

(D/SEC) 

(D/SEC»»2) 

(INCHES) 

(DEG) 

0. 

0.00 

-6.00 

0. 

-486.74 

0.00 

18.00 

5.00 

2.46 

-5.97 

12.94 

2919.14 

0.00 

18.04 

10.00 

4.93 

-5.87 

27.58 

2909.81 

0.02 

18.19 

15.00 

7.39 

-5.70 

42.11 

2797.49 

0.11 

18.43 

20.00 

9.85 

-5.45 

55.78 

2470.10 

0.28 

18.76 

25.00 

12.31 

-5.14 

67.43 

1910.03 

0.56 

19.16 

30.00 

14.77 

-4.78 

75.94 

1136.57 

0.95 

19.58 

35.00 

17.23 

-4.39 

78.87 

410.72 

1.42 

19.96 

40.00 

19.68 

-3.89 

165.98 

45885.45 

1.98 

20.44 

45.00 

22.10 

-2.34 

475.21 

75055.84 

2.62 

22.19 

50.00 

24.46 

0.99 

859.39 

75370.16 

3.33 

26.17 

55.00 

26.76 

6.16 

1192.31 

54718.84 

4.14 

32.50 

60.00 

28.96 

12.78 

1448.55 

43370.10 

5.05 

40.67 

65.00 

31.05 

20.41 

1584.59 

14783.20 

6.03 

49.63 

70.00 

33.00 

28.45 

1614.89 

-2245.84 

7.04 

58.42 

75.00 

34.79 

36.45 

1575.94 

-12969.93 

8.04 

66.53 

80.00 

36.41 

44.13 

1487.05 

-23031.97 

9.04 

73.63 

85.00 

37.87 

51.23 

1342.49 

-35251.10 

10.01 

79.44 

90.00 

39.16 

57.44 

1126.18 

-55018.74 

10.96 

83.60 

95.00 

40.31 

62.25 

775.59 

-84585.90 

11.87 

85.64 

100.00 

41.34 

65.40 

548.65 

-28691.19 

12.77 

85.55 

E-4 


TIME 

WH  AX  FOR 

WH  N f OR 

WH  MOMENT 

WH  RESIST 

WH  STROKE 

wh  bt  m 

(MS) 

(LBS) 

(LBS) 

(IN  LBS) 

(LBS) 

( INCHES) 

(IN/SEC) 

carrsas 

CXQBSCa 

znn  r,  to* 

cnanr 

ac  tots  name 

c ttr.naar. 

ia*c 

0. 

0. 

0. 

0. 

0. 

0. 

“99 

0 

0.  o. 

0. 

0. 

0. 

0. 

10.00 

0. 

0. 

0. 

0. 

0. 

0. 

15.00 

0. 

0. 

0. 

0. 

0. 

0. 

20.00 

0. 

0. 

0. 

0. 

0. 

0. 

25.00 

0. 

0. 

0. 

0. 

0. 

0. 

30.00 

0. 

0. 

0. 

0. 

0. 

0. 

35.00 

0. 

0. 

0. 

0. 

0. 

0. 

40.00 

614.84 

474.29 

1396.32 

424.36 

0.09 

47.04 

45.00 

955.71 

763.68 

2411.73 

654.96 

0.14 

59.19 

50.00 

1144.86 

935.04 

3088.70 

1101.88 

0.26 

0. 

55.00 

1183.34 

998.75 

3532.51 

1101.88 

0.26 

0. 

60.00 

1052.38 

982.22 

3718.27 

1101.88 

0.26 

0. 

65.00 

789.50 

871.21 

2976.67 

1101.88 

0.26 

0. 

70.00 

572.68 

717.67 

2211.17 

1101.88 

0.26 

0. 

75.00 

495.47 

624.31 

1832.77 

1101.88 

0.26 

0. 

80.00 

471.79 

559.83 

1592.05 

1101.88 

0.26 

0. 

85.00 

456.71 

499.70 

1370.11 

1101.88 

0.26 

0. 

90.00 

450.36 

445.02 

1170.11 

1101.88 

0.26 

0. 

95.00 

452.08 

395.47 

989.05 

1101.88 

0.26 

0. 

100.00 

458.94 

349.86 

823.90 

1101.88 

0.26 

0. 

TIME 

COL  AX  FOR 

COL  N FOR 

COL  MOMENT 

COL  RESIST 

COL  STROKE 

COL  ST  OEL 

(MS) 

(LBS) 

(LBS) 

(IN-LBS) 

(LBS) 

(INCHES) 

(IN/SEC) 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

5.00 

0. 

0. 

0. 

0. 

0. 

0. 

10.00 

0. 

0. 

0. 

0. 

0. 

0. 

15.00 

0. 

0. 

0. 

0. 

0. 

0. 

20.00 

0. 

0. 

0. 

0. 

0. 

0. 

25.00 

0. 

0. 

0. 

0. 

0. 

0. 

30.00 

0. 

0. 

0. 

0. 

0. 

0. 

35.00 

0. 

0. 

0. 

0. 

0. 

0. 

40.00 

693.16 

492.39 

1396.32 

107.26 

0.08 

61.60 

45.00 

1053.01 

775.57 

2411.73 

783.39 

0.70 

97.57 

50.00 

1272.79 

947.45 

3088.70 

870.91 

1.05 

66.33 

55.00 

1337.61 

1017.44 

3532.51 

500.00 

1.72 

211.40 

60.00 

1201.62 

1004.48 

3718.27 

838.23 

2.75 

117.29 

65.00 

910.68 

899.17 

2976.67 

760.12 

2.90 

0. 

70.00 

682.81 

759.29 

2211.17 

746.56 

2.92 

0. 

75.00 

607.83 

669.81 

1832.77 

746.56 

2.92 

0. 

80.00 

574.00 

600.98 

1592.05 

746.56 

2.92 

0. 

85.00 

548.82 

536.39 

1370.11 

746.56 

2.92 

0. 

90.00 

532.45 

477.12 

1170.11 

746.56 

2.92 

0. 

95.00 

524.21 

422.87 

989.05 

746.56 

2.92 

0. 

100.00 

521.18 

372.49 

823.90 

746.56 

2.92 

0. 

E-5 


TIME 

WH  AP  MOM 

WH  RES  MOM 

WH  ANGLE 

WH  ANG  VEL 

(MS) 

(IN-LBS) 

(IN-LBS) 

(DEG) 

(DEG/SEC) 

sse  = = s==  trsssss  r:sssar  assssss  :sss==s 


0. 

0. 

0. 

15.25 

0. 

5.00 

0. 

0. 

15.25 

0. 

10.00 

0. 

0. 

15.25 

0. 

15.00 

0. 

0. 

15.25 

0. 

20.00 

0. 

0. 

15.25 

0. 

25.00 

0. 

0. 

15.25 

0. 

30.00 

0. 

0. 

15.25 

0. 

35.00 

0. 

0. 

15.25 

0. 

40.00 

1396.32 

-623.99 

15.50 

0. 

45.00 

2411.73 

-1791.85 

15.97 

114.67 

50.00 

3088.70 

-2285.58 

16.17 

0. 

55.00 

3532.51 

-2285.58 

16.17 

0. 

60.00 

3718.27 

-2285.58 

17.27 

0. 

65.00 

2976.67 

-2285.58 

21.68 

0. 

70.00 

2211.17 

-1666.36 

26.07 

0. 

75.00 

1832.77 

-1979.04 

27.60 

0. 

80.00 

1592.05 

-1979.04 

27.60 

0. 

85.00 

1370.11 

-1979.04 

27.60 

0. 

90.00 

1170.11 

-1979.04 

27.60 

0. 

95.00 

989.05 

-1979.04 

27.60 

0. 

100.00 

823.90 

-1979.04 

27.60 

0. 

TIME 

COL  AP  MOM 

COL  RES  MOM 

COL  ANGLE 

COL  ANG  OEL 

(MS) 

(IN-LBS) 

(IN-LBS) 

(DEG) 

(DEG/SEC) 

0. 

0. 

0. 

15.25 

0. 

5.00 

0. 

0. 

15.25 

0. 

10.00 

0. 

0. 

15.25 

0. 

15.00 

0. 

0. 

15.25 

0. 

20.00 

0. 

0. 

15.25 

0. 

25.00 

0. 

0. 

15.25 

0. 

30.00 

0. 

0. 

15.25 

0. 

35.00 

0. 

0. 

15.25 

0. 

40.00 

0. 

0. 

15.25 

0. 

45.00 

0. 

0. 

15.25 

0. 

50.00 

0. 

0. 

15.25 

0. 

55.00 

0. 

0. 

15.25 

0. 

60.00 

28583.04 

-151.60 

16.35 

551.17 

65.00 

25099.61 

-14306.89 

20.76 

1198.79 

70,00 

20872.90 

-27155.57 

25.40 

607.16 

75.00 

18295.23 

-27155.57 

26.80 

0. 

80.00 

16362.79 

-27155.57 

26.80 

0. 

85.00 

14553.34 

-27155.57 

26.80 

0. 

90.00 

12896.64 

-27155.57 

26.80 

0. 

95.00 

11382.22 

-27155.57 

26.80 

0. 

100.00 

9978.92 

-27155.57 

26.80 

0. 

E-6 


TIME 

BAG  PEN. 

BAG  00L . 

BAG  PRESS. 

W/A  FORCE 

P.  FORCE 

(MS) 

(INCHES) 

(CU.1N. ) 

(PSIG) 

(LBS) 

(LBS) 

0. 

6.44 

4955.72 

-14.70 

0. 

0. 

5.00 

6.45 

4955.69 

-14.70 

0. 

0. 

10.00 

6.46 

4953.91 

-14.70 

0. 

0. 

15.00 

6.54 

4946.32 

-14.70 

0. 

0. 

20.00 

6.70 

4929.66 

-13.56 

0. 

0. 

25.00 

6.95 

4903.04 

-9.40 

0. 

0. 

30.00 

7.29 

4866.14 

-4.66 

0. 

0. 

35.00 

7.71 

4819.76 

-0.01 

0. 

0. 

40.00 

8.05 

4829.80 

3.22 

75.11 

628.19 

45.00 

8.00 

4833.03 

5.12 

122.66 

1006.16 

50.00 

8.11 

4821.96 

6.18 

152.74 

1211.82 

55.00 

8.07 

4828.26 

6.45 

165.42 

1270.15 

60.00 

7.65 

4874.70 

6.07 

158.37 

1197.56 

65.00 

7.45 

4949.09 

5.39 

126.26 

971.58 

70.00 

7.12 

5027.59 

4.57 

94.13 

748.49 

75.00 

7.27 

5033.32 

4.09 

82.02 

641.33 

80.00 

7.69 

5002.82 

3.82 

75.63 

590.68 

85.00 

8.13 

4967.16 

3.58 

69.71 

549.73 

90.00 

8.57 

4925.98 

3.38 

64.92 

518.07 

95.00 

8.99 

4879.53 

3.22 

61.23 

494.40 

100.00 

9.40 

4828.82 

3.08 

58.19 

476.10 

TIME 

CHEST  AP 

CHEST  SI 

HEAD  AP 

HEAD  SI 

(MS) 

(G’S) 

(G’S) 

(G’S) 

(G’S) 

0. 

-0.35 

0. 

-0.59 

0.16 

5.00 

-0.46 

0.12 

-0.10 

0.37 

10.00 

-0.47 

0.09 

-0.09 

0.34 

15.00 

-0.49 

-0.01 

-0.07 

0.22 

20.00 

-0.52 

-0.23 

-0.04 

-0.04 

25.00 

-0.57 

-0.57 

-0.01 

-0.44 

30.00 

-0.63 

-1.03 

0.02 

-0.99 

35.00 

-0.62 

-1.10 

-0.08 

-1.08 

40.00 

-10.55 

-1.81 

-3.51 

3.15 

45.00 

-16.30 

-2.11 

-4.71 

7.10 

50.00 

-18.69 

-2.24 

-6.99 

11.85 

55.00 

-17.95 

-2.05 

-8.98 

16.59 

60.00 

-17.59 

-13.13 

-9.00 

6.94 

65.00 

-14.58 

-15.86 

-13.77 

5.20 

70.00 

-11.53 

-15.98 

-16.12 

4.78 

75.00 

-10.04 

-15.88 

-18.13 

4.78 

80.00 

-9.43 

-15.29 

-19.71 

4.95 

85.00 

-9.15 

-14.04 

-20.85 

4.47 

90.00 

-8.95 

-10.89 

-21.46 

3.99 

95.00 

-8.40 

-4.81 

-21.25 

4.73 

100.00 

-7.39 

1.22 

-3.07 

7.36 

TIME 

EAUC 

EACC 

EAUR 

EACR 

EA 

<HS) 

(LBS) 

(LBS) 

(LBS) 

(LBS) 

(LBS) 

sssssss 

SSS3BSS 

SS33SCSS 

sesstsss 

szssesacs 

xssssss 

0. 

0. 

0. 

0. 

0. 

0. 

5.00 

0. 

0. 

0. 

0. 

0. 

10.00 

0. 

0. 

0. 

0. 

0. 

15.00 

0. 

0. 

0. 

0. 

0. 

20.00 

0. 

0. 

0. 

0. 

0. 

25.00 

0. 

0. 

0. 

0. 

0. 

30.00 

0. 

0. 

0. 

0. 

0. 

35.00 

0. 

0. 

0. 

0. 

0. 

40.00 

440.38 

111.17 

0. 

1000000.00 

0. 

45.00 

681.25 

811.98 

0. 

1000000.00 

0. 

50.00 

1147.24 

902.70 

0. 

1000000.00 

0. 

55.00 
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